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INTRODUCTION 


The free-energy minimization method of references 1 and 2 for determining thermo- 
chemical equilibrium properties of arbitrary gases and gas mixtures may be employed as 
a subroutine in a variety of real-gas computations. One example of such usage is refer- 
ence 3 where the equilibrium normal- shock and shock-tube solutions for carbon dioxide 
are presented. Unless the inputs to the free-energy minimization method are pressure 
and temperature, an iterative procedure is necessary which requires considerably more 
computer time. For real-air computations (ref. 4), this problem was avoided by employ- 
ing a magnetic tape containing thermodynamic properties of air over a wide range of 
conditions and by using appropriate tape search routines to interpolate for the desired 
thermodynamic properties. A similar tape has been generated for the thermodynamic 
properties of carbon dioxide. The corresponding species mole fractions are stored on 
two additional tapes. This report presents the data recorded on these tapes in the form 
of tables and charts. These data are for temperatures from 100 K to 25 000 K and pres- 
sures from 40 mN/m^ to 1 GN/m2. 

Unpublished tables for thermodynamic properties of pure carbon dioxide were gen- 
erated by Perry A. Newman of the Langley Research Center a decade ago. However, 
these unpublished tables are more restrictive in range of conditions than the present 
tables and do not include mole fractions. Charts for equilibrium properties of carbon 
dioxide are included in reference 5 for temperatures to 25 000 K. In the present study, 
the 10 species used in reference 5 (e - , C, C + , C ++ , O, 0+, 0 ++ , O 2 , CO, and CO 2 ) were 
supplemented by 6 additional species (o~, O^, O 2 , C", C 2 , and CO + ). Hence, the pres- 
ent results are believed to provide more accurate and detailed mole fractions for carbon 
dioxide than reference 5. Also, provision of tables as well as charts is a worthwhile con- 
venience, particularly in regions where curves on the charts overlap and for the 16 mole 
fractions. 


SYMBOLS 

a speed of sound, m/s 

Cp specific heat at constant pressure, kJ/kg-K 

Cy specific heat at constant volume, kJ/kg-K 


h 


specific enthalpy, J/kg 



p 

R 

s 

T 

W 

w 0 

Z 

>'E 

P 


pressure, N/m 2 

universal gas constant, 8.31434 kJ/kmole-K 
specific entropy, kJ/kg-K 
temperature, K 
molecular weight, kg/kmole 

molecular weight of undissociated carbon dioxide, 44.011 kg/kmole 

kmole of dissociated carbon dioxide per kmole of undissociated carbon 
dioxide, W Q /W 

isentropic exponent, Q ^ 

density, kg/m 2 


CONVERSION FACTORS AND CONSTANTS 

Conversion factors between the International System of Units (SI) and U.S. Customary 
Units (ref. 6) for the quantities presented in the present tables and charts are as follows; 

1 N/m 2 = 9.8692 x 10' 6 atm = 1.4504 x 10" 4 psi = 2.0885 x 10“ 2 lbf/ft 2 

1 kg/m 2 = 6.2428 X 10“ 2 lbm/ft 2 = 1.9403 x 10 -2 slug/ft 2 

1 J/kg = 1 m 2 /s 2 = 10.764 ft 2 /s 2 = 4.3021 x 10" 4 Btu/lbm 

1 m/s - 3.2808 ft/s = 2.2369 mph 

Physical constants appearing herein are 

p 0 = 101.325 kN/m 2 

W 0 = 44.011 kg/kmole 
aQ = 259.15 m/s 
p Q = 1.9626 kg/m 2 
T 0 = 273.15 K 
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STRUCTURE MEDIUM INTERACTION 


t W8 0 -17294 i 

FAILURE OF UNDERGROUND CONCRETE STRUCTURES 

SUBJECTED TO BLAST LOADINGS* j 


C. A. Ross 

University of Florida Graduate Engineering Center 
Eglin AFB, Florida 32542 

and 

P. T. Nash and G. R. Griner 
USAF Armament Laboratory 
Eglin AFB, Florida 32542 


This study presents the results of analytical predictions 
of response and failure of two-edges -free reinforced concrete 
slabs subjected to intermediate blast loadings. Approximate 
theoretical methods using stationary and moving plastic hinge 
mechanisms with linearly varying and time dependent loadings 
are developed. Reasonable agreement between analysis and ex- 
perimental observations was obtained. 


INTRODUCTION 

Failure of reinforced concrete 
structures, as well as other structural 
elements subjected to distributed normal 
loads, may be defined in terms of the 
degree of damage. In this study failure 
of reinforced concrete structures will 
be defined as a condition where actual 
separation or fracture of the reinforc- 
ing elements has occurred. This essen- 
tially means the structure will have 
been rendered almost completely incapa- 
ble of any load carrying ability. To 
reach this damage level the reinforcing 
elements at the failed position will 
have passed through the yield and ulti- 
mate stresses. For this type condition 
a reasonable assumption is that the ten- 
sile strength of the concrete is negli- 
ble when compared to the strength of the 
reinforcing elements. Failure in this 
c ( as/s is then based on the ultimate 
strength or elongation of the reinforc- 
ing elements. For plates and beams with 
transverse loads where bending is pres- 
ent, an idealized rigid perfectly plas- 
tic stress-strain assumption may be used 
to estimate the true moment carrying 

*7his study funded by the Joint Techni- 
cal Coordinating Group/Surface Target 
Survivability Vulnerability Program 
through USAF office of Scientific Re- 
search Grant 78-3592 in cooperation with 
the USAF Armament Laboratory, Eglin AFB, 
Florida, 


capacity of the structure. This assump- 
tion is justified for ductile steel re- 
inforcing elements whose true stress- 
strain relation shows a rather small 
elastic portion when compared to the 
overall stress-strain diagram. The 
"fully plastic" or "ultimate" moment as 
applied to statically loaded reinforced 
concrete beams and plates is summarized 
by Szilard [1]. The "fully plastic mo- 
ment" and "plastic hinge" for metal 
beams is covered in detail by Timoshenko 
and Gere [2] and plastic hinges and yield 
lines as applied to reinforced concrete 
beams are shown in some detail in [1] . 
Abrahamson, et al [3] used both station- 
ary and moving plastic hinges to describe 
the deformation process of thin metal 
beams subjected to various kinds of 
transverse dynamic pressure loadings. 
Jones [4] used this simple rigid-plastic 
method to successfully predict the large 
inelastic behavior of thin metal beams 
when subjected to impulsive loads and 
analytically verified the experimental 
work of Menkes and Opat [5]. 

Experimental evidence of Fuehrer 
and Xeeser [6] shows that buried rein- 
forced concrete beams exhibit a plastic 
hinge type deformation when subjected to 
a small underground explosion of close 
proximity. The results of this experi- 
mental program [6! show two kinds of 
failures II a team failure, with a 
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plastic hinge on each end and in the 
middle, with no appreciable concrete 
cracking elsewhere, associated with low- 
ered pressure loads, and 2) failure of 
reinforcing elements at sides and middle 
with large areas of cracked concrete 
over the interior section of the beam 
showing evidence of moving hinges asso- 
ciated with loadings larger than those 
of number 1 above. The failure associ- 
ated with complete shearing at the ends 
was not evidenced but should not be dis- 
counted as the loadings may not have ap- 
proached those necessary for complete 
end shearing. Based on the experimen- 
tally determined failure modes described 
above, a study was initiated to analyti- 
cally describe these type failure modes 
and predict reinforced concrete beam de- 
flections and failures. The approach to 
be used was chosen as the previously de- 
scribed plastic hinge methods and the 
remainder of this paper is concerned 
with development and results obtained by 
this analytical approach. 


ANALYSIS 

Based on a thod used by Abrahamson 
[3] , using two failure mechanisms as 
shown in Figure 1 for fixed end beams 
with a linear load distribution of Fig- 
ure 2, equations have been developed to 
predict deflection and failure of rein- 
forced concrete beams. The beam loading 
shown in Figure 2 is assumed linear spa- 
tially and every point on the beam has 
the same time variation given as 


f(t) = (l-t/T)exp(-at/x) (1) 


where t is time, a is a dimensionless 
decay constant and t is the length of 
the pressure'pulse in seconds. The load- 
ing is expressed in terms of a uniform 
loading denoted by P*. at the edges and 
a uniformly varying triangular load of 
magnitude P c at the center and zero at 
the edges. These loads are assumed con- 
stant across the beam width b. Using 
Figure 2 the beam loading, in force per 
area of the beam, at a distance x from 
the beam end is given as 

P(x,t) = (P E +P c |)f(t) 

for 0 sxsa and 0 stii (2) 

P (x,t) ** 0 for all x and t>T 



Figure la. Mechanism 1, showing station- 
ary hinges at ends and midspan M u «hinge 
moment per unit length. 
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Figure lb. Mechanism 2, showing station- 
ary hinges at ends and two traveling 
hinges. 
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t>T 



0<t<T 


(3) 


where f(t) is given by Equation (1). 

The equation of motion for mecha- 
nism 1 is given by 


where 0 is the angular rotation in radi- 
ans at the hinge moment as shown in Fig- 
ure 3, M u is the plastic hinge moment/ 
unit length, m is the beam mass per unit 
area, w is the weight per unit area, F 
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is edge fixity term, and n determines 
loading direction. F equals 1.0 for 
simply supported edges and 2.0 for fixed 
edges. The load direction indicator n 
equals +1, -1, 0 for beamn with loads 
from above, from below, and vertical 
beams respectively. The center point 
deflections and velocity are given as 


The equation of motion for the flat 
center portion is given by 

l - IT (p E + 7^ )_ 5r' 0<H - a ' °* tiT ' 

(5) 

'& ■ 0, t>T , 


da. 


A 


6a 


(4) where H is the hinge location. 


where a is the beam half span. 



area. 


Using the assumption of zero shear 
at the hinge the equation of motion for 
rotation of the outer portion of length 
H has been determined to be 



Equations (6) reduces to Equation (3) 
when H reaches a and a beam deforming 
initially in a mechanism 2 mode will re- 
duce to a mechanism 1 mode in the final 
stages of the deflection process. Equa- 
tions (5) and (6) are coupled through the 
continuity condition that 

A ■ H6 


and 


(7) 


A = H0 



The initial hinge location Hg may be de- 
termined by the initial condition that 


(A = H 0 0) t = Q (8) 

where H- is the initial hinge location 
for t = u 0. Applying the iritil condi- 
tion of Equation (8) results in an equa- 
tion from which the initial hinge loca- 
tion is determined knowing the applied 
loads and the hinge moment M u . The equa' 
tion for H Q is given as 


u 4 fj 

-§-[PE +p C<i^ "D-nwl * 3FM u . (9) 


Figure 3. Beam notation and deflection 
for mechanism 1. 


The static collapse load Pv S (at edge of 
beam) for a mechanism 1 mode is that load 
required to produce a bending moment M u 
at the ends or midpoint of the beam. For 
the loading assumed with K - Pq/P e , 


P ES 


12M„ 
a 2 (2K+3) 


(10) 
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Figure 4. Beam notation and deflection 
for mechanism 1. 



Figure 5. Beam cross section used for 
calculating M u . 


Solving Equation (10) for M„ and substi- 
tuting into the right hand side of equa- 
tion (9) results in 


y 2 (l+K(Y-l)-nw] - ■• 2 - K ?\ )F (ID 


where Y « Hg/a and X is the ratio of the 
applied edge load P E to the center load 
P E s required for static collapse. The 
solution for K - 0 results in Y - 1 for 
X ■ 3 which agrees with the results of 
Abrahams on , et al [3] for the uniformly 
loaded fixed end vertical beam. The 


solution of equation (12) then dictates 
the mechanism and type of equations of 
motion to be used. For 0<Y<a equations 
(5), (6) and (7) are used and for Y - 1 
equations (3) and (4) are used. Solu- 
tions where Y>1 the mechanism 1 solution 
is used. As X-*», Y-*0, therefore very 
large values of X are not r.xlowed in 
these solutions. 


The plastic or hinge moment is given 
by Szilard [1] as 


( 12 ) 


M 0 - 0.9 (d 2 qa R (l-0.59qor/a c )] 


where d is distance between reinforcing 
elements (Figure 5) , q is the reinforcing 
ratio in tension, o R is the yield strass 
of the reinforcing material and o c is 
the maximum compression strength of the 
concrete. Equation (12) was used as a 
basis for all calculations in the follow- 
ing sections. 

In this study the failure criteria 
was chosen as steel fracture and the lo- 
calized plastic strain to fracture ap- 
pears to be more eaeily estimated than 
the stress at failure. Since angular 
rotation 6 is the foremost independent 
parameter, strain in the reinforcing 
element will be expressed in terms of 
the rotation. In the usual manner 
strain will be presented in two parts, 
i.e. et strain due to change in length 
and Cb strain due to rotation. The de- 
formed length £ at the elastic hinges 
for homogeneous metal 1 ams [4] was as- 
sumed to be some multip .e or fraction 
of the thickness. In this case the de- 
formed length of the reinforcing elements 
at the plastic hinge will ba chosen as 
the distance between the elements of the 
reinforcing mesh. This assumption seems 
justified on the basis that, although 
the fr%e distance of -ui element in ten- 
sion is imbedded in concrete, contrac- 
tion due to Poisson effect will occur 
and allow elongation of the element be- 
tween the two transverse adjacent par- 
allel reinforcing elements fixed in the 
concrete. Neglecting higher order terms, 
change in the half span length may be ap- 
proximated for a mechanism 1 mode as 

r~~— 2 

6 * /A 2 +a 2 -a - a(l+i£,+ • • •-!) (13) 

2a 2 

. A 2 

' Ji 


where a is the halfspan length and A is 
the midspan deflection. If half of this 
change takes place at the end support 
and the remainder at the midspan, then 
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I - R ■ 


(16) 


6 is given in terms of e fc and i as 



(14) 


and the strain due to bending becomes 


The rotation 0 may be expressed 
0»A/a and the strain e t becomes 



(15) 


The strain e b due to bending at the 
plastic hinge may be expressed as d/2R 
where R is the radius of curvature and 
d is the distance between the reinforc- 
ing mats as shown in Figure 6 . The ra- 
dius of curvature is defined by 


e b * 0d/2i (17) 

The total strain e may be written as 

e = e t +e b * (18) 

A similiar equation may be determined 
for the mechanism 2 case where H the 
current hinge position is used in place 
of the midspan length a. 


REINFORCING 



Figure 6. Schematic of slab thickness showing notation, arrangement of reinforcing 
elements and the assumed hinge length. 
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EXPERIMENTAL TESTS 

Tests {61 were performed on several 
sizes of reinforced slabs witn two edges 
fixed and selected specimens from these 
tests are presented here for comparison 
with analytical results. 

The specimens were fabricated using 
6000 psi compression concrete reinforced 
with a steel wire 4 in x 4 in (10.16 x 
10.16 cm) mesh both top and bottom. Ad- 
ditional reinforcements were placed at 
the ends of the specimens as seen in the 
fractured specimen of Figure 13. The 
general schematic of test specimens is 
shewn in Figure 7 and a test specimen in 
place is shown in the photograph of Fig- 
ure 6. The legs of the specimen were 
placed against an upright concrete wall 
and were prevented from parallel outward 
motion at the wall by steel angles as 
shown in Figure 8. The explosive charge 
was placed on a line normal to and at 
various distances from the center of the 
slab. Pressure wave speeds were meas- 
ured at various distances from the 
charge . 

All tests were performed with the 
test specimen vertical. The specimens 
were then covered and packed with soil 
to a height approximately 3 ft (.91m) 
above the top of the specimen. This ex- 
perimental configuration eliminated the 
weight of overburden. The top of the 
upright specimen was covered to prevent 
backfill behind the test slab. 


RESULTS AND DISCUSSION 

The deflection or rotation of the 
slab may be determined using equations 
(3) or (5) and (6) , however the key to 
the initial response is found using 
equation ( 1 1 > . The initial position of 
the plastic hinge Y may be found using 
equation (11) and this position is nec- 
essary if a solution of deflection ver- 
sus time is required. 



FigureB. Photograph showing topsoil re- 
moved with test specimen in place. 



Figure 7. Schematic of test setup with soil removed. Test arrangement is comoletelv 
covered with soil before tests. * 


6 



Figure 10. Post-test photograph of 
Case 1. 


Figure 11. Post-test photograph of 
Case 2. 


However a simple Heck on the expected 
response may be determined, for example, 
on a fixed end beam, l -y letting Y-l, 

F^2, n=0, and solving for \=2K+3, which 
in the boundary between the mechanism 1 
md mechanism 2 initial response. There- 
fore for V < (2K+ J) the initial response 
is by mechanism L and for \ > (2K+3) the 
initial response is by mechanism 2. The 
use of equation (11) in the determination 
of the initial response mechanism re- 
quires the maximum pressure of the 


pressure loading, as well as the pres- 
sure distribution over the slab. The 
general shape of the pulse is necessary 
for a complete solution however general 
relations for certain cases similiar to 
those of reference |3| may be written 
using only the impulse and peak pressure 
of the pressure loading. The pressure 
loadings used were estimate from the 
measured pressure time curves of refer- 
ence |6). Well-defined pressure loadings 
are desired due to the requirement for 
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initial peak loading and total impulse , 
however, only estimates were available. 

Equations (1) through (12) have 
been coded and solutions of the three 
cases given below In Table I were pro- 
grammed. The analytical predicted re- 
sponse of Case 1 is shown in a scaled 
drawing of Figure 9 and the expez omen- 
tal final position and failure of this 
same slid) is shown in Figure 10. Ana- 
lytical prediction of Case 2 gave a max- 
imum midspan deflection of approximately 
1.0 inch (2.54 cm). Experimentally, one 
test of Case 2 showed about 1.5 inch 
(3.81 cm) midspan deflection, however, 
a second test, shown in Figure 11, indi- 
cates no deflection and minimal cracking. 
The analysis of Case 3 gives a mechanism 
2 initial response and the response- time 
relation is shown in the scaled drawing 
of Figure 12. The final experimental 
deflection for Case 3 is shown in Figure 
13. 

Using the failure criteria as pre- 
sented in the Analysis Section, the ro- 
tation for failure 8 m may be determined 
if the strain to failure is known. 
Post failure observations show that the 
distance between reinforcing elements is 
a reasonable estimate of the deformed 
hinge width. 

Fractures of the reinforcing ‘ele- 
ments occur either at the welds of the 


mesh or at localized points between the 
welds. Post test measurements of rein- 
forcing elements show a very localized 
reduction in area near the fracture 
point with very little or no plastic de- 
formation at positions away from the 
fracture point. Laboratory measurements 
on the mesh verified the post test meas- 
urements and deflection load curves in- 
dicate a yield/ultimate of 60,000-70,000 
psi (414 - 483 MPa) and a strain to fail- 
ure of 0.10 to 0.15. 

Using an experimentally observed 
strain to failure of 0.125 and equation 
(19) gives a rotation to failure of 0.22 
radians for a 28 in (.71 m) slab initial- 
ly deforming in a mechanism 1 mode. Ex- 
perimental evidence shows that for the 
28 in (.71 m) long slabs some reinforc- 
ing steel was fractured for all rotations 
above 0.20 radians. 

Mathematically, the distinction be- 
tween mechanism 1 and mechanism 2 is a 
line, however, for a real structure, 
where the hinge may only have to fall 
within half the hinge length of the mid- 
span to be classified as a mechanism 1 
mode. For Case 1, the predicted initial 
hinge location was found to be only a- 
bout 10% of the half span from the cen- 
ter and the response mode is distinctly 
mechanism 1 as shown in Figure 10. How- 
ever, for the Case 3 the free span is so 
short it is difficult to determine what 


TABLE I 


REINFORCED CONCRETE SLAB TEST, 2% STEEL, 
FIXED TWO EDGES. 



Effective Size 

Explosive Size 

Explosive 

PC 

PE 

T 


CASE 

ft (m) 

lbs (kg) 

Distance ft(m) 

psi (MP a ) 

psi (MP a ) 

Sec 

a 


3X2.33X.33 

8 

3 

140 

420 

.015 

6 


(.91X.71X.1) 

(3.63) 

(.91) 

(.97) 

(2.90) 




3X2.33X.33 

8 

4 

50 

190 

.015 

2 


(.91X.71X.1) 

(3.63) 

(1.22) 

(.34) 

(1.31) 



1 

3X1.33X.33 

8 

2 

380 

1520 

.015 

30 


(.91X.41X.1) 

(3.63) 

(.61) 

(2.62) 

(10.48) 


i 
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OPTIMIZATION OF REINFORCED CONCRETE SLABS 


J. M. Ferrittc 

Civil Engineering Laboratory 
Naval Construction Battalion Center 
Port Hueneme, California 


The Department of Defense uses reinforced concrete cells composed of 
concrete slabs to limit the effects of accidental explosions during 
hazardous explosives operations. An automated design procedure is 
discussed which considers the dynamic nouliuear behavior of the rein- 
forced concrete of arbitrary geometrical and structural configuration 
subjected to dynamic pressure loading. Optimum design of the slab is 
accomplished by use of an interior penalty function. The paper dis- 
cusses the optimization procedure and a discussion of the results is 
given. The results are compared with finite element analysis. 


INTRODUCTION 

The Department of Defense has numerous 
facilities engaged in the production of various 
types of explosives and munitions used by mili- 
tary services. In most cases the production of 
aammnition utilizes assembly line procedures. 
Projectiles pass through various stages of 
preparation: filling with explosive, fuzing, 

marking, and packing. Hazardous operations, 
such as the filling of the projectile case with 
an explosive in a powder form and the compaction 
of the powder by hydraulic press, are accom- 
plished in protective cells that are intended to 
confine the effects of an accidental explosion. 

Most of the existing production facilities 
were built in the 1940s. With few exceptions, 
the manufacturing technology and existing equip- 
ment represent the state of the art as of 1940. 


The production equipment was operated extensive- 
ly during World War II, again during the Korean 
conflict, and recently during the Southeast Asia 
war. Much of this equipment and the housing 
structures have been operating beyond their 
designed capacities |1). The Department of 
Defense is conducting an ammunition plant mod- 
ernization program (2) that is intended to 
greatly enhance safety in the production plants 
by protective construction, automated pro- 
cessing, and reduction of personnel involved in 
hazard operations. An automated procedure was 
required to give structural designers the capa- 
bility to perform rapid analysis of the struc- 
tural safety of blast-resistant construction. 

The design parameters interact in a complex way 
since the procedure is both nonlinear and dynam- 
ic. From a design point of view an optimization 
procedure was required to minimize cost and 
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maximize safety since blast-resistant construc- 
tion has been reported to cost 3 to 5 times as 
auch as conventional construction. Therefore, 
the first objective was to autoaate the analysis 
procedures for determining structural response 
of reinforced concrete slabs having a bilinear 
stiffness representation and subjected to blast 
shock and gas pressures. Concrete slabs are the 
basic eleaent forming sidewalls, roofs and 
floors of cells designed to confine the effects 
of accidental explosions. The second objective 
was to provide an optimum design procedure for 
laced (shear reinforced) and unlaced reinforced 
concrete slabs that will automatically produce a 
least-cost design for a given slab geometry, 
material properties, and explosive weight for 
both feasible and nonfeasible starting points. 

Analysis Procedure 


the design constraints; or, stated in optimiza- 
tion terms: Find such that M(?C) is a minimus 
aud 


g.(£) SO i = 1, 2, ... N 


where $ 
N 

8 

M 


vector of design variables 
number of design constraints 
vector of design constraints 
objective function 


Specifically for this problem the design vari- 
ables selected are areas of steel reinforcement 
and thickness of concrete. The design con- 
straints are minimum section properties and the 
flexural and shear limits. The objective func- 
tion consists of the costs of formwork, con- 
crete, and flexural and shear (lacing) rein- 
forcement. 


The procedure treats reinforced rectangular 
slabs with arbitrary boundary conditions, cross 
section properties, material properties and 
loading. The determination of the equivalent 
pressure load based on the quantity of explosive 
and cell geometry is based on Reference (3J. 

The computational procedure for the slab analy- 
sis utilizes reinforced concrete section proper- 
ties to determine resisting moment, from which 
using yield line analysis techniques, the static 
resistance of the slab is determined. The 
stiffness of the slab is calculated using elas- 
tic plate theory. Using the loading, resis- 
tance, stiffness and mass of the slab, the 
dynamic response of the slab is computed by time 
step iteration. The slab is modeled as an 
equivalent single degree-of-f reedom elasto- 
plastic system [4]. The design of the slab, 
based on its ultimate strength allows large 
nonlinear deformations to confine the effects of 
an accidental explosion. The deflections are 
limited in terms of support rotations. 

Structural Optimization 


Fixed Variables 

W = explosive weight 
H = wall height 
EL = wall length 

h = height of explosive above floor 

£ = distance of explosive from left 
side of wall 

R g - distance of explosive from wall 
I = reflection code 

f . = ultimate dynamic concrete strength 

dc 

f , = dynamic yield strength of reinforc- 

' in? steel 

6 = rotation criterion 
Design Parameters, X 


t 

c 

concrete thickness 

AVT 

area of vertical reinfccing 
steel tension 

AVC 

area of vertical reinforcing 
steel compression 

AHT 

area of horizontal reinforc- 
ing steel tension 

AHC 

area of horizontal reinforc- 
ing steel compression 

, AS 

area lacing steel 


The optimization problem consists of find- 
ing the least-cost structure that satisfies all 


12 




Constraints. g (X) 


6(X) = 6(6), maximum deflection 

V(X) S VC for 6 i 2 deg, maximum shear 
t 2 12, minimum thickness 

minimum steel rein- 
forcement tension or 
compression 

The methodology [5,6] selected used the 
unconstrained minimization approach. The prob- 
lem is o. verted to an unconstrained minimiza- 
tion by constructing a function, $, of the 
general form 


AV i 0.0025 bd I 
AH i 0.0025 bd I 


♦(X, r) = M(X) + P[ gl (X), .... g n (X), r] 

For this problem the interior penalty function 
technique was selected. This methodology is 
suitable when gradients are not available, and, 
because the method uses the feasible region 
(region where all constraints are satisfied), a 
useable solution always results. The objective 
function is augmented with a penalty term that 
is small at points away from the constraints in 
the feasible region, but increases rapidly as 
the constraints are approached. The form is as 
follows: 


♦ (X, r) = M(X) + r 2 —^r 

1 = 1 gj(X) 


+ (AV + AH) (EL • H)C $ 
+ (A s )(EL • H)C l 


where C 

c 

= 

cost of concrete ($/cu ft) 

C 

s 

s 

cost of horizontal and vertical 
reinforcement ($/cu in.) 

C L 

= 

cost of lacing reinforcement 
($/ cu in.) 

A 

s 


area lacing reinforcement 
($/ cu in.) 

AV 

— 

vertical steel (tension and 
compression) 

AH 

= 

horizontal steel (tension and 
compression) 



? 1 


where r = penalty parameter. 


The program requires a starting point in 
the feasible region before optimizat' in can 
proceed. This is accomplished automatically by 
the program by incrementing the design variables 
until a feasible point is reached. 

An algorithm which comprises the steps most 
commonly used is as follows: 

1. Given a starting point, X Q , satisfying 

all gj(X) > 0 and an initial value for 

r, minimize A to obtain X . 

min 

2. Check for convergence of X . to the 

min 

optimum. 


where M is to be minimized over all £ satisfying 
g (S) S 0, j = 1 ... N. Note that if r is 
positive, then, since at any interior point all 
of the terms in the sum are negative, the effect 
is to add a positive penalty to M(X). As the 
boundary is approached, some g (X) will approach 
zero, and the penalty will increase rapidly. 

The parameter, r, will be made successively 
smaller in order to obtain the constrained 
minimum of M. 

Objective function 

Cost = M = H ■ EL • t • C 
c c 


3. If the convergence criterion is not 
satisfied, reduce r by r *■ rc, where 

c < 1. 

4. Compute a new starting point for the 
minimization, initialize the minimiza- 
tion algorithm, and repeat from step 1. 

The logic diagram for the interior penalty 
functions technique is shown in Figure 1. 

The minimization for $(X, r) shown in 
Figure 1 is accomplished by a method developed 
by Powell using conjugate directions [5,6]. 

Powell's method can be summarized as 
follows: Given that the function has been 
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Figure 1. Logic diagraa for interior penalty 
function technique. 


ainiaized once in each of the coordinate direc- 
tions and then in the associated pattern direc- 
tion, discard one of the coordinate directions 
in favor of the pattern direction for inclusion 
in the next ainiaizations , since this is likely 
to be a better direction than the discarded 
coordinate direction. After the next cycle of 
ainiaizations, generate a new pattern direction, 
and again replace one of the coordinate direc- 
tions. 

Figure 2 is a logic diagraa for the uncon- 
strained ainiaization algorithm. The pattern 
aove is constructed in block A, then used for a 
ainiaization step (blocks B and C) , and then 

stored in S (block D) as all of the directions 
n 

are up-nuabered and Sj is discarded. The direc- 
tions S Q will then be used for a ainiaizing step 
just before the construction of the next pattern 
direction. Consequently, in the second cycle, 
both X and Y in block A are points that are 
ainiaa along S Q , the last pattern direction. 

This sequence will iapart special properties to 
S n+ j = X - Y that are the source of the rapid 
convergence of the aethod. 


Figure 2 shows a block requiring a one- 
diaensional ainiaization of a* of the function 
4($ + a S^) . The one -diaensional ainiaization 
uaes a four-point cubic interpolation. It 
finds the ainiaua along the direction S^, where 
$ is the coordinate of the previous ainiaua. By 
trial and error it finds three points with the 
aiddle one lets than the other two. It aakes a 
quadratic interpolation, and then a cubic inter- 
polation. If the actual function evaluated at 
the new interpolated point is not sufficiently 
close to that of the preceding point or if it is 
not sufficiently close to the interpolated 
function, then another cubic interpolation is 
SMde. 



Figure 2. Logic diagraa for minimization 
of *(X). 
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Discussion of Results 


•how* the constraints and objective function 
considering for this example that reinforcement 
The objective function is segaentally in tension and compression is the same and that 

linearly dependent on the design variables; vertical reinforcement is related to horixontal 

however, the constraints are both linearly and reinforcement (thus reduced variables for a 

nonlinearly related to the design variables. 2-dimensional presentation). The shear stress 

The aininuai area of steel aud ainiaua concrete and deflection are nonlinearly related to the 

thickness are linear constraints. Figure 3 thickness of concrete. The shear stress is 



Art* of lloritonttl Reinforcing Steel, AH (sq in 7ft) 


Figure 3. Design space, 2 sides fixed. 
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further noted not to be convex. This will be 
discussed later. Figure 3 shows the useable 
region bounded by flexure, shear, and minimum 
steel constraints. The optimum least-cost 
solution is shown. This specific example solu- 
tion considers an unlaced section; thus, the 


maximum shear constraint is active. Laced 
sections eliminate the shear constraint. If the 
number of sides supported were increased from 2 
to 3, the design space would change as shown in 
Figure 4. There are two regions that are use- 
able areas. Obviously, the lower one offers the 



Area of Horizontal Reinforcing Steel, AH (sq in./ft) 


Figure 4. Design space, 3 sides fixed. 
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There it clearly a complex interaction of 
conatrainta ahowing the effect of ahear. Unfor- 
tunately, the optimal aolution found by the 
prograa depends on the starting point selected. 
The prograa converges on the cloaeat relative 
optima. Several alternative starting points 
should be used to verify a questionable optiaua. 
Revising the design paraaetera could possibly 
shift the constraints such that only one useable 
solution would appear. However, a slight in- 
crease in shear stress (10%) can significantly 
reduce coat by allowing the near-optiaua non- 
feasible solution to be accepted, Figure 5. 

The dual-space problea of finding a useable 
aolution is peculiar to unlaced concrete slabs 
only because lacing (shear reinforceaent) aodi- 
fiea the shear constraint. Nonautoaated optiaua 
design for unlaced conditions is alaiost impossi- 
ble when one considers thj complexity of the 
design space and the large nuaber of iterations 
required when an initial solution is not feasi- 
ble. 

Cost data used in the program can be se- 
lected by the user. However, the data herein is 
based on work by Picatinny Arsenal on contract 
with Aassann and Whitney (7). Table 1 shows a 
coapariaon of unlaced and laced concrete walls 
for various boundary conditions . The example 
considers a 15-foot-high by 12-foot-wide (4.6 a 
x 3.7 a) wall subjected to a 200-psi (1,400 kPa) 
10 -as triangular loading function. In all cases 
the laced concrete (12-degree rotation) ia less 
expensive than unlaced (2-degree rotation) 
designs. For the unlaced sections ultimate 
deflections expressed as support rotation capa- 
city auat be limited to 2 degrees whereas laced 
sections are allowed to reach 12 degrees by 
virtue of the added concrete confinement pro- 
vided by the lacing. 

COMPARISON OF RESULTS TO FINITE ELEMENT 
TECHNIQUES 

To evaluate the accuracy of the approxima- 
tions awde in the analytical technique described 


ibove, a coapariaon was aade to a finite element 
prograa, IKSLAB {8]. Prograa INSLAB ia a finite 
element prograa having plate element with 12 
degree: freed' An elasto-plastic material 
model was uaef throughout this work. The model 
was defined by tho elastic mdulus, yield modu- 
lus, equivalent yield stress assuming rectangu- 
lar section (yield moment divided by thickness 
squared) and Poisson's ratio. Yielding is 
determined by principal moments. Uniform pres- 
sure loading is input and mass is lumped at the 
modes based on material density and integrating 
the mass over the element interpolation func- 
tion. The solution procedure uses time step 
integration asauming a constant acceleration 
between successive time increments. 

A reinforced concrete slab supported on 
three sides was analyzed for a 400 psi (2,800 
kPa) pressure load with 1.17 msec duration. The 
slab measured 4 ft by 4 ft with 6-inch thickness 
(1.2 m x 1.2 m x 15 cm). Horizontal and verti- 
cal reinforcement was 0.18 sq in. /ft (3.8 sq 
cm/m) both sides. A comparison of peak dis- 
placement from both the finite element and 
approximate solutions is shown in Figure 6. 

The finite elesient analysis was performed 
using the computer program INSLAB, [8], The 
mesh consisted of 80 elements each having 12 
degreea of freedom. A time increment of 0.01 ms 
was used. The displacement difference at maxi- 
mum value is about 13 percent. 

Both techniques used an elasto-plastic 
analysis. The approximate solution uses a 
cracked moment of inertia averaged with a gross 
moment of inertia to determine stiffness. The 
finite element method uses a gross moment of 
inertia based on section thickness. However, 
the modulus of the concrete was reduced to 
account for reduction in stiffness resulting 
from cracking. The rise time of the load for 
the finite element solution was increased to 
allow a number of steps before peak load. This 
effect is minimal since the loading is impulsive 
in nature and ia thought to result in about a 5 
percent reduct. on in peak diaplacement. 
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uKtblc area 



Area of Horizontal Reinforcing Steel, AH (»q in./ft) 


Figure S. Reviaed deaign spice, 3 sides fixed 



Deflection (in.) 



Figure 7 shows the deflection contours 
obtained by the finite .. lement analysis and the 
yield line obtained froa the approximate solu- 
tion. The agreement is very good. The yield 
line marks regions of the plate assumed to 
displace as units. The contour lines tend to 


run parallel to the support and ver* teal yield 
line. This is the same mechanism predicted by 
the yield pattern since the slab deflects most 
along the vertical yield line varying linearly 
to the nondeflecting support. 


free 



Figure 7. Deflection contour of 4 ft x 4 ft (1.2 m x 1.2 m) slab. 


20 




Table 1. Coaparison of Optiaua Solutions 


(200 psi; 10 ms; wall, 12 ft L x 15 ft H) 
(1,400 kPa; 10 ms; wall 3.7 m Lx 4.6 m H) 



N 

Side 


Shear 

Reinforcement/ 

Deflection 

Cost 

($) 

I 

N=2 

7777777/ 


unlaced 2 degrees 
laced 12 degrees 

3,290 

2,289 

i 

h 

N=3 

\ 

unlaced 2 degrees 

2,753‘ 

7/777773 

V//////i 

y 

laced 12 degrees 

2,019 


N=4 


unlaced 2 degrees 

2,001 U 

J 

7T7TT77) 


laced 12 degrees 

1,958 


'Shear capacity slightly exceeded. 


CONCLUSION 

An optiaua design procedure is given fcr 
reinforced concrete slabs using the internal 
penalty function approach. The solution has 
been found to converge rapidly and at ainiaua 
coaputer cost. Results agree well with more 
refined eleaert analysis. This prograa is in 
use by the Navy in the design of facilities to 
resist dynaaic pressures, such as froa acciden- 
tal explosions. 
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A NUMERICAL COMPARISON WITH AN EXACT SOLUTION FOR THE TRANSIENT 
RESPONSE OF A CYLINDER IMMERSED IN A FLUID 


M. E. Giltrud, D. S. Lucas 
Naval Surface Weapons Center 
White Oak, Silver Spring, Maryland 20910 


The transient response of an elastic cylindrical shell 
immetsed in an acoustic media that is ingulfed by a plane 
wave is determined numerically. The ntchod applies to the 
USA-STAG5 code which utilizes the finite element method 
for the structural analysis and the Doubly Asymptotic 
Approximation (DAA) for the fluid-structure interaction. 
The calculations are compared to an exact analysis for two 
separate loading cases: a plane step wave and an 

exponentially decaying plane wave. The results of the 
comparisons are most favorable which supports the 
j applicability of USA-STAGS for fluid-structure interaction 
. problems . 


INTRODUCTION 

Fluid structure interaction 
problems have recently become quite 
important to the Navy. Fundamental 
knowledge of the fluid-structure 
interaction problem is a basis for a 
more thorough understanding of damage 
to submerged structures from underwater 
shock loading. The difficulty in 
developing analytical models to predict 
the response of submerged structures to 
underwater explosions is twofold. 

First, the structure must be suitably 
modeled in order to determine its 
elastic-plastic response. This 
certainly is within the capability of 
many structural analysis computer 
codes. Secondly, the loading effects 
are altered according to the state of 
motion of the structure in the fluid. 
Thus, the problem in determining the 
loading is one in which the state of 
motion of the fluid and structure are 
coupled. 

Analytical solutions heve been 
obtained for simple elastic structural 
geometries with simple temporal and 
spatial waveforms. However, in some 
cases the results have to be generated 
numerically since the resulting 
solut.on is quite complex. In addition, 
analytical solutions seldom address 
the problem of internal structure in 


any detail. Analytical investigations 
have resulted in the adoption of 
methods for uncoupling the structure 
from the fluid. 

These procedures, known as 
surface approximation techniques, have 
been compared by Geers [1]. It is 
apparent from his analysis that the 
Doubly Asymptotic Approximation (DAA) 
[1-3] is the most accurate for early 
and late times in predicting the 
behavior of submerged shells. The DAA 
also affords a smooth transition 
between early and late time response. 
More recently, the DAA has been used 
in conjunction with several linear 
elastic finite element computer codes 
to predict the response of submerged 
targets to underwater shock loading 
[4, 51. 

Since the close-in underwater 
detonation of high explosives causes 
large local plastic deformation, a 
structural analysis computer code 
capable of performing such analyses 
for general shell structures is 
required. The STAGS [6, 7] (Structural 
Analysis of General Shells) code has 
been selected since it is highly 
efficient in performing such analyses. 

The STAGS code has been combined 
with a DAA code, USA (Underwater Shock 
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Analysis) , with the resulting code 
called USA-STAGS. In the following 
sections we discuss the application of 
USA-STAGS to the problem concerning 
the response of a linear elastic 
cylindrical shell immersed in a fluid 
that is engulfed by a plane acoustic 
wave. The results of the USA-STAGS 
calculations are compared to an exact 
analysis by Huang [8]. 

EXACT ANALYSIS METHODOLOGY 

Huang's analysis focuses upon the 
problem of the transient interaction of 
a plane acoustic wave with a circular 
cylindrical shell immersed in a fluid. 
For the circular cylindrical shell, 
Huang states that "this paper presents 
an exact formulation by simple 
transformations of the basic 
differential equations governing the 
coefficients of the series solution for 
the shell deflections and the wave 
pressure for each "mode* into a 
Vol terra integral equation of the 
second Kind." The integral equations 
are then solved by a step-by-step 
integration scheme. This method is 
also very efficient for parametric 
studies, and the entire computation 
process is quite straight forward, and 
can easily be controlled to have a high 
degree of accuracy. 

Exact series solutions are presented 
in this work for the cases of a 
neutrally bouyant shell submerged in 
water subjected to step and 
exponentially decaying incident waves. 
The complete pictures of the total 
motion and stress responses of the 
shell are revealed. The material 
properties used were those of steel for 
the shell and water for the fluid. 

The case which has been chosen to model 
on USA-STAGS is that for a thickness 
to radius ratio of 1/31, which 
corresponds to the dimensionless 
parameter M = 2. [8] 

In the derivation of the shell 
equations, the relative changes in 
length and shear of the shell middle 
surface and the change of curvature and 
twist have been accounted for as in 
Junger [9], Thus, the shell equations 
U3ed are quite general and are exact 
for the elastic behavior of the shell 
within the Kirchoff assumption of shell 
theory [10]. The geometry of the 
problem for the circular cylindrical 
shell and the non-dimensionalized 
velocity responses are shown in Figures 
1 and 2 . 


USA-STAGS ANALYSIS PROCEDURE 

The USA-STAGS code makes use of 
the finite element method for the 
structure and the DAA for the fluid- 
structure interaction. The coupled 
system of equations are 


M x + K x = - T A (P e + P T ) (1) 

% — X 'V X — o —l 


A P_ + p e A M_ _1 A P = 

p c A (T T x - x T ) 

x x — 


Variables appearing here are defined 
as: 

M Structural mass matrix 

K Linear or nonlinear structural 

% stiffness operator 

x Structural displacement vector 

x Acceleration vector 

T A transformation matrix 

^ (rectangular) which relates 

displacement variables in 
fluid to those ir> structure 

A Fluid element area matrix 

^ (diagonal) 

Fg Vector representing scattered 

fluid pressure 

Pj Vector representing incident 

fluid pressure 

Mp Fluid mass matrix 

p Fluid density 

c Speed of sound in fluid 

Xj Vector of fluid particle 

acceleration due to the 
incident wave 

Step by step integration of these 
equations carried out with a staggered 
solution procedure described below: 

(1) estimate x at next step by 

extrapolation of previous 
value 


(2) 

solve 

(2) 

for 

5s 

(3) 

solve 

(1) 

for 

X 
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This solution procedure is condition- 
ally stable. Therefore, the fluid 
equation (2) has been altered to 
achieve unconditional stability [11]. 

In order to predict the transient 
response of a cylindrical shell 
immersed in a fluid excited by an 
acoustic wave, a finite element model 
of the cylinder has been constructed. 
Since the plane strain response of the 
cylinder is desired, a finite element 
model of arbitrary length could be 
used. In fact, a finite element model 
consisting of one element along the 
length is adequate. However, in order 
to model the added mass effects 
(described by Mp in equation (2) ) a 

longer model is required. The actual 
model used for the analyses has a 
length to diameter ratio of 10. 

RESULTS AND DISCUSSION 

Two analyses have been performed 
with USA-STAGS and the results have 
been compared to the exact analyses by 
Huang. The first analyses concerns a 
step pressure wave that is shown in 
Figure 3a and the second concerns an 
exponentially decaying pressure wave 
that is shown in Figure 3b. The 
standoff distance from front of the 
cylindrical shell to the point source 
of 100 diameters is used. This 
effectively reduced the sphericity of 
the wave to nearly a plane wave. 

The integration procedure used for 
the staggered solution is the Park's 
linear multi-step method [12]. Since 
the integration procedure for the 
coupled equations is unconditionally 
stable, a time step is chosen to 
achieve accuracy in the solution. A 
time step of about 1/50 of a transit 
time has been used for all analyses. 
(The transit time is the time required 
for the acoustic wave to traverse the 
shell . In this case it is 

1.67 x 10 4 sec.) 

The result of the USA-STAGS 
calculations are compared to the eract 
analyses in Figures 4 and 5. Figure 4 
gives the velocity response of the 
shell at the horizontal centerline for 
the step wave loading. The calcula- 
tions compare very well with the exact 
analyses for the 1.75 transit time 
shown. The USA-STAGS results show good 
agreement for both frequency and 
amplitude. 

The velocity response comparisons 
for the exponential wave are shown in 
Figure 5. USA-STAGS results tend to 
oscillate more than the exact analysis. 


This can be accounted for by the 
observation that the exact solution 
contain contributions from the first 
seven modes. The USA-STAGS calcu- 
lations contain contributions from 
higher modes. These higher modes are 
excited by the rapidly decaying 
incident pressure wave (the duration 
of the loading is 0.06 transit time). 

SUMMARY 

A comparison of an exact solution 
for the velocity response of a 
cylindrical shell immersed in a fluid 
with that of a solution generated by 
a fluid analyzer and structural code 
(USA-STAGS), for plane shock wave 
loading has been presented. Good 
agreement has been obtained between 
the two results for step and 
exponential pressure wave loading. 

These comparisons are supported by 
previous investigations. 
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CASE STUDIES IN DYNAMICS 


FOIL SYSTEM FATIGUE LOAD ENVIRONMENTS 
FOR COMMERCIAL HYDROFOIL OPERATION 


D, L, Graves 
Boeing Marine Systems 
Renton, Washington 


Complexity and severity of the hydrofoil fatigue loads environment in 
the open sea Is without parallel In terms of flight vehicle environments. 
The random nature of wave orbital velocities, periods and heights plus 
boat heading, speed and control system design must all be considered to 
assess the structural fatigue requirements. Major nonlinear load events 
such as hull slannlng and foil unwetting are key contributors to the 
fatigue environment. Full scale "rough water" load tests, field experi- 
ence plus analytical loads work on the Model 929 JETFOIL commercial hydro- 
foil provided a major Increase In knowledge of the hydrofoil loads environ- 
ment, The problem of developing an overall sea environment for design Is 
discussed. Structural fatigue damage Is shown to be dominated by the most 
severe and Infrequent load environments which are also the most difficult 
to define statistically. State of the art analytical approaches were 
found to be limited In providing satisfactory predictions of the complex 
fatigue load environments, 


Nomenclature 


Symbol Definition 



a 

L AA 

L AS 

l a 

X 



Exceedance frequency for load level X 

Average wave encounter frequency 
measured from moving craft 

Significant wave height, equal to the 
average of 1/3 highest waves 

Total dynamic lift load on front foil 
semi span 

Nominal steady state lift on front 
foil semi span 

Maximum expected lift on front foil 
semi span 

Incremental load factor, (Lp-Lp )/Lp 

s* s * m 

Total dynamic antisymmetric lift load 
on aft foil 

Total dynamic Symmetric lift load on 
aft foil 

Nominal steady state ltft on aft foil 

Maximum expected lift on aft foil 


Incremental antisymmetric foil load 
factor, (L w -L a }/L^ 

Incremental symmetric foil load factor, 

(L «- \.l\ 


N Fatigue life, cycles 

P(X^>X) Probability that peak value X^ exceeds 


Symool Definition 

T , Average stationary wave period of 1/3 

9 highest waves, seconds 
V w Orbital wave velocity 

o Standard deviation 

u Circular frequency, rad/sec 

I. Introduction 


The general problem of developing design 
criteria for a random loading environment Is a 
familiar one. However, In the specific case 
of a hydrofoil operating In the open sea, the 
problem has received minimal theoretical deve- 
lopment, In terms of practical design experi- 
ence, the state of the art of hydrofoil fatigue 
environments has been nonexistent. Even though 
displacement hull ships have been a cannon form 
of transportation for centuries, It Is only 
within the last several years that the problem 
of determining displacement hu.l statistical 
behavior In a random sea has received serious 
attention In the technical literature. For the 
hydrofoil craft, the need for accurately 
assessing fatigue load environments Is more 
urgent. Weight and performance are both criti- 
cal and the loading environments more severe 
than for displacement hulls. For these design 
considerations the hydrofoil Is closer to a 
fighter airplane than a displacement hull ship. 
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It Is accurate to categorize the fatigue load 
environments of a modern sea going hydrofoil 
craft with submerged foils as being more 
severe and complex than that of most, If not 
all, other flight vehicles. At follborne 
speeds the sea Imposes a continuous spectrum 
of random gusts on the foils and struts. If 
an active control system which operates p 
surfaces is used, It Is continuously Inter- 
acting with the environment to control and 
stabilize hull motions. The wetted area on 
the struts Is In a state of d. ..tmlc transition. 
Cavitation Is an Important hydrodynamic 
phenomena which affects the load behavior. In 
large seas, major nonlinear load events occur 
when the foils approach or break the water 
surface or the hull slams Into a wave. Not 
only the magnitude but the frequency of these 
load events must be predicted to provide 
adequate load c-lterla for the structural 
design. 

The Model 929-100 JETFOIL Is currently the 
only commercial hydrofoil In operation which 
employs submerged foils and operates on a 
dally basis In open seas at various locations 
throughout the world. Operational experience 
to date has Included routes In Hong Kong, the 
Pacific Ocean off Hawaii, Caribbean waters off 
Venezuela, the Sea of Japan and the English 
Channel. Full scale "rough water" load tests, 
coimterclal field experience and analytical 
simulation studies of the JETFOIL have pro- 
vided a major Increase In understanding and 
appreciation of the hydrofoil fatigue ''oad 
environment. 


II. Vehicle Configuration and Design 


Figure 1 Illustrates the basic JETFOIL hydro- 
foil craft. The hull Is of aluminum construc- 
tion co achieve low weight while the foils and 
stmts are made of a high strength steel alloy. 
The JETFOIL has an upside down "T-Tall" for 
the front foil system. The aft foil system 
consists of port and starboard foil segments 
joined together by a large center strut and 
"capped" on each end by an outboard strut. 
Propulsion is accomplished by two waterjet 
pumps which receive water from a duct Inlet . 
on the aft center strut. Cruise speed while \ 
follborne Is over 40 Knots. An automatic ( 
control system essentially flies the JETFOIL, 


FIGURE 1 MODEL 929' 


Activation of full span flaps on the foils Is 

accomplished by the control logic to stabilize 

lull motions and maintain a nominal foil depth 

(6.5 feet). Onboard accelerometers and gyros 

provide hull motion Inputs to the control logic. 1 

A sonic height sensor on the bow Is ustd to , 

sense the water surface and thus provide a 

means of controlling the depth of the submerged f 

foils, A steerable fror* trut Is used for turn 
coordination and also reacts to lateral acceler- ! 

at Ion feedback commands. Roll stability Is 1 

maintained by differential flap deflections on ] 

the aft foil. | 

Two control modes are used: a platform mode and 1 

a contour mode. In the platform mode the i 

control system has strong acceleration feedback | 

gains and attempts to nullify boat accelerations ; 

and pitch angle, Irrespective of the seas. In i 

the concour mode the feedback gains ar much 1 

lower and the boat has a greater tend: ■ to ! 

contour the wave profiles. Load environ. ?ts ; 

are fundamentally affected by the control - i ! 

mode. j 

Ill, Full Scale Testing Conditions 

and Instrumentation 

The major portion of the hydrofoil load behavior 
was obtained from strain gage Instrumentation 
placed on the foil structure during "rough 
water" testing of a production JETFOIL cra f t. 

Analytical tools were available for examining 
the statistical load behavior but could not be 
used to confidently establish design criteria. 

The complexity of the environment and the 
absence of sufficient data to verify the 
Impl If led analytical representations necessi- 
tated a strong test oriented approach. Limited 
data available prior to the full scale test 
program Indicated that the analytical estimates 
of foil loads seriously underestimated the 
actual environments. 



JETFOIL CONFIGURATION 
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Figure 2 Illustrates the basic strain gage 
Instrumentation used In the full scale testing. 
After Installation a comprehensive static 
loading program was conducted to load calibrate 
the gages. With the appropriate combination of 
strain gages at a particular foil section, 
shear, bending moment and torsion loads could 
be determined. A multiplicity of gage loca- 
tions allowed determination of distributed 
loads. 

Many of the gages, particularly those on the 
foil surface, were exposed to high speed water 
flow when the boat was follborne. This severe 
environment necessitated careful consideration 
of the strain gage installation technique to 
insure a good test life. Figure 3 illustrates 
the grge Installation used on the foil surfaces 
The strain gage was a weldable gage with an 
integral waterproof stainless steel jacketed 
lead which was secured with retention straps. 

A flexible, strong epoxy was faired over the 
gage and lead. After curing it was sanded to 
provide a smooth contour with the steel surfac 
of the foil. 


Follborne testing in rough water was accomplish- 
ed In the Pactftc Ocean and the Straits of Juan 
de Fuce near Neah Bay, Washington, Total test- 
ing activity was of several days duration per- 
formed at three separate calendar periods and 
encountered a diverse range of sea conditl.ns. 
Operation In sea state five with significant 
wave heights up to 12 feet was accomplished. 

Data In at. sea headings (head, bow, beam, 
quarter and follow) was obtained, 

A fundamental part of the data analysis was not 
only to define basic loads behavior but to 
associate the load information to specific sea 
environments. To accomplish this, the sea 
environment was measured by an onboard wave 
height instrument. This Instrument gave a 
reading of inertial wave profile versus time by 
subtracting measureo boat motion from the sonic 
height sensor measurement taken it the bow of 
theJETFOIL. The wave profile was then 
statistically analyzed to determine the generic 
parameters generally used for sea state 
characterization (e.g., H ^ and T ■ ) . 



FIGURE 2 FOIL SYSTEM STRaIN GAGE LOCATIONS FOR ROUGH WATER LOADS TESTING 
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FIGURE 3 EXTERIOR FOIL STRAIN GAGE INSTALLATION 


IV. Load mechanisms 


Because a hydrofoil ship consists of a normal 
displacement hull plus a foil system which func- 
tions very similar to an airplane wing, the 
variety of loading possibilities becomes quite 
nunerous. The foil systems individually ex- 
perience many types of loading. However, when 
follborne In very large seas, load interactions 
can also occur which involve both the hull and 
foil systems. The primary foil load environ- 
ments can be categorized into three basic load 
mechanisms: wave loads, hull slamminq loads, 
and foil unwetting loads. These phenomena 
produce the basic beam bending and torsion 
fatigue loads whicn dictate the structural 
design. Detailed attention is also required 
for local pressure loading which is required to 
analyze panel tending fatigue environments. 

Have Loads 

The basic wave loading phenomena Is produced 
by orbital wave velocities causing a change in 
angle of attack on the foil and struts, Be- 
cause the foil Is not submerged at an infinite 
depth, free surface effects also result In load 
fluctuations as the f 11 depth changes. As the 
craft's control system attempts to stabilize 
hull motions It Is constantly moving the flaps 
which produce major modification of the loads 
directly produced by the waves. In essence, 
the foils, struts and flaps are experiencing a 
gust environment each time a wave passes. 
Frequency of this wave load phenomena is depend- 
ent on sea state and boat speed but Is typically 
around 1.0 Hertz. Figure 4 Illustrates the 
basic wave environment Imposed on the struts and 
foils. In describing the wave loading phenomena 
It is convenient to distinguish the loads as 
either synmetrlc or antisymmetric. 

The symmetric loads affect only the foils. As 
Illustrated In Figure 4 they result from head 
or following sea wave velocity components. 

Since the craft s control system operates on 
acceleration feedback, the flaps deflect In an 
attempt to null symmetric lift forces. The 
flaps are continuously operating as gust 


alleviators. In the process, however, torsion 
loads are increased because of the difference 
in center of pressure locations between flap 
lift and angle of attack lift. Likewise, a 
severe environment is imposed on the flap mech- 
anisms. 

Antisymmetric wave loads primarily result from 
beam sea orbital velocity components. The beam 
sea components cause side loads on the struts 
which in turn result in antisymmetric lift on 
the foils. Some antisyronetric lift is caused 
by hydrodynamic flow induction between the 
struts and foils but by far the biggest anti- 
symmetric lift contributor is the antisyimetric 
aileron action of the aft flaps. It is this 
aileron action which provides the only source 
of roll stability while the craft is foilborne. 
As the strut side loads produce a rolling moment 
about the craft's center of gravity, the aft 
flaps deflect and product a corresponding anti- 
syimetric lift which negates the induced roll- 
ing moment. 

Hull Slamminq Loads 

When the seas become relatively severe, many of 
the waves encountered will contact the hull. A 
moderate interaction of the wave surface with 
the hull of the hydrofoil is referred to as a 
crest and is dominated by buoyancy forces on the 
hull. This produces a moderate increase of foil 
loading which is not readily distinguishable 
from wave loading. A severe impact of the hull 
results in hydrodynamic slamming pressures which 
dramatically affect the foil loads. Hull slam- 
ming is infrequent compared to wave loading 
frequency but causes large loading fluctuations. 
The hull motions induced by slarrcnlng create 
angle of attack changes on the foils The con- 
trol system and flaps tend to increase the load- 
ing changes because they are operating to coun- 
teract the hull fo-ces by decreasing lift on the 
foils. This type of loading mechanism can be 
particularly severe on the front foil where 
large negative flap angles are concurrent with 
negative angle-of-attacks. Figure 5 illustrates 
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HULL SLAMMING FOIL UNWFPING 




FIGURE 5 FOIL HULL SLAMMING AND UNWETTING FORCE 



foil loads which typically result fra# a hull 
slam event. Strut side loads can also be sig- 
nificant during a slam because of lateral hull 
motion. The typical load cycle Illustrated 1n- 
volves a large negative ltft fluctuation and a 
following positive increnental lift which nor- 
mally occurs after the hull slam as the craft 
stabilizes Itself. Full load reversal of the 
front foil was observed from many hull slams 
during sea trials. 

Foil Unwetting Loads 

When the foil comes sufficiently close to the 
water surface, a ventilation or unwetting of 
the top foil surface will occur because the low 
flow pressures on the upper surface cannot be 
sustained In the vicinity of atmospheric 
pressure. This event, as for hull slaiming, 
occurs In larger seas and although Infrequent, 
causes large load fluctuations. Because condi- 
tions are never completely symmetric, unwetting 
generally occurs on one semispan of the foil as 
illustrated in Figure 5. The semispan which 
unwets loses lift and the other semispan will 
subsequently experience an Increase In lift 
because the flaps are deflected down to offset 
the Initial lift loss and downward hull motion. 
This uwsymmetrlc lift results In large bending 
moments on the strut. For the JETFOIL this 
loading event was more critical to the strut 
than to the foil. In the extreme event where 
the foil completely breaks the water surface 
the lift of both semispans obviously goes to 
zero. This is called a broach. Although 
possible, it is very difficult to achieve un- 
wetting of the aft foil. Foil unwetting is a 
fatigue concern primarily for the front foil 
and strut. 

Local Pressure Loads 

Although so called local panel pressure loading 
primarily results from the basic wave loading 
mechanism, separate consideration of this pro- 
blem Is needed. Factors affecting local 
pressure include wave orbital velocities, flap 
motions, fluctuating waterline on the struts, 
and cavitation. Another consideration for the 
water duct Inlet pressure was Inlet unwetting 
which occurs when the water duct breaks the 
wave surface In large seas. Inlet unwetting 
Is not frequent but large pressure changes are 
involved. 

Near the waterline on the struts wc,ve surface 
movement causes a large number of pressure load 
cycles due to cyclic unwetting and rewetting of 
the flow. This Is a symmetric pressure cycle 
on the strut and causes no net beam loading. 
Magnitude of these pressure cycles is up to 
15 psl and Is of primary concern for the panels 
at midchord and near the trailing edge. Induc- 
tion effects from the foil and local angle of 
attack from beam sea orbital wave velocities 
become Important for areas located a few feet 
below the waterline. These later effects re- 
sult In antlsynmetrlc pressure loading of the 


strut a' (.. tmarily affect the panel* near the 
leadirK The centerline strut panels are 

subject Internal pressure changes of the 
water duct as well as the above external pres- 
sure environments. 

Another local phenomena is cavitation which re- 
sults when the water pressure drops below the 
vapor pressure causing an air cavity, this 
phenomena Is most prevalent on the upper foil 
and flap surface. During rough water operation 
the wave and flap dynamics result In rapid 
changes of pressure which result In "make and 
break" cavitation areas on the upper foil sur- 
face. Test data showed the local panel cavita- 
tion loading behaved in a quasi-steady state 
manner being limited to a certain value when 
cavitation occurred over the panel. The cavi- 
tation environment over major surface areas 
and flaps was also important to overall loading 
distributions. 

V. Development of Fatigue 
Load Design Spectrums 

Data Analysis 

Due to the wide variety of sea conditions and 
total time duration of the full scale testing 
activity. It was possible to build a good 
statistical data base relating foil loads to 
specific environment conditions. The problem 
was reducing the large number of loading cycles 
(on the order of 10,000 for each parameter) 
obtained from testing down to some manageable 
parametric description. Figure 6 illustrates 
the basic procedure used to statistically 
analyze the data. Data samples typically of 
three to five minutes duration were digitally 
analyzed to give the following parameters: mean, 
standard deviation, average of 1/3 (X 1/3) and 
1/10 (X 1/10) highest peaks. Detailed informa- 
tion on the distribution of the entire popula- 
tion of load cycle events could also be obtain- 
ed In the form of cumulative probability 
distributions or exceedance plots. The desired 
objective was to establish a simple parametric 
statistical distribution which would satisfac- 
torily approximate the distribution of the test 
data. If this could be done, then the only 
statistical parameters needed to provide a 
complete characterization of the load behavior 
were the standard deviation and the average 
load cycle frequency. This objective was 
generally met for all sea conditions where the 
occurrence of major nonlinear lead events was 
infrequent. The load cycles arising from hull 
slamming and foil unwetting were much larger 
in magnitude and distinctly more lonsymmetrlc 
In character than the wave loads. Thus, It 
became necessary to hand count and tabulate the 
large nonlinear loud events Individually, A 
parametric description of these loads was not 
feasible, thus necessitating the development 
of empirical tabular data to characterize 
their frequency and magnitude, fortunately, 
the wave loads could usually be characterized 
by the familiar Rayleigh distribution 
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known to be valtd for the distribution of peaks 
of a random variable from a linear, narrow-band* 
ed. Gausian process. In those cases where a 
Rayi«1gh assumption was not satisfactory, using 
a para.-'etrlc combination of the Rayleigh and 
exponential distribution 


P(X t >X)-e" X/ff (2) 


The technique of spectral moments (2) Is pre- 
ctse only for stationary conditions. For non- 
stationary conditions a numerical approach was 
requtred to correctly translate the wave spec- 
trim) equation Into a moving system and provide 
an accurate correlation between Tstg and f^. 

Reduced statistical load parameters were then 
cross plotted against Tslg and Hslg. Loads 
for sea conditions not specifically encountered 
in testing were defined by utilizing the appar- 
ent data trends and incorporating trend data 
from analytical simulations. Thus, a data base 
was formulated which could define load spee- 
trums for any design environment. 

Design Environment 


provea to be useful. 

Sea state characterization was obtained by 
calculating Hslg and Tslg from the wave height 
instrument data: Calculation of Hslg from the 
wave height Instrument was a straight-forward 
process, however, determining Tslg was more In- 
volved because Tslg is a stationary parameter. 
Since the JETFOIL provided a nonstationary 
platform for the wave height sensor, some trans- 
formation of the data was required to estimate 
the generic Tsig description of the wave's 
periodicity. The wave height sensor data was 
analyzed to define f for the moving craft 
system. This could "be related to Tslg by using 
standard statistical relationships developed 
from theoretical wave spectrums Cl ) , Because of 
the nonstationary measurement condition, a pre- 
cise mathematical relationship between f and 
Tslg was not possible. 



FIGURE 6 LOADS TEST DATA STATISTICAL ANALYSIS 


Selecting a specific design environment had to 
be done somewhat arbitrarily as no traditional 
guidelines existed. Tabulations of visually 
estimated wave characteristics from many areas 
around the world were available! 3 ) , v 4 I . How- 
ever, the question of what area to select as 
the representative design condition was a 
troublesome one. What routes the JETFOIL craft 
would eventually encounter could cover almost 
any coastal or Inter-island area in the world. 
The problem was eventually managed by examining 
the wave data from several "rough sea" areas in 
the world and then enveloping the wave height 
and period data. This defined a so called 
"composite sea state" for structural design. 

Some of the most severe sea conditions were 
then truncated out because they were sufficient- 
ly rough to prevent sustained foil borne opera- 
tion of the JETFOIL craft. Commercial opera- 
tion In extremely high seas was obviously un- 
feasible. Figure 7 shows the long term wave 
height distribution for this composite sea 
model. The complete model included a large 
number of Tsig and Hsig combinations, each with 
an allotted time percentage of the total opera- 
tional design life. 



FIGURE 7 COMPOSITE SEA MODEL HAVE 
HEIGHT DISTRIBUTION 
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Wave Loading Spectrums 

Figure 8 provides front foil wave load exceed- 
ance data for several sea conditions. The 
abscissa is a normalized load magnitude for 
fully reversible symmetric load cycles. The 
load value used for normalizing the data is 
based on the expected maximum rough water load. 
The three bow sea spectrums for specific wave 
heights show the dramatic increase of load from 
a moderate to a rough sea. For example, the 
exceedance rate of a load value of 0.1 went 
from 1.5 for Hsig equal 4.0 feet to 350 for a 
Hsig of 9.0 feet. This is indicative of the 
sensitivity of fatigue damage to rougher seas 
which will be more clearly shown in Sectton VI. 
Also shown is a design spectrum which is based 
on the weighting of the composite sea of Figure 
7. 

An appreciation of the severity of the hydro- 
foil fatigue environment is gained from a 
comparison with commercial jet airplane gust 
loads. For example. Figure 8 shows an 
exceedance rate of 24 per hour at a load 
factor of 0.125 for the composite sea 
environment. A typical design exceedance 
rate at the same equivalent gust load for 
a commercial jet liner is on the order of 
a few per hour. This does not consider 
the hull slamning and unwetting loads 
experienced by a hydrofoil which can be more 
damaging in terms of fatigue than the wave 
loads. 


The change in character of the curves at low 
exceedance rates is due to the partial use of an 
exponential distribution for larger seas. In 
low to moderate seas the Rayleigh distribution 
was valid, but in higher seas a few of the wave 
load parameters exhibited a greater frequency of 
high level loads than predicted by a Rayleigh 
distribution. This Increased activity of high 
ragnltude wave loads is attributed to the major 
"rough sea" nonlinearities which increase the 
loading energy (e.g., free surface depth effects 
and cavitation effects which tend to limit flap 
effectiveness). Assuming a portion of the load 
cycle events to be exponentially distributed 
provided a higher rate of large magnitude loads 
and better test correlation. 

Figure 9 shows the sensitivity of aft foil anti- 
symmetric lift to sea heading. In general, head 
and bow seas were more severe than beam, quarter 
and follow seas for all load parameters. Anti- 
symmetric loads showed the most pronounced and 
surprising differences due to heading changes. 
From an overall point of view the bow sea Is the 
worst heading because it has the highest anti- 
synmetric lift and strut side loads and also has 
symmetric lift loads only moderately less than 
those of a head sea. The antisymmetric load 
levels indicated for a head and follow sea are of 


particular Interest because theoretically, one 
would predict an antisymmetric load level of 
zero, Theoretical wave energy models such as 
the Bretschnelder (1) spectrum are based on one 
dimensional seas with long crested waves. This 
condition Is never achieved in real seas because 
some multidirectional characteristics are always 
present. For example, a typical sea Is wind 
blown waves from one direction (e.g., beam) on 
top of an underlying swell from another direc- 
tion (e.g., head). This multidirectional proper- 
ty of real seas is an important factor in the 
load behavior. However, another suspected fac- 
tor is the un symmetric behavior of cavitation 
patterns which occur on the foils. Cavitation 
characteristics are very se-sitive to the foil 
contour and because of manufacturing tolerances 
a truly synmetric foil is never achieved. This 
results in unsymmetric hydrodynamic loading in 
head and follow seas. Both of these effects 
present analytical difficulties for a load 
simulation approach, the latter being particu- 
larly difficult to predict. In any event, the 
conslaeration and accurate identification of 
load effects related to sea heading is very 
important. A real hydrofoil route must consider 
a mixture of at least two different headings. 

For the JETF0IL a 50% bow, 50% quarter was 
considered to be the most realistic for a design 
condition. 

FigurelO illustrates the sensitivity of foil 
loads due to the control system mode. In the 
platform mode, the control system has higher 
feedback gains which provide more effective 
"gust alleviation" due to increased flap motions 
counteracting the orbital wa ve velocity forces. 

In the contour mode the feedback gains are much 
weaker and thus gust load alleviation is less 
effective causing nigher loading. Although the 
higher feedback gains do provide lift load alle- 
viation, foil torsion loads and flap loads are 
increased. 

Hull/Slamming and Unwetting Loads 

Figure 11 Illustrates event frequency for hull 
slamning and foil unwetting loads. The curves 
are basically empirical but the predicted 
frequency rates in the low seas have a low 
confidence level due to a lack of sufficient 
long term statistics. Using a simplified 
analytical approach provided some useful trend 
data for the hull slamming frequency but the 
foil unwetting event could not be reliably esti- 
mated in any analytical manner. Without proper 
testing and strain gage Instrumentation the 
severity of these loads could not have been 
established. 

As for displacement hulls, head and bow seas are 
nuch more severe for hull slamming. Only in the 
roughest seas does the JETF0IL experience sig- 
nificant slamning rates in a beam, follow or 
quarter heading. As well as wave height, the 
slamming rate is also very sensitive to Tsig. 

For steep waves with a short Tsig. the rate of 
slamning may be high while long swells with a 
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large Tsi 9 will cause a relatively low rate of 
slamming. Figure 11 does not Illustrate the 
sensitivity to Tslg but Is based on a typical 
sea for the given wave height. The frequency 
of foil unwetting loads Is not as sensitive to 
sea conditions as slamming. Effect of sea 
heading is particularly different. Foil un- 
wetting loads were as common in a following or 
quartering sea as head or bow. 

VI. Fatigue Damage 

The prediction of structural fatigue damage is a 
complex subject in its entirety. However, a 
basic understanding of the structure's fatigue 
sensitivity to different types of load spectrums 
is fundamental to accurate load analysis of the 
hydrofoil craft. By necessity, simplification 
of the fatigue loading histories is required to 
put the total fatigue analysis effort into a 
manageable form. An understanding of actual 
fatigue damage sensitivity puts the load analy- 
sis requirements into proper perspective. 

Fatigue damage or life after a complex loading 
pattern is generally predicted by using the 
linear cumulative damage theory, more widely 
known as Miner's rule(S). The formal statement 
of the linear damage theory is that the summa- 
tion of all fractions of life, or damage equals 
one at failure. 

L Vl+li* «1 (3) 

N i n 2 

where Ni are the single level fatigue lifes at 
each given level of loading and n^ are the actu- 
al number of load cycles at the specified level. 
If the damage number is less than one, it indi- 
cates the fraction of life “used up" for that 
particular loading spectrum. The Ni can be cal- 
culated in convenient equation form by using 
linear S-N and constant life plots (Goodman Dia- 
gram) illustrated in Figure 12. These linear- 
ized curves are approximations to the actual 
complex fatigue life properties and can be con- 
veniently based on test data at an appropriate 
characteristic life (e.g., 10 5 cycles). 


Thus: 

„ =10 (5-!oyYi) 
”i iU logs 

(4) 

where 

V f *i 
f A-f , \ 

(5) 


y-v ) 
\ %/ 


This technique is used to define the numbers in 
Table 1 which show the fatigue dan, age sensitivi- 
ty to sea conditions of a typical area on the 
bottom surface of the front foil. The tabu’ated 
numbers are arbitrarily normalized to the maxi- 
num net damage number which is illustrated for a 
Hsig of 9.0 feet. The comparisons vividly demon- 
strate the sensitivity to sea conditions. Oper- 
ational requirements were based on a conriercial 
route which consisted of a takeoff and landing 
every hour and a heading mix of 50% bow, 50% 
quarter. The takeoff/landing cycle was defined 
by the maximum expected wave load added to the 
lg steady lift load. A ground-air-ground cycle 
is similarly defined for commerc'al airplanes. 

For short range routes ot airplanes this cycle 
typically dominates the fatigue damage. Fcr 
hydrofoil operation Table 1 indicates that it is 
insignificant unless only moderate seas are en- 
countered. In fact, the comparisons point out 
that, for a route which encounters “rough" seas 
a significant portion of the time, the moderate 
seas are a negligible contributor to the overall 
fatigue damage. These net damage numbers show a 
sea of 6.5 feet significant wave height is over 
9 times as damaging as a 4.0 feet significant 
wave height and 9,0 feet is 56 times as damaging! 
Also evident is the major effect of the infre- 
quent but large loads from slamning and foil un- 
wetting, These large loads are important con- 
tributors to fatigue damage for nearly all foil 
and strut structure but are most dominate for 
the front foil system. Irrespective of location, 
all primary foil system structure has the same 
qualitative fatigue life sensitivity to sea state. 


These comparisons vividly point out the basic 
fatigue analysis problem for hydrofoils 
operating in the open seas; structural fatigue 
damage is heavily influenced by those load 
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FIGURE 12 LINEARIZED FATIGUE LIFE CHARACTERISTICS 
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events which are the most Infrequent and comp- 
licated to predict and those sea environments 
for which It Is the most difficult to obtain 
accurate statistics. From the hydrofoil 
design point of view, a great need exists to: 
(1) obtain better statistics on sea enclron- 
ments and (2) develop reliable analytical 
simulation techniques which represent basic 
load behavior and accurately predict load event 
frequency for a wide range of environments. 

Although the subject of flaw growth will not be 
discussed In any detail, the author wishes to 
point out that flaw growth considerations are 
important to the overall hydrofoil fatigue 
damage problem as well as classical fatigue. 
This Is particularly true for welded structure. 
The qualatatlve fatigue weighting of sea 
environment discussed above Is also applicable 
to flaw growth. For sizing some structural 
members flaw growth Is the dominant considera- 
tion. 


LOAD PARAMETER 

h sig, feet 

mm 

6.5 

ESI 

TAKE OFF/LANDING CYCLE 
(1 PER HOUR) 

0.008 

0.0143 

0.025 

WAVE LOADS 

0.0099 

0.0551 

0.162 

UNWETTING AND HULL 
SLAMMING LOADS 

0 

0.0954 

0.813 

DAMAGE SUMMATION 
OF ALL LOAD CYCLES 

0.0179 

0.165 



TABLE 1 FRONT FOIL NORMALIZED FATIGUE 
DAMAGE COEFFICIENTS 

VII, Analytical Simulations 

Although the emphasis on test data to this point 
suggests that analytical efforts were nil, this 
was not the case. Two digital computer simula- 
tion models were utilized. The primary one 
which provided a basic statistical analysis 
tool was a flve-degree-of-freedom (heave, sway, 
pitch, yaw and roll) pertubation model based on 
the classic spectral response analysis method 
for linear, gaussian, random processes 

$ 0 (<o)»|H(j W )| 2 S t (o>) (6) 

where S 0 (u) is the output spectrum, Si(u) Is 
the Input spectrum and H(jw) Is the complex 
frequency response function for the chosen out- 
put variable. Si(w) Is defined from a theoreti- 
cal sinusoidal wave spectrum model'- 1 ’ , H(Jw) 
is computed from the hydrofoil craft model which 
may Include any degree of sophistication in 
hydrodynamic lift and drag characteristics of 
the foil surfaces, control system logic and 


structural properties as long as the model is 
mathematically linear. It 1s Implicit that the 
sea model characteristics are also linear, 
Obvtously, this model could not be used to 
analyze the Inherently nonlinear slamming and 
fall unwetting loads but was potentially valid 
for foil wave load analysis. Original load 
magnitudes using this model were seriously low 
when correlated with test data. By better 
parameter Identification of items like free 
surface effects, structural f'exiblllty and 
flap effectiveness the load estimates were 
increased to provide better correlation with 
the wave load data for bow and head seas. How- 
ever, some of the estimates were still de- 
ficient, particularly for larger seas. Predic- 
ted load behavior for other headings correlated 
poorly with test data. A particularly divergent 
comparison is for the antisymmetric foil loads 
in head and follow seas (Section V). 

Further Improvements are possible with this 
approach, For example, multidirectional sea 
spectras could be employed and piecewise 
linearization of nonlinear hydrodynamic para- 
meters and wave characteristics Is feasible, 
tony further potential refinements require 
empirical data for a realistic definition of 
parameters. 

The second model was a complex time domain 
simulation whtch Included nonlinear hydrodynamic 
force representations and also a slamning model 
to predict pressure distributions on the hull 
due to wave impact. This model was of some 
help in characterizing foil load behavior 
during slamming events. However, Its complexity 
and cost was prohibitive far use as a tool to 
provide comprehensive statistical data. No 
theoretical representation of the foil unwetting 
event was available. 

VIII. Concluding Remarks 

The hydrofoil fatigue load environment from 
rough water operation is complex, severe and 
difficult to predict. Primary loading of the 
foil system structure involves three distinct 
mechanisms: wave loading from orbital velocities 
and foil depth changes, hull slamming and foil 
unwetting. Loading distributions and local 
panel pressures are complicated by cavitation 
and constant changes of both waterline location 
and wetted area of the struts. Major non- 
linearities occur for all of the load mechanisms 
which limits the usefulness of current state-of- 
the-art theoretical approaches to prediction. 
Probabilistic descriptions of wave height, wave 
period and sea heading are required to predict 
the spectrum of loading events. A lack of 
accurate statistics on sea environments for : 
given area Is a general problem. This prob'em 
is compounded for load spectrum predictlor 
because the most severe and infrequent seas 
cause the majority of the fatigue damage. 

These considerations necessitated a strong 
dependence on test data to obtain a timely 
definition of loads criteria for the JETFOIL. 
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A comprehensive amount of test data Is required 
to define the spectrum of hydrofoil load events. 
The large number of wave load cycles can gener- 
ally be represented by the familiar Rayleigh 
distribution and In cases where a Rayleigh 
distribution Is Inadequate, a parametric combi- 
nation of Rayleigh and exponential distributions 
can be used. The complicated nature and distri- 
bution of hull slamming and foil unwetting loads 
necessitates a completely empirical definition. 
Because of Implementing a large scale test 
program conducted In several types of sea 
environments, the JETFOIL fatigue environments 
were adequately defined. However, the problems 
and expenses Involved with obtaining high 
confidence fatigue load spectrums from testing 
points out the critical need for more sophisti- 
cated but efficient analytical approaches to 
predicting the loads. This Is especially 
Important for preliminary design. 

The state-of-the-art theoretical approaches to 
analyzing ship responses In a random sea are 
basically valid for the typical wave loading 
environment of the hydrofoil where the non- 
linear loads of hull slamming and foil unwetting 
are infrequent. However, further refinements 
in these linear modeling techniques are needed. 
Larger waves Introduce significant nonllnearl- 
tles and multidirectional seas can invalidate 
the typical one dimensional sea representation. 
From the viewpoint of the hydrofoil structural 
analyst, there Is a need for more direct and 
statistically accurate measurements of sea 


environments. Experimental measurements of 
wave orbital velocity distributions In many 
types of sea environments Is needed for verifi- 
cation and upgrading of the current energy 
representations of the sea surface. This Is 
particularly Important In sea environments 
which are confused In direction and often 
common In coastal areas where hydrofoils 
predominantly operate. Future developments 
In modeling techniques for hydrofoil loads In 
a random sea should enhance the ability to 
concisely define the complicated load environ- 
ments. 
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EVALUATION OF ROTOR-BEARING 
SYSTEM DYNAMIC RESPONSE TO UNBALANCE 


Roger E. Thaller and David W. Ozimek 
Aeronautical Systems Division 
Wright-Patterson Air Force Base, Ohio 


A complete Investigation of the vibration environment within 
air conditioner rotating machinery referred to as an Air 
Cycle Machine (ACM) was needed to effectively Increase ACM 
reliability. To assist in the selection of design changes 
which would result in improved ACM performance, various 
design modifications were incorporated into a baseline ACM 
configuration. For each design change, testing was conducted 
with the best balance achleveable (baseline) and with various 
degrees of unbalance. The purpose was to establish relation- 
ships between unbalance (within the context of design changes) 
and the parameters associated with design goals. The results 
of rotor dynamics tests used to establish these relationships 
are presented herein. 


BACKGROUND 


Overview 

Air Cycle Machines (ACMs) from Air Force 
Ground Support Air Cycle Air Conditioner Carts, 
used for cooling and heating tactical aircraft 
on the ground, continue to experience an 
unacceptable rate of failure. At times these 
failures have reduced the availability, and 
have seriously Jeopardized the mission, of 
various aircraft. ACMs have been returned for 
failure investigations and/or maintenance over- 
haul following rotor seizure or any condition 
which, by experience, indicated that rotor 
seizure was imminent. Imminent seizure was 
determined from such conditions as (1) lubri- 
cating oil turning dark, (2) smoke exhausting 
from the compressor discharge and (3) metallic 
particles appearing in the oil. Vibration was 
considered to be a possible cause of failures. 

Air Cycle Machine 

As shown in Figure 1, the ACM is basically 
a rotating assembly that consists of a turbine 
wheel (C) and fan (E) fastened to opposite 
ends of a common shaft (D). In general, high 
temperature and pressure air passes through the 
nozzle (B) and is expanded through the turbine 
wheel to a lower temperature and pressure air 
as energy is extracted in the form of work and 
heat. Work is done on the rotor. Work trans- 
formed into heat passes through the angular 


contact ball bearings (H). The ACM rotor con- 
sists of the shaft (D), compressor wheel, tur- 
bine wheel, oil slingers (N), and the inner race 
of each bearing. The remainder of the ACM 
structure with the exception of the seals and 
bearing balls, can be considered the stator. 

The rotor is housed in one section of the stator, 
the cantilevered cartridge (J). The cartridge 
is fastened to the turbine end of the main 
housing (0). The bearings, which provide a 
dynamic interplay between shaft and cartridge, 
are preloaded through action of the center 
spring (P) against each sleeve (Q). The bear- 
ings are cooled by action of the oil system 
which circulates MIL-L-6085 oil. The oil sys- 
tem consists of Dacron wicks (K), capillary 
tubes (L), the ACME threaded portions of the 
shaft (M), and oil slingers. Oil is fed to the 
shaft by capillary action through wicks and 
tubes and by pumping action generated by the 
screw threads. Oil is then forced through the 
bearings by the centrifugal action of the oil 
slinger. 

The ACM is a primary component of the Air 
Conditioner system. A general schematic, of 
this system Is provided as Figure 2. The sys- 
tem input and output air conditions dictate 
the ACM input and output air conditions which 
govern the rotor dynamic response. Table 1 
provides typical measured values for the 
governing input and output air conditions. A 
maximum ACM operating speed of 525 Hz (31,500 
rpm) is attained at the 36°C db/SO^C wb (97°Fdb/ 
87°F wb) ambient temperature condition. 
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Fig. 1 - Cross Sectional View of Air Cycle Machine 


MOISTURE 



Fig. 2 - Simplified Schematic of Air Cycle Air Conditioner 
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TABLE 1 


AIR CONDITIONER PERFORMANCE DATA* 


Ambient 

Bleed Air 
Into Air 


ACM 

Turbine 


Air Conditioner 

Temperature 

Conditioner 



INLET 

OUTLET 

O 

utlet 

db/wb 

PRES 

TEMF 

PRES 

TEMP 

PRES 

TEMP 

TEMP 

Relative 
Humid! ty 

-c/'c 

kPa 

°C 

kPa 

°C 

kPa 

*C 

°C 

X 

-3. 9/1. 7 

310 

171 

276 

6.7 

165 

-10 

10 

45 

10/7.2 

307 

193 

265 

14 

165 

3 

11 

45 

27/18 

286 

213 

238 

28 

10 

1.7 

11 

55 

36/30 

262 

221 

214 

38 

62 

13 

14 

53 

38/29 

255 

224 

207 

43 

59 

14 

17 

60 

54/24 

248 

232 

172 

52 

48 

15 

18 

52 


♦All values within this table are normalized for atmospheric pressure 
of 101 kPa and a delivery pressure of 14 kPa. Delivery airflow Is set 
at 0.S3 kg/sec. 


Rotor Bearing and Cartridge System Dynamic 
Performance 

The dynamic behavior of the Air Cycle 
Machine rotor was theorized to be the major 
cause of ACM deterioration and ultimate fail- 
ure. Rotating systems tend to whirl or whip 
at frequencies which are characteristic of the 
r^tor system design. Operation at or around 
these frequencies could lead to destructive 
results. These vibration type phenomenon can 
be either synchronous or non-synchronous in 
nature. Synchronous whirl is characterized by 
a slow buildup of transverse vibration with 
the fundamental frequency equal to a particu- 
lar running speed of the rotor system. This 
particular running or synchronous speed is 
called the critical speed. Synchronous whirl 
can be primarily attributed to mass unbalance 
of the rotating system which could result from 
such factors as 'nequal load distribution, 
unsymmetrical stiffness of shaft or bearlng(s), 
oil friction in bearings, and hysteresis 
damping in the shaft. Non-synchronous whirl 
is characterized by self-excited vibration 
induced by various elements in the rotor sys- 
tem. Such elements as seals, aerodynamic 
pieceparts and fluid film bearings can induce 
these vibrations thereby causing rotor insta- 
bilities which directly result In the non- 
aynchronous whirl condition. 

Rotor systems containing rolling element 
bearings, such as this ACM, are particularly 
sensitive to system balance when operating at 
or within a range of one or more system 
critical speeds. Previous development tests 
on an earlier generation of this ACM estab- 
lished that at least one first critical speed 


would be very near the maximum ACM operating 
speed. The interrelationships between these 
critical speeds and rotor balance were similarly 
speculative. 

aCM Improvement Project 

An extensive ACM Improvement Project was, 
therefore, undertaken by thp PRAM Program Office 
of the Air Force to thoroughly investigate fail- 
ures by -nalytlcally and experimentally identify- 
ing ACM deficiencies and by determining necessary 
corrective actions. One of many ACM designs 
currently used in the Air Force was selected for 
study. The project continues to involve a com- 
bined effort between Air Force management, 
engineering and logistics organizations. Two 
Air Force contracts were awarded to meet project 
objectives: Mechanical Technology Inc. (MTI - 

Latham NY) was require-! to perform a majority 
of the analytical efforts and United Aircraft 
Products Inc. (UAP - Forest OH) was required to 
simultaneously perform a complete test effort in 
conjunction with independent analytical studies. 

The principle objectives of the ongoing Im- 
provement Project continue to be to provide an 
ACM which: 

1. physically and functionally interface 
with a given Air Conditioner Cart, 

2. meets or exceeds performance and 
reliability requirements of the governing speci- 
fication, 

3. is compatible with the current high 
pressure bleed air sources used by the Air Force, 
and 
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4. 


decreases system life cycle cost. 

The project was divided Into four tasks. A 
long term study and test effort was determined 
necessary to Identify both ACM failure cause(s) 
and optimum corrective actions. This effort 
was divided into three tasks in order to estab- 
lish a baseline (Task 2), identify optimum 
design change(s) (Task 3), and demonstrate the 
acceptability of the optimum design through 
reliability testing (Task 4). To protect exist- 
ing operational ACMs from further failure during 
this long term effort, a short term test effort 
was required to size an orifice which would 
reduce the ACM maximum operating speed. This 
effort (Task 1) is r.ot discussed herein. 

In Task 2, baseline performance tests were 
conducted on two ACMs (S/N 312 and 313) to pro- 
vide engineering baseline data as well as to de- 
fine ACM deficiencies. Analytical studies were 
simultaneously performed to establish correla- 
tion between theory and test results. Details 
of these Task 2 efforts that did not relate to 
unbalance were presented at the 24th Annual 
Technical Meeting of the Institute of Environ- 
mental Sciencles. The results that relate 
response to unbalance are pi sented herein. 

In Task 3, design changes wc.e analytically 
assessed and selected for test. The resulting 
hardware changes were then subjected to a com- 
plete testing sequence. After identifying an 
optimum design approach, the optimum design was 
further tested in Task 4. A portion of the 
Task 3 and 4 results are presented herein. 


ACM DYNAMIC RESPONSE TESTS 


Test Procedures 

At the completion of the Task 1 interim fix, 
the baseline tests were initiated. The same 
basic test instrumentation was used for all the 
tests from baseline through reliability testing. 

A more complete description of the instrumenta- 
tion and processing technique used throughout 
these efforts Is presented in the following sec- 
tion. Two ACMs (S/N 312 and S/N 313) were 
tested; each with 13 proximity probes, two 
accelerometers and a variable quantity of ther- 
mocouples and pressure sensing devices. Three 
of the proximity probe sensors provided the most 
significant information required for rotor un- 
balance response analyses. Referring to Figure 
1, Y1X (U) provided a measure of shaft rotation; 
Y1Y (R) and Y5Y (T) provided measures of verti- 
cal displacement from the shaft at the com- 
pressor end and turbine end respectively. Data 
was obtained during acceleration of the rotor 
shaft. For each acceleration run, the data 
analysis produced plots of Y1Y, Y5Y and the 
phase angles between the pairs Y1X-Y1Y and 
Y1X-Y5Y. 

During the tests, unbalance was generated 
by adding calibrated threaded nuts to the com- 
pressor end of the shaft for baseline and de3lgn 


change tests. One nut was used for each un- 
balance condition beyond the baseline unbalance. 
Calibrated nuts were installed on the turbine end 
of the sh. .t during the first part of the reli- 
ability tect. After this first part, the ACM 
was returned to a baseline unbalance. The ACM 
was then allowed to become unbalanced during 
operation. 

During che tests, data were obtained during 
acceleration and deceleration of the rotor shaf 
using a preselected range of ACM input/output 
air condidoui- . Data were also acquired du'lng 
ACM operation at a constant speed. 

Baseline Tests 

In order to provide a starting point or 
baseline, the existing ACM configuration was 
chosen for test and analyses. The configuration 
conslste ! of an elastomeric O-ring and 120N 
(27 lb) bearing preload. The 0-rlng was located 
between the cartridge and diffuser housing 
saddle located at the comp' ssor eno of the ACM. 
Preload was provided by a spring located between 
the shaft and cartridge at the center of the 
rotor system. The spring provided load to the 
bearings by transmitting the necessary force 
through sleeves mounted between the spring and 
each bearing. 

Baseline test resul.s showed *• cenpected 
differences in the j.oca iu-s of both the first 
and second critical speed between the two test 
ACMs. Differences of at least 50 rps (3000 rpm) 
were recorded. The resjlting analysis indicated 
that the ACM first critical speed Is primarily 
a function of cartridge stiff’ and the second 
critical speed is primaril-' -ction of 
bearing rotor stiffness. 

Baseline performance analyses 

also identified deficlenc c pertairlng to the 
rotor-bearing and cartridge system of the ACM. 
Analysis of the baseline data Indicated that as 
a minimum the following three design goals be 
established: 

1. Minimize shaft excursions tfcjlting 
from synchronous rotor response. 

2. * Move the machine second critical speed 
15Z or more above the maximum operating speed 
(maximum operating speed Is approximately 525 rps 
[31500 rpm]). 

3. Minimize shaft excursions resulting 
fiom subsynchronous rotor response. 

Design Change Tests 

The baseline test results ltd to the de- 
velopment of design changes, each design change 
was an attempt to achieve at leust one of the 
desired design goals. A 178N (40 lb) bearing 
preload was determined necessary to sufficiently 
increase the second critical speed well above 
the maximum operating speed. Adding either the 
dual 0-ring or the C-rlng was to aid in damping 
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the rotor vibration thereby decreasing rotor 
excursions over the range of ACM speeds (in par- 
ticular the first critical speed). An Increase 
in cantilevered cartridge support, accomplished 
by changing the four bolt circle flange configu- 
ration to seven, was to add stiffness to the 
cartridge, thereby further reducing rotor excur- 
sions while increasing the first critical speed. 

To achieve the design goals, a variety of 
ACM configurations was tested. Each configura- 
tion was a combination of design changes as 
identified in Table 2. The first configuration 
was anal>tlcally determined to provide the opti- 
mum combination of these design changes. The 
second configuration was determined to be the 
optimum design from a stiffness point of view. 
Neither configuration demonstrated the proper 
combination of desired response. Therefore, two 
additional combinations were tested in Phase III. 
The Phase IIIB configuration demonstrated the 
type of performance desired. Therefore, this 
configuration was further tested in another ACM 
(Phase 1VA) to determine if the dynamic perfor- 
mance (in particular the rotor signature) was 
repeatable. Complete repeatability could not be 
obtained. However, the Phase IV configuration, 
tested to determine if the O-ring in *’ j base- 
line configuration, had served a useful purpose, 
demonstrated far superior performance over any 
previous configuration. (Baseline test data 
indicated that the 0-ring had no effect on the 
performance of one of the baseline ACMs.) The 
Phase IVB design was therefore selected for the 
lengthy reliability test. 

Reliability Test 

To demonstrate that the optimum design was 
capable of continual and reliable on-off opera- 
tion in a variable envi-onment. the optimum 
design was subjected to a sequential reliability 


test. One ACM was to operate in a cyclic on-off 
fashion under a variety of ACM input and output 
conditions until an MTBF of 435 hours was 
achieved. At least 11 on-off cycles of each 
combination of input-output conditions identi- 
fied in Table 3 were to be tested (a total of 
528 cycles). The results of this test as well 
as baseline and design change tests were analyzed 
to determine dynamic response to unbalance. 


TEST INSTRUMENTATION AND DATA PROCESSING 


Test Instrumentation 

The rotor dynamics investigations were 
conducted by AFFDL using the instrumentation 
shown in Figure 3. AFFDL used proximeters with 
signal conditioning equipment between proxi- 
meter output and tape recorder input. The 
Automatic Gain Ranging Amplifier Units (AGAUs) 
immediately downstream from the proximeters 
were used to increase signal voltage levels 
out of the background ambient noise levels and 
carry the signals over relatively long cables 
to the test monitoring station. The second 
series of AGAUs was used to increase or de- 
crease signal levels to an acceptable range 
for analog tape input. Time code was recorded 
on tape along with all data signals for use in 
data editing and analysis. 

Signals from each pai r of the shaft and 
proximity probes were continuously displayed 
on dual beam oscilloscopes in order to visually 
check rotor orbits. Signals from each of the 
remaining probes were also selectively moni- 
tored. In addition, for the AFFDL test portion, 
selective signal activity aver the entire 
frequency (speed) range of the ACM was dis- 
played on the oscilloscopes. 


TABLE 2 

DESIGN CHANGE TEST CONFIGURATIONS 


Test 

Phase 

ACM S/N 

Configuration 

i 

313 

178N (40 lb) bearing preload; dual elastomeric 
O-rings between cartridge and diffuser housing. 

ii 

313 

7 bolt flange mounting portion of cartridge; 
metallic C-ring between cartridge and diffuser 
housing. 

IIIA 

312 

Increased cantilever support as in Phase II 
combined with 178N bearing preload. 

IIIB 

313 

178N bearing preload and metallic C-ring 
between cartridge and diffuser housing. 

IVA 

312 

Same as IIIB. 

j IVB 

313 

178N bearing preload and clearance between 
cartridge and diffuser housing. 
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TABLE 3 


RELIABILITY TEST CONDITIONS 
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Fig. 3 - Air Cycle Machine Rotor Dynamics Test Instrumentation 
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Fig. 4 - Air Cycle Machine Data Processing Methods 
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For this phase of testing, data were 
recorded as each of two ACMs (both of the same 
basic design) was accelerated then decelerated 
over a range of rotor unbalance conditions. Com- 
pressor air inlet temperatures were also varied. 
Data from the acceleration and deceleration runs 
were used to locate the critical speeds or syn- 
chronous whirl ftequencies. Maximum rotational 
speeds of 583 rps (35,000 rpm), well above 
maximum operating speed, were monitored during 
this exercise. Each ACM was also subjected to a 
series of constant speed runs at a series of un- 
balance and compressor air inlet temperature 
conditions. The ACMs were simply held at a con- 
stant speed as data were recorded for a short 
(10 - 15 second) time interval. Speeds for this 
exercise were selected from a range of speeds in 
the area of the critical speeds observed from 
acceleratxon/deceieiation runs. Data from con- 
stant speed runs were also used to located any 
nonsvnchronous frequencies which could poten- 
tially result in additional rotor excitement. 

Data Processing 

The bulk, of the data processing effort was 
conducted by AKFDL. Processing of the test 
data was performed according to the scheme 
shown in Figure a. The analog time histories 
were processed using a real-time spectrum 
analyzer operated in the peak mode. For each 
recorded time history, consecutive segments of 
data were converted to displacement spectra 
(with a 2 Hz bandwidth) and the maximum value 
at each frequency was updated. The final 


spectrum, consisting of the highest values at 
each frequency, was plotted. Y1X provided a 
measure of shaft rotation; other sources pro- 
vided measures of shaft displacement from the 
axis of rotation. In particular, Y1Y and Y5Y 
provided measures of vertical displacement of 
the shaft at the compressor end and turbine 
end respectively. 

The analog signals acquired from each 
proximity probe during rotor shaft accelera- 
tions were also low-pass filtered at 1C00 Hz 
and digitized. Narrow band spectra were pro- 
duced using a hardware lmpelemntation of the 
Fast Fourier Transform. Displacement spectra 
were acquired for each of the 13 sensors used 
during test and cross power spectral densities 
were obtained for each nair, l.e., Y1X-Y1Y and 
Y1X-Y5Y. The bandwidth for each spectrum ras 
10.18 Hz. For each Interval of .7858 seconds 
during the shaft acceleration test condition, 
a complete set of spectra (displacements and 
cross-power spectral densities was obtained. 
Each spectrum is the average of the outputs of 
lo Fast Fourier Transforms. 

The spectra were processed using a Fortran 
IV program in the following manner. The Y1X 
spectrum was Inspected to determine the shaft 
cpm during each .7858 second Interval. The 
frequency corresponding to that shaft rpm was 
stored. For each of the sensors, the software 
produced a plot of the displacement and the 
phase with respect to Y1X versus shaft speed. 
Figure 5 Is a typical plot. 


Fig. 5 - ACM Rotor Dynamics Data from Analog-to-Digital Processing Methods (Phase 1IIA) 




Fig. 5B - Rotor Turbine End Activity 
Probe Y5Y) 
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ACM DYNAMIC RESPONSE RESULTS AND ANALYSIS 


Analysis of test results involved careful 
Inspection of the curves of displacement and 
phase versus frequency. Figure 5 is a typical 
set of curves. Dlsplacenent peaks and large 
phase changes were used as Indications of criti- 
cal speeds. For each test unbalance condition 
critical speed, the dlsplacenent at tbat speed 
and the imbalance was then identified and tabu- 
lated. Tables 4 through 8 provide these results. 
For particular unbalance conditions and critical 
speeds, repeated acceleration runs resulted in 
different displacement values for the sane cri- 
tical speed. In this case, a range of values 
is presented. Each table presents the data for 
a particular configuration. For the baseline 
test using ACM S/N 312 and for Phase I of the 
Design Change Test, the data with respect to 
unbalance is not sufficient for tabulation. For 
the first critical speed, compressor end dis- 
placement is presented; for the second critical 
speed, turbine end displacement is presented. 
Infv/citlon is provided in this manner since 
both analytical predictions (using Che MTI pro- 
gram titled Unbalance hesponse of a Flexible 
Rotor) and baseline test data established that 
the first critical speed is most sensitive to 
response at the compressor end of the rotor sys- 
tem. Similarly, the second critical speed is 
most sensitive to turbine end response. 

After analyzing baseline data, project de- 
sign goals became more visible and well-defined. 
For instance, two of these design goals were to 
minimize shaft excursions resulting from syn- 
chronous rotor response and to move the machine 


second critical speed 151 or more above the 
maximum operating speed. Fulfillment of these 
goals required review of the data to determine 
the following: 

1. The actual relationship between dis- 
placement and a range of unbalance for each 
tested configuration. 

2. The actual relationship between cri- 
tical speeds and a range of unbalance for each 
tested configuration. 

3. An optimum configuration which would 
provide the minimum displacements over the range 
of unbalance for the range of rotor speeds. 

4. The relationship between displacement, 
unbalance and operating time for the optimum 
configuration. 

Displacement and unbalance were expected to 
vary in a direct manner. The manner in which 
critical speeds varied with unbalance could not 
be predicted, however. With these goals and 
expectations in mind, the test results from 
each task were analyzed and evaluated. 

Baseline Test 

Table 4 presents baseline test results with 
respect to changes in unbalance. Changes in 
unbalance caused the most significant changes in 
compressor end displacement. The only signifi- 
cant change in turbine end displacement is 
related to a relatively small change in un- 
balance. There was no discernible critical 
speed relationship with unbalance. 


TABLE 4 
BASELINE DATA 

ROTOR DYNAMIC RESPONSE RESULTS 
ELASTOMERIC O-RING BETWEEN 
CARTRIDGE AND DIFFUSER HOUSING 


UNBALANCE 

FIRST CRITICAL 7 
1 

SECOND CRITICAL 

COMPRESSOR/TURBINE 

COMPRESSOR END 
DISPLACEMENT 

SPEED 

TURBINE END 
DISPLACEMENT 

SPEED 

N-m X 10 6 

3 

mm p/p X 10 

rps 

3 

mm p/p X 10 

rps 

7/7* 

8.4 

250 

38 

467 

14/7.8 

5.1 

267 

76 

500 

42/18 

28 

250 

76 

483 

56/28 

25 

267 

b6 

483 

150/63 

80 

283 

NO DATA 

NO DATA 


♦BASELINE CONDITION 
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Design Change Teat-Phase II 

The results for Phase II are shown In 
Table 5. With each change In unbalance there 
Is a significant change In coapreaaor end dis- 
placeaent be' not In turbine end displacement. 
First critl speeds for various unbalances do 
not exhibit a trend. There Is an apparent 
10 rps (600 rpm) Increase In second critical 
speed for each Increase In unbalance. Also, 
second critical speed Is increased by 45 rps 
(2700 rpa) when ACM unbalance is returned to 
the baseline condition. No analytical rationale 
for these trends is available. 

Design Change Test - Phase II!A 

Phase IIIA results are shown In Table 6. 

The compressor end displacement varies with un- 


balance In the expected direct relationship. 
Turbine end displacement does not Increase with 
Increased unbalance. There are two anomalies in 
the data. The displacement 0.0026 mm p/p (1.04 

mils p/p) fc. the 7.1xl0“*-7.1xl0“ 6 N-m (0.001- 
0.001 or In.) baseline unbalance la greater than 
the 0.017 am p/p (0.69 alls p/p) for the 

slightly worse 7.1 + xlO* 6 -7.1 + xlO~ 6 N-m (0.001 + 

-O.OOl^or in.) baseline unbalance. This anomaly 
consistently appears In the remaining tests. 
Secondly, the greatest unbalance Is related to 
thr smallest turbine end displacement 0.084 am 
p/p (3.31 ails p/p). The only significant 
change In critical speed (a decrease In first 
critical speed to 377 rps (22,600 rpa)) coin- 
cided with the largest increase in displacement. 


TABLE 5 


DESIGN CHANGE TEST - PHASE II 
ROTOR DYNAMIC RESPONSE RESULTS 
LARGER MOUNTING FOOTPRINT AND C-RING - S/N 313 


UNBALANCE 

FIRST CRITICAL 

SECOND CRITICAL 

COMPRESSOR/TURBINE 

COMPRESSOR END 
DISPLACEMENT 

SPEED 

TURBINE END 
DISPLACEMENT 

SPEED 

N-m x 10^ 

3 

am p/p x 10 

rps 

HD p/p x 10 3 

rps 

3. 5/10* 

29 

347 

110 

543 

56/25 

54 

382 

100 

560 

150/56 

92 

363 

110 

570 

3. 5+/ 10+ 

21 

385 

130 

615 


•BASELINE CONDITION 


TABLE 6 

DESIGN CHANGE TEST - PHASE II 
ROTOR DYNAMIC RESPONSE RESULTS 
LARGER MOUNTING FOOTPRINT, C-RINC 
AND INCREASED BEARING PRELOAD - S/N 312 


UNBALANCE 

FIRST 

CRITICAL 

SECOND CRITICAL ] 

COMPRESSOR/TURBINE 

COMPRESSOR END 
DISPLACEMENT 


SPEED 

TURBINE END 
DISPLACEMENT 

SPEED 

N-m x 10 6 

ED p/p X 10 3 


rps 

mm p/p x 10 3 

rps 

7/7* 

26 


460 

100 

568 

32/25 

32 


420 

100 

583 

130/70 

no 


375 

84 

373 

7+/ 7+ 

17 


465 

100 

567 


•BASELINE CONDITION 
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TABLE 7 


DESIGN CHANGE TEST - PHASES IIIB AND IVA 
KOTOR DYNAMIC RESPONSE RESULTS 
C-RING AND INCREASED BEARING PRELOAD 


UNBALANCE 

FIRST CRITICAL j 

SECOND CRITICAL j 


COMPRESSOR END 


TURBINE END 


COMPRESSOR/TURBINE 

DISPLACEMENT 

SPEED 

DISPLACEMENT 

SPEED 

N-m x 10 6 

■a p/p x 10 3 

rps 

mm p/p x 10 3 

rps 

PHASE IIIB - S/N 313 I 

7/7* 

14 

342 

97 

638 

42/28 

17 

305 

110 

653 

140/92 

53 

300 

150 

675 

7+/7+ 

9.4 

345 

110 

652 

| PHASE IVA - S/N 312 j 

7/7* 

58 

412 

62 

575 

25/14 

41 

<35 

84 

595 

46/21 

75 

392 

110 

602 

7+/ 7+ 

41 

432 

* 

587 


♦BASELINE CONDITION 


Design Change Test - Phases IIIB and IVA 

Two phases of design test. Phases IIIB and 
IVA, were run with the same configuration. The 
data for these phases are presented In Table 7. 
There Is an anomalous point In 'he data for 
Phase IVA - the compressor end displacement for 

the 4.6x10 3 -2.1xl0 3 N-m (0.0065-0.003 or In.) 
unbalance condition. Except for this point, 
both compressor end and turbine end displace- 
ments are related to unbalance In the typically 
expected direct manner. The Phase IIIB results 
showed a high second critical speed (and a 
relatively low compressor end displacement at 
the first critical speed). Even though a criti- 
cal speed - unbalance relationship could not be 
ascertained, the Phase IIIB configuration was 
selected as a potentially acceptable solution 
since a majority of the project goals (and 
the aost of any configuration tested thus far) 
were aet. This acceptability became more ques- 
tionable when the Phase IVA results did not show 
complete similarity with the Phase IIIB results. 

Design Change Test - Phase IVB 

The Phase IVB configuration was tested to 
determine the effect of the 0-rlng on the base- 
line test results obtained from ACM S/N 3l2 and 
S/N 313. This configuration demonstrated sur- 
prising results with nearly constant low levels 
of rotor displacement regardless of unbalance. 

In addition, critical speeds appeared to have 
no more than a small random variation with 
unbalance. These findings were sufficient to 


result in the selection of the Phase IV 
configuration for the reliability test. 

Reliability Test 

Data from the last design change test and 
from the reliability test are summarized In 
Table 8. For Design Change Test, Phase IVB, 
there Is no significant increase in displace- 
ment. Both compressor end and turbine end 
displacement are not sensitive to unbalance. 

The 3.5xl0'^-1.4xl0'^N-m (0.0005-0.002 oz 
in.) and 3.5x10 -2. 8xl0 _:> N-m (0.0C 05-0. 004 oz 

in. ) reliability test unbalance were Imposed on 
the machine as all unbalances for previous tests 
were imposed on ACMs. The machine was then 
returned to a baseline unbalance. For the series 
of tests thereafter, the unbalance was not 
deliberately changed; the ACM rotor was permitted 
to become unbalanced during operation. Measure- 
ments made after the last test showed that the_^ 
ACM did, become significantly unbalanced 2.1x10 . 
1.5x10 N-m (0.003-0.021 oz in.). The 1.4xl0 -;> 

-6.3x10 3 N-m (0.C02-0.009 oz in.) unbalance is 
an estimate. 

The imposition of unbalance did not Increase 
displacement at either end of the rotor. When 
the ACM was returned to a baseline unbalance, the 
displacement at both ends decreased; the turbine 
end displacement decreased significantly. When 
the machine was permitted to become unbalanced, 
test data showed no Increase in displacement at 
either end of the shaft. 
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TABLE 8 


DESIGN CHANGE TEST - PHASE IVB AND RELIABILITY TEST 
ROTOR DYNAMIC RESPONSE RESULTS 
INCREASED BEARING PRELOAD AND CLEARANCE 
BETWEEN CARTRIDGE AND DIFFUSER HOUSING - S/N 313 


UNBALANCE 

FIRST CRITICAL 

SECOND CRITICAL j 


COMPRESSOR END 


TURBINE END 


COMPRESSOR/TURBINE 

DISPLACEMENT 

SPEED 

DISPLACEMENT 

SPEED** 

N-m x 10* 1 

3 

mm p/p x 10 

rp6 

3 

am p/p x 10 

rps 

DESIGN CHANGE TEST - PHASE IVB 

7/7* 

17 

305 

47 

600 

14/11 

20 

333 

45 

600 

39/21 

22 

283 

47 

605 

7+/ 7+ 

15 

335 

49 

605 

RELIABILITY TEST 

7/7* 

30-34 

283 

43 

550 

3.5/14 

34 

312 

37-39 

525 

3.5/28 

28 

337 

34-41 

517 

7/7 

18 

275 

8.7 

550 

14/63 (EST) 

23 

257 

10-11 

550 

21/150 

19-23 

260 

6.4 

550 


‘BASELINE CONDITIONS 

“Displacement was increasing at upper limit of test run; critical speed was 
estimated to exist beyond this speed. 


There was no discernible trend in firs: 
critical speed. Displacement increased at the 
upper limit of all but one of the acceleration 
runs. The second critical speed could not be 
accurately determined. Figures 6 and 7 are 


plots of the compressor end and turbine end dis- 
placements respectively. Each figure includes 
one curve of displacement versus rotor rps for 
each unbalance condition. 


rm»r mf 



Fig. 6 - Reliability Test Rotor Dynamics Data 
from Real Time Analyzer Processing 
Methods. Compressor End Activity 
(Probe Y1Y ) 



Fig. 7 - Reliability Test Rotor Dynamics Data 
from Real Time Analyzer Processing 
Methods. Turbine End Activity 
(Probe Y5Y) 
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Review of these figures begins with a com- 
parison of the 335-550 rps (20,000-33,000 rpm) 
portions of curves IA and II of Figure 6. This 
comparison shows the only data for which Increase 
in unbalance did not actually decrease displace- 
ment. The remaining data of Figures 6 and 7 
(curves 1A, II and III) showed that increased 
unbalance reduced displacement at both ends of 
the shaft. When the machine was allowed to 
become unbalanced during operation, the data 
(curves IV, V and VI of Figures 6 and 7) showed 
no significant increase in displacement at 
either end of the shaft. These results con- 
firmed the selection of the Design Change Test 
Phase IVB configuration. 


SUMMARY AND CONCLUSIONS 


An effort was made to determine the rela- 
tionship of the rotor dynamic response to 
unbalance. In each of the tests, the critical 
speeds (synchronous whirl conditions) varied 
randomly. Similarly, no specific relationship 
could be determined between non-synchronous 
whirl conditions and unbalance. During Phases 
II and IIIA of the Design Change Test (as in the 
Baseline Test), there was a direct relationship 
between unbalance and compressor end displace- 
ment due to first critical speed. There was no 
discernible relationship between unbalance and 
turbine end displacement due to second critical. 
During Phases IIIB and IVA, both of the dis- 
placements were directly related to unbalance. 
For the Phase IVB test, there was no discernible 
relationship between either of the displacements 
and unbalance. The second critical speed was 
relatively high and was not significantly 
altered by changes in unbalance. Although U' - 
balance relationships were difficult to assess, 
the Phase IVB configuration, by providing the 
most significant reduction in rotor displacement 
through a typical range of unbalance, was 
selected for Reliability Test. 

The first part of the Reliability Test 
involved unbalance Imposed on the machine. Con- 
centrating only on the compressor end and tur- 
bine end displacement at first and second criti- 
cal speed respectitely, there was no discernible 
relationship with unbalance. Considering the 
entire speed ranxe, displacements at both ends 
of the shaft decreased when unbalance was 
increased, After the machine was returned to a 
baseline unbalance condition the rotor operated 
without any controlled, imposed unbalance. 
According to the predictions, rotor unbalance 
Increased. However, contrary to expectations, 
shaft displacements and critical speeds remained 
nearly unchanged. Further analytical efforts 
are required in order to determine the rationale 
for these unusual trends. 
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Introduction 


The cos Angeles Division of Rockwell Interna- 
tiona.’. Corporation has been engaged since the mid 
1960 ’s in development of the manned penetrating 
attack aircraft, the B-l. A research, develop- 
ment, test, and evaluation program initiated 
during May 1970 continues. Currently, three of 
the four test aircraft are on flight status. 


Structural acoustics considerations have played 
major roles, perhaps to an unprecedented extent, 
in establishment of the design conditions for this 
aircraft. Because the engines are mounted in 
nacelles under the .ling roots IFigure 1) and are 
large duct-afte'buming turbofans, the majority of 
aft fuselage and empennage structure was sized for 
acoustical fatigue prevention. The vehicle has 
three l 5- foot- long, deep weapons bays, required to 
withstand cavity pressure-oscillations during 
weapons drops throughout the flight envelope. Out 
lines of the three sets of weapons bay doors can 
be seen in Figure 2. The aft structure and the 
weapons bays house eauipment in potential ly severe 
dynamic environments. 


The elaborate and critical dependence of the 
stiuctural design of the airframe on acoustical 
fatigue requirements has provided severe tests of 
the availability of data bases, of laboratory 
techniques for adding to them, of methods for 
specifying requirements and combining them with 
those for other load conditions, and of techniques 
for their detailed implementation in the design. 

This paper considers a selected subset of the 
structural acoustics work performed in support of 
the B-l design, the results of which may have 
utility for future programs. Specifically, some 
results of measurements of engine noise on the 
airframe surface during static ground runs are 
presented, and the treatment of those results in 
preparation for a structural redesign is discus- 
sed, Wi tii regard to design for acoustical 
fatigue, the somewhat unusual approaches to the 
acquisition in laboratory tests of design data 
for two configuration concepts are reveiwod. Some 
coasiderations related to mechanical fasteners in 
acoustical fatigue applications are explored. 
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Mode 1 for Acoustic Loads 


Figure i shows the aircraft in flight, with 
the variable -diameter nozzles of the General 
Electric F101, 30,000-pounU-class engines open to 
the afterburner position. The figure and Fig- 
ures 2 and 3 suggest the spacial relationship 
between the aft structure of the airframe and the 
iioise field of the engines. Engine noise ranged 
from contributory to dominant as a design condi- 
tion for the fuselage aft of the engine nozzle- 
plane and for the empennage. Because of the 
attendant weight and the cost of deta’ls appro- 
priate for fatigue resistance, and because an 
opportunity to fine-tune component thicknesses 
would exist for much of the structure for A/C-4, 
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an elaborate static sound survey for these sur- 
faces was conducted on A/C-l. Measurements were 
made at 280 locations on the structure for sev ral 
combinations of number of engines and engine 
thrust-setting. Figure 3 shows the aft fuselage 
with the microphones mounted on pads of closed- 
cell foam. The diaphragms face aft and are about 
2 inches from the surface. 

The results for four-engine, niaximum after- 
burner ope rat on wete used in construction of a 
model for acoustic pressure on the structure. 

Some of the steps taken and results produced in 
tha r process, as it dealt with the aft fuel bay 
and aft equipment bay, follow. 


Overall pressure levels, raised 1 decibel (dB) 
to correct to sea level from the elevated measure- 
ment site, Edwards Air Force Base, and 2 dB for 
thrust loss related to the ambient temperature, 
were plotted for eight water and (roughly) buttock 
lines (Figure 4). The amount of data thins beyond 
station 1600, as fewer than eight microphones pro- 
vided coverage on a roughly 20- by 20- inch grid at 
station planes having smaller cross sections. 

Some attrition occurred; its effect is particu- 
larly evident at station 1640. 
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The resulting smoothed values were cross- 
plotted for station planes at 50-inch intervals 
(Figure S) . Interactive smoothing between the two 
sets of curves provided the results shown here. 
These levels then were transferred to a fuselage 
diagram, and smooth contours were drawn. Fig- 
ure 6 is the result. Figure 7 shows the overall- 
pressure- level field on the lower surface of the 
horizontal stabilizer; the contours were devel- 
oped in a similar process. 
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Frequency Distribution 

It frequently is assumed for engine-noise pre- 
dictions that spectra at any station plane have a 
shape roughly independent of distance above, 
below, or to the side of the wake, provided those 
distances are not excessive and that special cir- 
cumstances, such as surface or ground reflection, 
do not dominate. (See Reference 1, for ex amp le . I 
Results of work with this large data base support 
this concept. 


Overall Pressure Levels 
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Samples from the 15-second records were 
reduced in one-third octaves from 10 to 8,000 
hertz (Hzj. Figure 8 presents live (from a total 
of 17) third- octave plots for roughly equal spaces 
along the fuselage. Each of the values shown 
(as dB below the overall level) is the average 
taken from the signals of all microphones at the 
indicated fuselage station. The characteristics 
and suitability of these averages are discussed 
in the following paragraphs. The distances from 
the wake to the surface vary from about 10 to more 
tli an 100 inches. The closest microphones to the 
ground (those forward, at the bottom centerline) 
were 110 inches above it. 

These average third-octave- level distributions 
(involving, of course, logarithms of pressure and 
frequency) seem remarkable in that they can be 
approximated tightly in each case by three 
straight lines (although a gentle S-shaped curve 
would yield a slightly tighter fit for reg.on 
III). The lines of Figure 8 are least-square- 
error fits. 


A model of frequency distribution was con- 
structed for use, with the overall contours of 
Figure 6, in tne A/C-4 design process. For this 
purpose, each of the 17 station-plane one- third- 
octave distributions was represented by three 
least-square-error lines. Figure 9 presents 
smoothed longitudinal fits to values calculated 
from the linear representations for several third- 
octave bands. The curves for the 200 and 630 Hz 
bands, which define the region K range of the 
model, approximate the points to within ±1 dB. 
Those for the 125 and 1,250 Hz bands, which are 
from the high-pressure portions of regions I and 
III, respectively, also fit the values from tbs 
linear equations tightly. The curves for the 
31. S and 2,500 Hz bands, roughly the middle bands 
for regions I and III, respectively, manage a 
close fit at the cost of some distortions in 
their shapes. This degeneration of shape at 
frequencies remote from region II continues 
to the curves for the extreme 10 and 8,000 Hz 
bands. 
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The curves of Figure 9 permit construction of 
a three- line representation for design purposes 
that surely provides a better estimate of the 
frequency distribution at an arbritrary point on 
the surface than would interpolations based on 
measurement results nearby microphone locations. 
This can be asserted .v.th particular confidence 
for frequencies in and near region II, the range 
of interest for acoustical fatigue analysis, as 
will be shown. 

Several observations can be made about the 
representations of Figure 8 and the 12 others not 
shewn: 

1. The five shown are typical of the 17. 

2. The individual band levels "add" to zero 
dB (to within one-half dB) , as is required by the 
meaning of "overall." 

3. The three- line representation is superior 

in fit to that provided by any single gentle curve. 

4. The efficacy of the three-line fit is not a 
result of the averaging process that leads to the 
band averages; the same three- segment property is 
found for the distributions of the single-micro- 
phone signals, although, of course, each has more 
scatter than the average. 

5. The trace of the frequencies of intersec- 
tion between the lines representing regions I and 
II rises to 200 Hz at about station 147S and re- 
mains roughly constant aft of there. The corre- 
sponding trace of intersections between the lines 
representing regions II and III rises to a peak of 
about 1,100 Hz in the vicinity of station 1400 and 
declines smoothly aft of there to about 500 Hz; at 
station 1680, the region II and III lines effec- 
tively have merged. 

More needs to be said about the primary in- 
gredient of the frequency-distribution model; i.e., 
its treatment of the distribution at a station 
plane as an invariant, and its representation for 
that purpose by the band-by-band averages of 
levels (re the overalls) from all microphones at 
that station plane. It is not obvious on physical 
grounds that this should lead to an acceptable 
representation. Figure 10 presents the results of 
an evaluation of adequacy: 

1. Consider, for example, the point plotted in 
Figure 8 for the 10 Hz band for fuselage station 
1342. That value is the simple mean of the levels 
(in dB below the corresponding overall levels) for 
the signals from seven microphones spaced around 
the fuselage, from the top of the fuel bay to its 
lottom centerline, at that station plane. 

2. The standard deviation of those seven 
levels about the plotted mean was calculated, and 
from it the confidence interval at the 90-percent 
level was computed. Similar confidence intervals 
for the 10 Hz band for the other four stations 
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represented in Figure 8 were compute’., and the 
average for the five was fort _d. gurv 10 pre- 
sents that average and those for he other 29 
bands. Here, the averages are weighted for the 
number of contributing signals at each station 
(i.e., seven or eight). The figure shows a series 
of averaged inter a. of roughly ±1 dB from 63 to 
1,250 Hz; this range includes the frequencies of 
importance for acoustical fatigue performance of 
the fuel and equipment bays. 

3. Also plotted in Figure 10 are thp confi- 
dence intervals at the yo-percent 2»vel resulting 
n-om an e v actly parallel wi igh ted- overage treat- 
ment for five groups of sewn or eight third- 
octave frequency distributions for a "white no- =e" 
(produced '<y a signal generator) processed by the 
same data-reduction program. Th: comparison 
demonstrates that the uncertainties shown for the 
pressure model are attributable to the treatment 
of the frequency distributions as invariant af a 
station plane, not to the data sampling and reduc- 
tion processes. 

Acoustical Fatigue Design 

The high intensities of pressure sources, 
large areas of involvement, and multiplicity of 
incorporated structural concepts make the B-l a 
candidate for history's most complex aircraft 
design with respect to acoustical fatigue. None- 
theless, the practice of this discipline had 
achieved suificieet maturity that, with very few 
exceptions, concepts, thicknesses, and detail 
requirements could be established on the basis 
of analytical modifications to existing, largely 
published information. 

Some of the concepts in this eutc-goiy of 
requirements specification by analysis ineorpo 
rated into the design are skin and r rame (or skin 
and rib), flat anu curved; honeycomb, full -depth 
and panels, flat and curved (aluminum, titanium, 
fiberglass and graph it e/epoxv ) ; beaded or corm- 
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gated inner-skin panels (aluminum and diffusion- 
bonded/superplastic- formed t tanium); quartz fiber 
glass panels of various constructions; ported 
plates on beam suppoits: large access panels 
supported by trusses; cantilevered access panels; 
and integrally stiffened skins (aluminum and 
titaniiml . 


The effect of the straps is to replace, in 
part, the inertial reaction to the rocking 
motion of the frames by tensile forces in the 
straps. This reduces skin deflections and 
stresses through the partial restraint of frame 
motion, but at the cost of increased web load 
and of concentrated loads at the strap- to- inner 
cap intersections and in the web stiffeners 
adjacent to those intersections. The sketches 
of Figure 12 suggest these consequences. Fig- 
ure 13 shows six straps (horizontal in the 
photograph) riveted to the inner caps. 


Two important exceptions to design based on 
analysis arose. The two configurations wer" to be 
used over large areas of the airframe subjected to 
high acoustic load and, thus, had potential fer 
costly and disruptive maintenance requirements. 
While these considerations entered into, they did 
not govern the decision to require tests. What 
did prove determinant was that each represented a 
significant and not readily analyzable departure, 
in ways important to dynamic behavior, from struc- 
tural concepts for which performance data existed. 
These exceptions are considered in the following 
paragraphs. The test programs for them exhibit 
seme unusual approaches to design support and 
validation. 


One of the configurations for which test 
results were required is used for the aft fuel 
bay. (See Figure 6 for its location.} This 
20-foot-long section is constructed of skins 
supported by lightweight, built-up frames on 
8.3-inch centers. The stricture forms the 
tank surface and carries the structural loads 
in this area. The bay is essentially enptv of 
fuel during ground roll and takeoff, the time 
during which it experiences significant acoustic 
load; it is filled to rebalance the aircraft as 
the wing is swept aft. The configuration differs 
from conventional skin -nd frame structure in 
that it has longitudinal frame-to- frame straps 
that react fuel-slosh loads and radial frame-web 
stiffeners that divide the webs into small shear 
panels. Figure 11 is a sketch of one of the 
frames, that used for the test specimens. Loca- 
tions of the slo«h str-ps and web stiffeners are 
indicated. 
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The effect of the stiffeners is to compel 
reaction to the rolling motion of the outboard 
fra'ie caps, induced by skin deflections that are 
uns> iimc-tric about the frame, to occur primarily 


at the web stiffeners rather than continuously 
along the web. fhe result is locally concen- 
trated toad in the outboard frame cap, in the 
stiffener at its intersection with the cap, and 
in the fasteners that tie them together, fig- 
ure 14 shows a frame segment with some of its 
stiffeners installed. 
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The Structural Model . The fuel -hay design for 
the first three test aircraft was based on requi ce- 
ments generated by a structural model; i.e. , a 
series of equations that accounts for its vide 
variety of curvatures, temperatures, and sound 
pressures. The coefficients of the model, which 
Jeterained component thicknesses and iasterer re- 
quirements, were based on considered est invites, 
and the details were executed with the informed 
judgment used for configurations not subject to 
test, [Tha judgment include. I the conviction Taut 
the stiffener details shown in figure 14, ini- 
tially strongly desired for static strength pur- 
poses, could not prove successful.) Requirements 
generated by the model made clear that acoustical 
fatigue miu K 1 dominate the design of the fuel hoy. 
In view of the uncertainty about the effects of 
tiro slosh straps and web stiffeners on the per- 
formance of the configuration, an absence of test 
results for use in airworthiness evaluations of the 
three test aircraft and in the design revisions for 
•VC-4 would have compelled use of coefficients 
that were substantially more const votive than 
those chosen. An unnecessarily heavy and not 
assuredly reliable structure would have resulted, 
Opf inti tat ion through and element siring by test 
was essential. 


1. The tests evaluated the structural model, 
used for design of the airframe and the specimens, 
rather than the fuel bay itself. 

1 . Opt i.ni tat ion was sought through an evolu- 
tionary technique of developmental testing, in 
whidi deliberately underst rerig th substructure 
components and mechanical fas toners were installed 
initially and upgraded as the test revealed rela- 
tive fatigue strengths. In this process, evalua- 
tion and refinement c.f design details also took 
place. 

5. As a necessary concomitant to the evolu- 
tionary method, "step testing," the technique in 
which testing proceeds in steps of time at pro- 
gress ively higher sound levels, was employed. 

The nrthod was developed for use with sinusoidal 
sources; in this instance, a random- in- t ime 
source was used. 

It is the purpose of this discussion about the 
fuel bay to exhibit the reasoning that led to these 
approaches and the techniques of their execution. 
The equations of the model (with the exception of 
that shown in Figure 27, which illustrates the 
genre) are not presented; without sketches of the 
components and details, presentation of which is 
beyond the scope of this paper, they would not be 
useful. 

The experimental design for the fuel-bay tests 
evolved concurrently with development of the 
structural model. Because of the broad ranges of 
acoustic load, curvature, and surface temperature 
in the airframe application and the need for evolu- 
tion of a configuration optimised for minimum 
weight, it was clear that a developmental test of 
the model (with specimens that incorporated all 
pertinent details), rather than a proof test of 
the hay itself or a portion of it, would be appro- 
priate. Figure 15 indicates the relationships 
among the experimental design, the structural 
model, and the design of the airframe. 
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The laboratory tests performer! for this purpose 
were unusual in concept in three respects: 


Fval nation of t he Structural Model , The 
evaluation of the mrslel concerned two assumptions 
about how the configuration functions; once these 


r>o 







were confirmed, attention could shift to the 
refinement of constants for thickness requirements 
and to the performance oi details. One of the 
assumptions is that an exist mg hut not firmly 
established correction for curvature to skin 
thickness and, in turn, to the thicknesses for 
substructure elements tall of which are sired in 
proportion to the skin) is valid over the wide 
range of curvatures of the fuel bay. this meant 
that the specimens had to provide independent 
tests for a variety of curvatures. The other 
assumption is that all substructure elements other 
than fasteners are critical in bending rather than 
in tension or transverse shear. \ description of 
the specimens and a statement of the criterion 
for demonstration of the validity of the model 
follow. 

The cross section at fuselage station li'5 was 
selected for the test specimens because of its 
wide range of curvatures. iSee Figure lt». ! Die 
specimen of Figure 15 is from the bottom of this 
section; it is symmetrical about the bottom 
centerline of the fuselage. It has a radius of 
curvature of t-t> inches. Figure 1' shews the 
side specimen mounted in the test window; it 
represents half of the fuselage, from upper to 
lower centerline. In that orientation, a region 
having a 69- inch radius of curvature was exposed 
to acoustic load. A second orientation exposed 
the portion between this and the bottom region; 
it has a 5t> inch radius of curvature. A longeron 
divides the two regions. Parts of a heating sys- 
tem used to maintain the skin at 509 ^ F during 
the test are evident in Figure 15. The specimens 
comprise six duplicate frames (those to be tested! 
and two stronger ones at each end to terminate the 
specimen and to provide shear -panel terminations 
for the slosh straps. (See Figure IS.) 
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IVmonstrat ion of the validity of the model re- 
quired that any element type of the skin or sub 
structure that failed in any region would begin to 
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experience essentially identical failures at 
roughly the same load/ time conJuti.it ion in all 
regions, i.e., independent of the curvature of 
the skin. The 1 1 -~* t results did show uniformity 
among the regions in load level lor, and reasuiv 
(tl’le cuenpat ability in time to, failure for each 
element type. 

Opt mi; at ton through Ivolntien . The web stift 
oners and slosh straps create redundant load paths, 
thus rendering analysis of thickness intenvl.it ion 
ships for the structural model Jtvidedlv approx 
imate. The refinement of the-e relationships 
through developmental testing followed this 
sequence : 
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1. The skins and outer frame caps, made of a 
titanium alloy in the original design, were sized 
by an analytical estimate. The acceptance of a 
fixed thickness relationship between these two (as 
a function of frame spacing) was a concession to 
the long lead times for material acquisition and 
fabrication for the specimens, a characteristic of 
parts made from this material. The choice was 
hedged by simultaneous procurement of additional 
specimens that differed from the first only in 
having outer caps of increased thickness. 

2. The remaining subs . ructure components, all 
of which are of alumimm alloy, were sized to be 
as balanced as analytical limitations would per- 
mit, but were hedged to be inferior to the skin 

3. A similar inferiority of fasteners was 
sought. This was achieved through substitutions 
for the airframe fasteners by selection among 
size, material, spacing, and type. 

4. The test for each region was initiated at 
an acoustic level substantially below the design 
objective for the specimen; this load was imposed 
for a major portion of the d 'sign lifetime. The 
process was repeated for pressures increased in 
3-dB increments. 

5. Repairs or replacements were made as compo- 
nents and fasteners exhibited damage. 

6. Once patterns of damage occurrence had 
been established for each conponent and fastener 
application in each of the three regions defined 
by curvature, a resized and redetailed substruc- 
ture that incorporated the rough optimization 
implied by the damage interrelationships was sub- 
stituted fer the six test frames and the straps 
in both the bottom and side specimens. (See 
Figures 14 and 18.) 

7. Testing was continued, supported by addi- 
tional repairs, until a full lifetime of exposure 
at the maximum level capability of the facility 
had been imposed for each of the three regions. 

The near-optimum relationships, absolute thick- 
ness requirements, and refinements of details 
established from the results of this procedure 
were used in a revision of the structural model 
that governed the A/C -4 redesign. The redesign, 
of course, made use of the measurement-based loads 
of Figures 6 and 9. 

In an unanticipated result, surface- tenperature 
tests on A/C-l established that aluninun (rather 
than titanium) alloy skins and outer caps could be 
used for the fuel bay. The A/C- 4 design incorpo- 
rates this change at a substantial cost saving and 
slight weight penalty. The substitution was sized 
for acoustical fatigue by a routine analytical 
adjustment of the test results for allowable 
fatigue stress. Because of the certainty that 
some loss of optimization resulted from the inabil- 
ity to establish a permissible stress level in the 
aluminum skins and caps without allowing increases 
in their stiffnesses over those of the titanium 


Conponents they replaced, it was decided uiat 
proof tests of the revised design would be con- 
ducted. The bottom specimen has been tested at 
design load plus a 3-dB margin. After about 
10 hours of the 26-hour objective, a slosh strap 
failed. The test has not resumed, and the side 
specimen has not been tested. 

Step Testing. The step- test method works well 
with sinusoidal sources. A full life ime of expo- 
sure to stresses at levels that do damage in a 
random- load application can be applied for a 
single mode in a sinusoidal- load test in about 
15 minutes. Seldom are more than four modes 
recognized as potential contributors to signifi- 
cant damage. Thus, loads in a series of several 
steps, with interruptions for periodic inspections , 
can be applied in a few days. (See Reference 2, 
for exanple.) However, with randan- in- time sources, 
the use of which constitutes a major inprovement ■ 
in testing realism, the accumulation of damage 
takes place in real time. A lifetime of exposure 
to the loads of ground runs at high thrust and 
takeoffs for the B-l, for exanple, is 26 hours. 

Only an inpe-rative that transcends time and cost, 
in til is situation the revelation of balance among 
regions with different curvatures and the quest 
for optimization, can justify multiple load steps 
of this duration in a large, costly to use 
facility. 

Stabilizer Structural Boxes 

The other configuration for which evaluation in 
laboratory tests was required is that used for the 
main structural boxes of the horizontal and verti- 
cal stabilizers. This structure has thick aluni- 
num skins to provide stiffness foT flutter resis- 
tance, multiple lightweight spars of welded tita- 
nium sine-wave-web and cap construction, and light- 
weight alviunun sheet-metal ribs. Figure 19 is a 
structures diagram for the horizontal stabilizer, 
and Figure 20 displays the details of the sine 
wave spar. Minimum thickness on th- ^rframe for 
manufacturing operations and handlir 0.020 inch 

for the cap and web; the latter is i 1 in two 

steps by chemical milling to a minimun. f 0.012 
inch. With few exceptions, these minimum thick- 
nesses are adequate for static strength requirements. 
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1. The program involved a dccid-xily unusual 
example of simultaneous tests of several boxes in 
a large facility. These were made the elements of 
two larger specimens, so that the ribs that formed 
interfaces between the boxes could be evaluated 
along with the spars. The virtue of economy was 
co-nilemented, in Fact overtaken, by necessity. 

Hach specimen included individual boxes of 
widely varying strengths Therefore, as for the 
fuel bay tests, stepped loading was used. Mechan- 
ical neutralization of damaged boxes permitted 
continuation of the test. 

3. rhis program of limited and not obviously 
adequate scope provided all the information neces- 
sary to support qualitative design requirements 
for the airframe. 


((full 20. O'tMli •' '« wokfod titanium um wm ip« Tbo mottnol n •■■oolod 
SA1-4V alloy. 

In acoustical fatigue, the significant motion 
for this configuration is that in which adjacent 
skin panels move out of phase with one another, 
attempt irg to roll the spar cap and how the web. 

( See figure ’1.1 ITie sine wave web provides stiff, 
broad-based support to the cap: in spite of the 
thin gauges, the rolling mot ion is resisted effec- 
tively. Nonetheless, because attempts at static 
stress analysis of the spar for this motion were 
not rewarding, no strong preconception ol acousti- 
cal fatigue adequacy could he sustained; it was 
necessary that tests be performed (.11 to provide 
an evaluation of the airworthiness of the stabili- 
zers of the first three test aircraft and (2) to 
supply quantitative requirements for the A/ 1 ' - 4 
redesign or, alternatively, to establish what 
additional vests would be necessary for the latter 
purpose. 
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Hie test plan developi.il for ibis purpose is 
noteworthy in three respects: 


Fxperimer.tal Design . Figure 22 illustrates the 
scheme from which the coiqnonent combinations were 
selected. Results from tests whose component 
sizes are governed by a fully- filled cube of 27 
combinations (to account fer three values each of 
three independent variables) would be required in 
principle to establish the functional relation- 
ships among allowable acoustic load (for one life- 
time, for example) and the spar- related indepen- 
dent variables: spar spacing, and cap and web 

thicknesses. Limitation to even this large number 
(.1) assumes scatter sufficiently small that the 
relationships are not obscured, (2) does not pro- 
vide for the required information regarding rib 
and fastener damage, and (3) requires tliat damage 
not related to the spars does not intervene. In 
these tests exogenous variables, un shunned gaps 
between the rigid skins and other components, 
assured wide scatter. Spar damage traceable to 
prior rib failures, a combination referred to 
below as the rib/ type- IV syndrome, occurred; it 
proved to be of primary importance for airworthi- 
ness evaluations of the stabilisers. 
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Iii recognition of the uncertainties about per- 
formance, of part and test costs, and of time 
limitations, it was decided that a preliminary 
series of tests would be conducted to evaluate the 
18 combinations of Figure 22. A limited objective 
of the establishment of minimum capabilities for 
ihe spar components and of functional relation- 
ships for rib thicknesses and fasteners was inher- 
ent in this choice. 



Test Spec linens . Figure 23 is a sketch of one 
of the two specimens. This specimen conprises 
three major beams (simulations of the front and 
rear beams of the airframe) , the six boxes of 
Figure 22 that have 6.8 inch spar spacing, and 
eight sets of ribs. Because the rib elements at 
the open ends of the specimen could be replaced 
easily and these inside only with great diffi- 
culty, the former were sited initially at 0.02£ 
inch at one end and 0.032 inch at the other. One 
line of inner ribs was set at 0.040 inch, and the 
other at 0.050 inch. The other specimen con- 
tains 12 boxes, four with 3.0- and eight with 
3.8- inch spar spacings; it has four major beams, 
five lines of ribs, and similar rib-thickness 
arrangements (Figure 24). The maximum aspect 
ratio for individual skin panels for each of the 
two specimens is 3.5 li.e., effectively infinite). 
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Test Procodin Table 1 lists the levels and 
durations for the step- test procedure. Figure 24 
shows the second specimen after substantial damage 
had occurred. Eight of the nine internal ribs had 
failed; external supports had been applied over 
them on the upper and lower skins. The end Tibs 
had been upgraded. Some spar damage had been 
detected at the time the photograph was taken; for 
that reason, stiffeners were to be applied over 
the damaged boxes themselves. 
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Test Results. The results from the tests of 
these two specimens provided all information nec- 
essary fer the redesign for A/C-4. Table 2 sun- 
marizer the types of damage. It can be inter- 
preted more readily with the help of Figures 20 
and 25. 


TaMa l TypaasM 


Skins - no damage 


Front and rear beams - no damage 


Spars - damage in four modes: 

Type 1 - scalloped cracks of caps 

Type II - web- land cracks adjacent to cap 
welds 

Type III - web- field cracks* 

Type IV - web splice-plate cracks adjacent 
to cap welds 

Ribs - cracks in rib segments and clips 

Fasteners - losses in five applications: 

2219-TS1 regular head rivets through lower 
skin into spars and ribs 
2219-T42 (Inconel 600 stem) Cherry bulbs 
through upper skin into spars and ribs 
2219-T81 protruding head rivets in 
substructure 
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1. Type I spar damage is a primary failure. 
Prevention of its occurrence on the airframe 
requires caps of more than minimus thickness in 
only a few areas of combined high zcoustic load 
and wide spar spacing. A number of incidents of 
type I damage occurred in the tests as a con- 
sequence of unshimned gaps between the caps and 
the skins; fastener- induced preload in the rap 
results. A shim requirement that limits the mis- 
match to 0.010 inch for the airframe precludes 
damage attributable to this condition. 

2. Type IV damage is derived, or secondary. 

Only when the rib segments at both sides of a web 
splice-plate have failed is the plate free to move 
enough to fail. The failure is on the plate at 
the plate-to-cap weld; thus, it is similar to the 
type II failure on the web proper. Type IV fail- 
ures in the test soon were followed by secondary 
failures of the adjacent webs in either the type II 
or III mode. 

3. A few cases of primary types tl and III web 
damage (i.e., damage not incident to the rib/ 
type- IV syndrome) occurred. These were used to set 
minimum thicknesses for the webs. A slight in- 
crease beyond the minimum thickness for fabrication 
over a small area of the horizontal stabilizer 
resulted. 

4. Fortunately, both on their own account and 
in view of their role in the type IV spar-damage 
syndrome, a great deal of information about thick- 
nesses for the rib segments was obtained. Fig- 
ure 26 presents these results. Here, "primary" 
refers to the first failure of any segment in one 
of the six-segment rius. "First secondary" refers 
to the first occurrence of failures of an adjacent 
pair of segments, the damage mode essential to the 
type IV syndrome. The first secondary failures 
collapse reasonably well to the function noted 

in the figure; the primar failures do not. The 
dashed line and "design line" have the slope of 
that function. Note the weak dependence of allow- 
able pr jure on spar spacing. For the A/C-4 
. design, the acoustical fatigue requirement exceeds 
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that for strength for approximately half of the rib 
segments of the horizontal stabilizer and for a 
few of the vertical stabilizer. 
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S. A substantial number of fastener failures, 
in all five classes of application noted in 
Table 2, yielded adequate information for design. 
For A/C-4, roughly half of the fasteners in 
the horizontal stabilizer and some in the verti- 
cal stabilizer are governed by the results of 
these tests. 

The wide variety of modes of failure, primary 
and derived, assured that only limited information 
for requirements could be obtained for the caps, 
web lands, and web fields. The design values for 
these conqjonents were established from a combina- 
tion of primary failures and undamaged parts. 

They are, therefore, somewhat conservative over 
the range of spar spacings represented in the test. 
That these conservative requirements led to only a 
small weight penalty for A/C-4 meant that further 
tests were unnecessary. 

Design Requirements . The design values for 
A/C-4 established on the basis of these test 
results form a structural model for the configura- 
tion, equally applicable to the horizontal and 
vertical stabilizers. It is a series of 15 equa- 
tions, the subjects of which are the components 
and fasteners listed in Table 2. IWo examples 
from that description are considered nerer 

1. The thicknesses of bent- flange sheet -metal 
rib segments, in inches, are controlled by: 

Log 10 t R » dB^/27.6 + D. 143 !.og 10 V - 6.839, 


CD 


3.0 £ K £ 8.0 


where W is the spar spacing in inches. Ihis 
equation corresponds to the "design line" of Fig- 
ure 26. The safe range of applicability is 
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limited to small extrapolations beyond the 3.8- 
and 6.8-inch exti ernes of spar spacing represented 
in the specimens. The design pressure is stated 
as a one- third-octave level; as Figures 7 and 19 
show, it need not be constant over the length of 
a rib. 

2. Once a requirement for rib thickness has 
been established, it must be satisfied for the 
equations that describe the fasteners that attach 
the rib to other components. The limiting rib 
thickness for which a particular blind fastener 
of diameter d is applicable for joining the skin 
and ribs is set by: 

t R s (3.35 x 10' 1 x Fp/F R x l/s) 1/2 x d (2) 

where s is the rivet spacing and 1 is the lever 
arm defined in Figure 29. Any consistent set of 
units is acceptable. The factor F/F is a dimen- 
sionless ratio of figures of merit in fatigue for 
the fastener and the rib. The ratio is near unity 
if the fastener and coaponent are of alloys of the 
same material, e.g., aluminum or titaniun. If 
these two are mixed, the ratio is approximately 
four (or one-quarter) . 

Other Tests 

It should be clear from the foregoing that the 
test methods used in the B-l acoustical fatigue 
development program place emphasis on the produc- 
tion of damage. There is no more definitive evi- 
dence of convergence toward optimum design than 
failures for each component and fastener 
application. 

In contrast to this dependence on step-by-step 
destruction, two proof tests were performed with- 
out recourse to stepped loads. One of the sub- 
jects was a large two-box specimen that represent- 
ed the two .limiting load-to-strength conditions in 
a' graphite/epoxy- composite version of the horizon- 
tal stabilizer. The load was taken inmediately to 
the design level (including required margin) for 
the design life. No detectable damage occurred. 
This configuration, selected fcr A/C-S and beyond, 
would have provided a substantial reduction of 
weight from that of its metallic predecessor at 
similar cost. The other proof test was of a 
diffusion-bonded/superplastic- formed titanium 
access panel in the form of a corrugated sandwich. 
It demonstrated that the local details character- 
istic of this fabrication process do not degrade 
the acoustical fatigue advantages for which sand- 
wich panels and annealed 6A1-4V titanium alloy are 
known. 

Mechanical Fasteners 

Currently, only ve-y limited information is 
available regarding design requirements for 
mechanical fasteners for acoustical fatigue 
prevention. 


Tested Configurations 

The test programs for the fuel bay and stabi- 
lizer structural boxes provided definitive require- 
ments for fasteners appropriate for those con- 
figurations. Equation 2 is an example of this 
information for the stabilizer configuration. For 
the fuel bay, seven application situations exist 
for each of the individual frame assemblies; 
several other applications occur in the bay 
assembly, along longerons and intercostals , at 
longeron/frame intersections, and for skin 
splices. 

Figure 27 provides an example of a design 
requirement for A/C-4 established through the 
extension by analysis of information from the fuel 
bay tests. It is for the web- to- inner- cap rivets, 
except for those that also pass through the stiff- 
ener ends. (See Figure 14.) The inboard cap is 
restrained in translation and rotation by the 
slosh straps; unsymmetrical skin motion bows the 
web. This induces tensile load in the web-to-cap 
fasteners that is experienced principally by the 
fasteners at the ends of those stiffeners that tie 
to the inner cap inmediately adjacent to the inter- 
sections of the cap and straps. Thus, the spacing 
of web stiffeners is significant because it con- 
trols the allocation of reaction to roll of the 
frame between the stiffeners and the web. Note 
the similarity of form of this equation to that 
of Equation 2. even though the situations described 
are substantially different. This constant is 
larger principally because the fasteners at the 
stiffener ends carry most of the load. 
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The equation specifies the requirement in terms 
of a limiting allowable thickness of the upstand- 
ing leg of the outboard cap; all substructure com- 
ponents and fasteners are controlled from this 
quantity; it incorporates the effects of applied 
load, frame spacing, curvature, surface tempera- 
ture, and material properties. Two fastener types, 
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both having protruding heads, are permitted. The 
XF-6 is a 3/16-inch-diameter 2219-T81 aluminum 
rivet, and the AMP-5 is a 0.164- inch-diameter 
Cherry bimetal titanium rivet. The allowable 
range of upstanding leg thickness for each type of 
fastener depends on the rivet spacing, stated in 
the equation in inches and in the diagram as multi- 
ples of the fastener diameter. 

Other Configurations 

Somewhat surprisingly, almost no information 
about fasteners existed to support the design by 
analysis for the many concepts for which tests 
were not conducted. This necessitated formulation 
of analytical criteria, for a variety of fastener 
types and materials, applicable to the fastening 
of the skin and substructure of each of several 
conf igurations . Functional relationships similar 
to that of Figure 27 resulted that probably are 
valid, but the lack of data to provide constants 
of proportionality dictated a conservative bias 
in their use. 

Some information on skin fasteners for extend- 
ed skin and rib structure (Figure 28) was included 
in Reference 3, an indispensable source of data 
for acoustical fatigue analyses. However, the 
stated functional relationship, derived in the 
reference, is not that plotted. The information 
is useful only if the example superimposed on the 
nomograph reflects measurement results; if it 
does, it supplies a valid constant for the pro- 
portionality. That uncertainty aside, the in- 
formation of the figure is inherently limited in 
utility: 
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1. It is applicable only to skin and rib or 
skin and frame structure. 

2. It has no substructure information except 
that which can be produced by analytical extension. 


3. It apparently applies to in-phase motions 
of adjacent skin panels. 


Figure 29 presents a derivation, parallel to 
that idiich produced the proportionality in Fig- 
ure 28, for a skin- fastener relationship for 
out-of- phase - adjacent skin panels, a motion prob- 
ably substantially more significant than in-phase 
motion for acoustical fatigue damage. 
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Figure 30 suggests in an idealized form a 
laboratory procedure that would be useful in 
filling this void in essential design information 
by confirming the validity of, or supporting re- 
vision to, the functional relationships exempli- 
fied by that of Figure 29 and by providing the 
constants of proportionality. This discussion 
is intended to emphasize the principles involved 
and the« variables that require consideration. 

The procedure is cast in terms of skin and rib 
extended structure to provide a concrete exanple. 
Its key objective is to produce failures of fast- 
ener and skin (or fastener and rib, if desiredl 
in sufficiently close load/time coincidence 
that, with small adjustments through use of a 
fatigue curve and use of stress versus pressure 
data from the test, they will reveal the re- 
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quired conditions of compatibility between skin 
thickness and fastener strength. Thus, for a 
particular structural concept, the requirement 
for the fastener would depend only on its own 
properties, the geometry of the application, and 
the thickness of the skin, because that thickness 
is presumed to carry all necessary information about 
support spacing, load, curvature, and temperature. 

Significant variables for the investigation 
could include support spacing, angle versus tee 
supports (and, thus, one versus two rows of 
fasteners) , compatible versus rigid supports, skins 
chemically milled between supports or not, shear 
head versus tension head fasteners, and blind 
versus driven installations. Cost could be mini- 
mized through the use of specimens designed to be 
remanufacturable and through simultaneous loading 
of several specimens in a random-load facility. 

As suggested in Figure 30, the test might begin 
with rivets installed that were thought to be 
understrength (here symbolized by XD-4). Load 
would be applied in a series of increasing levels 
for perhaps 3 hours each until the rivets begin to 
fail. They then would be replaced by rivets of 
the next larger size (or greater strength). Load- 
ing in steps would proceed until the skin or the 
replacement fastener fails. That then would be 
replaced by the next size of skin or fastener. In 
this way, a combination of failures of rivets of 
two sizes (or strengths) and a skin, or vice 
versa, would provide information about the func- 
tional dependence and for the constant of 
proportionality. 

The following list suggests, beyond the fore- 
going particularized example, those more general 
objectives that could influence and, unfortu- 
nately, Sharply expand the test plan: 

1. Determine the character of the rivet fail- 
ure (e.g., shaft tension, head shear, head pry- 
ing); a difference in functional dependence of 
fastener strength on diameter, between the square 
and the cube, is involved. 

2. Determine the contributions of symmetrical 
and unsymnetrical motion. 

3. Learn to what extent the K values are 
"universal” constants, 

4. Find the equivalent strengths of important 
fasteners: aluninum rivets (shear and regular 
heads), bimetals, Huck bolts, Hi-Loks, and the 
blind fasteners: Jo-Bolts and Cherry locks and 
bulbs. 

5. Determine the appropriate "1" values for 
tees and angles. 

6. Evaluate the effects of skin splices and 
end bays. 


7. Establish the utility of "progressive 
remanufacture" of specimens as a cost-control 
technique. 

8. Repeat for or include substructure 
applications. 

9. Repeat, in part, for other than skin 
and rib structure. 

That there is a speculative tone to ranch of 
this discussion regarding fasteners reflects 
the state of knowledge in this aspect of design 
for acoustical fatigue prevention. 

Concluding Observations 

The primary objective of this paper is to 
exhibit some of the major characteristics of the 
design-support program for structural acoustics 
for the B-l. Emphasis has been placed on the 
approaches and techniques and on the points of view 
that led to them, with the anticipation that they 
may be found suggestive of methods applicable to 
future programs. To the extent that these tech- 
niques deviate from conventional practice, dis- 
cussion may ensue that could prove useful. 

As the contents of the paper suggest, there has 
existed a narrow focus in use of the results of 
analyses and tests toward airworthiness evalua- 
tions, the update of design requirements, and the 
support of weight- and cost- reduction efforts. 
Substantial quantities of information from analy- 
ses and measurements were generated in the pro- 
gram that have yet to be exploited; the engine- 
noise data for nonafterbumer operations and for 
a single engine are examples. Had the production 
program continued, it would have provided the 
opportunity and funding for these investigations 
of less inmediate applicability. 
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EXPERIMENTAL INVESTIGATION OF DYNAMIC CHARACTERISTICS 
OF TURBINE GENERATORS AND LOW-TUNED FOUNDATIONS 


S. P. Ying and M. E. Forman R. R. Drumm 

G ilbcrt /Commonwealth Penmylvania Power and Light Company 

Jackson, Michigan 49201 Allentown, Pennsylvania 18101 


Two 820 MW turbine-generator (T-G) units with low-tuned concrete foundations were 
investigated experimentally during both startup and normal operating conditions. The 
dynamic response of the turbine generators and their foundation systems was characterized 
by resonance curves obtained during the startup transient and by vibration mode shapes 
at the normal operating condition. No severe resonance at the fundamental frequencies of 
the system was observed. The deflection or displacement of the foundation pedestal does 
not lead to unacceptable T-G vibration. The experimental results are compared and dis- 
cussed with some previous theoretical studies. 


INTRODUCTION 

The definitions of the so-called “high-tuned” and “low- 
tuned” foundations for turbine-generators (T-G) are based on 
the fundamental frequency of the natural resonance of the 
system. This resonance is either higher or lower than the 
normal operating speed of the turbine generator. In partic- 
ular, these resonant modes themselves refer to the elastic 
vibration of the turbine-generator and its foundation pedestal. 
They do not include the mode t'lat relates to the vibration of 
the turbine-generator and its supporting structure as a body 
resulting from the softness of the foundation subsoil. 

The traditional high-tuned foundations are constructed of 
large solid concrete blocks. The low- tuned foundations 
possess either steel [1,2] or concrete pedestal columns [3]. 
This paper discusses measured dynamic responses for two 
low-tuned concrete foundations designed and constructed for 
twin 820 MW turbine-generators. 

Theoretically, the advantage of a low-tuned foundation 
lies in the reduction in bearing loads and forces which arc 
transmitted to the sub-founda.ion and the supporting soil. A 
practical advantage is the additional free available space 
under and around the low-tuned foundation pedest The 
larger free available volume under the turbine-generator 
permits a more generous layout of piping and auxiliary 
equipment. Moreover, the cost of low-tuned foundations is 
less than that of high-tuned foundations [4]. 

Because of these technical, practical and economical 
advantages, low-tuned foundations are quite popular in 
Europe and have recently been introduced in several U.S. 
power plants. The T-G foundations in this investigation are 
the first two built in this country. 


The foundations are part of Pennsylvania Power and Light 
Company’s Martins Creek Station, Units 3 and 4. The two 
820 MW oil-fired units are designed to operate at 1,000°F and 
2,400 psi steam conditions. Unit t Has been in commercial 
operation since October, 1975 with Unit 4 operating com- 
mercially since March, 1977. The plant is local i along the 
Delaware River approximately 75 miles north of Philadelphia. 

STEADY STATE OPERATION 

The concrete foundation pedestals shown in Figure 1 were 
designed and built with large amounts of reinforcing steel. A 
program for the experimental investigation of the dynamic 
characteristics of the two 820 MW units was divided into two 
phases. The first phase consisted of dynamic measurement of 
the entire turbine generator and foundation pedestal during 
steady state operation for both units. Acceleration data 
provided the necessary information for vibration charac- 
terization. 

The accelerations in orthogonal directions at each loca- 
tion were measured by accelerometers and recorded on 
individual channels of an FM tape recorder through charge 
amplifiers. This data acquisition system was calibrated in 
the field using a General Radio Model 1557-A vibration 
acceleration calibrator. The simultaneous use of multiple 
accelero-mctcrs permitted inspection of several adjacent or 
distant points with equal case. The data of a selected point, 
was recorded twice and each record contained at least two 
new data points. The repeated data point was used as a phase 
reference in rime domain. Approximately 500 data points 
were collected during the two unit survey. Acceleration 
frequency spectra were obtained from the data stored on 
magnetic tape by utilizing a Fast Fourier Transform (FFT) 
analyzer (Nicolet Scientific Corporation, Model 660). 
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f igure I, Low-tuned T-G foundationlpi dcstal (dimensions in meters J. 

Figure 2. Typical foundation pedestal acceleration during steady state, 800 MW {column 
3 from high pressure turbine, height 7.2 in. direction Xgl 
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Figure 3A. Mode shapes for the steady state operation of a low tuned foundation pedestal 
at 60 Hz in the longitudinal (Xj) direction. 


Acceleration data of the foundation pedestals for both 
units were investigated in the frequency domain. Figure 2 
shows a typical spectrum of acceleration in decibels (dB) 
with a reference acceleration of 10-6 x gravitational acceler- 
ation (g). The dominant vibration frequency is 60 Hz, which 
was detected at most measurement positions. Oscillation 
phase relations in the time domain were studied by using two 


band-pass filters and a dual channel oscilloscope to determine 
the oscill* <on modes at the dominant frequency of 60 Hz, 
corresponding to the T-G speed of 3,600 rpm. The oscillation 
displacements were calculated by using the double integra- 
tion in Nicolet 660 from the recorded accelerations on e.- h 
pedestal column and the platform in each direction. Then, 
the displacements were plott d as a function of distance. 


Figure 3B. Mode shapes for the steady state operation of a low tuned foundation pedestal 

at 60 Hz in the transverse (Xg) and vertical (Xj) directions ( Xj magnitude only). 
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These standing waveforms represent a mode of the forced 
oscillation resulting from the ce'tnfugal force generated by 
the rotor of the turbine-generator. 


STARTUP TRANSIENT 


Tlieoretically, the oscillation o the foundation |>edcstal 
is symmetric with respect to the tut. -tie axis. Figures 3A and 
3E depict the vibration mode in orthogonal directions. Xj. 
Xs. and Xy (see Figure 1), at one side ol the T-G axis of a 
unit, rhe standing waves on the other side and on the othet 
unit are similar; nevertheless, the magnitudes may varv as 
much as 6 dB. 


The second phase of this investigation included the con- 
tinuous monitoring of the turbine-generator and foundation 
pedestal at selec-d locations during a cold startup of a 
turbine -generator from aero to 3.600 rpm. These data were 
c >cd lor the transient vibration investigation of the low-tuned 
foundation. 


Most of the spectra obtained from the T-s bearings also 
show a peak at 60 Ha. A primary cause of the u0 Ha forced 
vibration is the imbalance ot the rotating system. Figure 4 ts 
a typical spectrum of such a vibrat: The data wc.c com- 

pared wilh the acceptable criterion loi general machinery - 
commercial ( 5 ] . Most ol the measured dynamic displacements 
of the TG bearings meet the maximum criterion. 


Tile accelerations on a turbine-generator and its foun 
dation pedcstel during cold startup were measured at seven 
locations, which were selected based oil the steady state data, 
and continuously recorded on individual channels uf FM tape 
recorders. Frequency spectra ot dynamic accelerations were 
obtained from the data stored on magnetic tapes by utilizing 
the narrow hand FIT analvzcr. 


On the other hand, a mis aligned coupling, which might 
result from the foundation, can generate l.urmonics at 1 20 Hz 
(see Figure 5). However, most 120 Hz components in the 
frequency spectra were found to be negligible, ar.d ah vihi 
tions possibly caused by the misalignment satisfied the 
maximum acceptable limit. Fi-.ally. observed !n v 
components at 30 Hz. a subharmonic frequency oi 
running speed, resulted from self-excited s-ibraimns Ucsacw 
of boundary conditions at the bea ri ng or its reL J steel 
elements. 


Figures 6 and 7 arc acceleration spectra mi sured at a 
pedestal column and a T-G bearing, respectively, during 
startup. As sl.isbtii in the spectra, the first mode oscillation 
frequency of each speefum is the T-G speed below 3.600 rpm. 
In contrast, the First inode oscillation in those spectra ob- 
tained during steady state operation (such as those shown in 
Figures 2, 4. and 3) is 60 Hz. its harmonics or its sub- harmonics. 


Dyna nic displacements of the first mode oscillations at 
the monitoring positions were calculated from (heir acceler- 
ations. The Fust inode displacements ire plotted vs. frequency 
'TG speed) for each posttiou (except (he front bearing of the 
high pressure (HP) turbine) during "a'tup in Figures 8 and 9. 


Figure 4 Typical hearing .tcce/er.trioit during .•ready stare. Jf 00 3/ It’ fgt neo-foi rear hear- 
ing. direction \'y. 
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F/jpur-.f. Typicjl bearing jceeUrtotion .iurittg strj.lv st.nr. $00 MU' (generator front branny, .iirrction .Vj 
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Typical acceleration spectrum during sunup. In.’t' 
9 ft rr. high, .lire. Iron X,. first rr.sit frc.)u,nc\ .’7 II 
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These graphs represent remnant cui.es at live locations. 
The three curves in Figure 8 are the resonant turves detected 
(t the rear bearing ol the !oa> pressure (LP) turbine, and the 
front and rear bearings of the generator. 


The data m Figure 9 were obtained from the filth and 
sixth columns counted from the HP turbine end. approxi- 
mately 9.(i m (Jl. S It' above the nut. The LP turbine rear 
bearing and the generator from bearing are at the top of the 
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Figure 7 Typical acceleration spectrum during startup, 3105 rym ilow pressure turbine 
rear bearing, direction Xg. peak frequency 51.75 H:). 


fifth column, and the generator rear bearing is at the top of 
the sixth column. The peak near 40 Hz detected on the fifth 
column also appears on the LP turbine rear bearing, and the 
peak at 27 Hz obtained on the sixth column is also shown on 
the generator rear bearing. 

The data obtained from the front bearing of the HP 
turbine in a vertical direction do not explicitly show resonant 
behavior The vibration accelerations increased gradually as 
the T-G speed increar d. 

Because the centrifugal force exerted on the turbine 
generator and : .s foundation during startup is smaller than 
that at a normal operating condition, the maximum accel- 
erations at the resonances are smaller than or comparable to 
those at operating condition. No severe resonance at the 
natural frequencies of the T-G and foundation was observed. 
Table I compares the accelerations in dB at the resonances 
during startup and at a normal operating condition. 

During startup of the unit an accelerometer was mon- 
itoring the vibration at the fifth column, same height (9.6m). 
of the T-G foundation-pedestal of rhe other unit. It was 
found that there was no noticeable change in violation 
amplitude at the monitoring location during startup. 

A more detailed study was conducted for a r-onitoring 
position on the sixth column (from the HP turbine end) of 
the startup unit. In addition to the first mode data, the 
second and third mode data were also analyzed and plotted 
in Figure 10. As shown in the figure, the • .tance fre- 
quencies of the first three modes are 27 Hz, 5i - 1 . and 70 Hz. 
There is no resonance at or near 60 Hz. All these modes were 


detected in the direction (X]) parallel to the T-G axis. The 
resonance frequency of another vibration mode detected on 
column No. 5 in the horizontal direction (X 2 ), transverse to 
the T-G axk> is about 40 Hz as shovm in Figure 9. 

DISCUSSION 

Natural frequencies of a large structure calculated from an 
Eigenvalue equation usually involve hequencies from very 
low to very high values, associated with many modes of 
vibration. However, the number of dominant resonance 
frequencies is quite limited. The data presented in the pre- 
ceding section indicate that there are only one or two domi- 
nant peaks in each resonant curve for these 820 MW units 
with the low-tuned concrete foundations. 

Wilson and Brebbia have estimated the dynamic response 
of 500 MW turbine-generators with low-tuned steel founda- 
tions [2]. Their calculated results of vertical displacements 
also show that there are only two dominant peaks in the 
frequency domain. Although the experiment described in this 
paper is not the same as their werk, the rewnance behaviors 
found by the two groups arc similar. 

In considering the entire T-G foundation-pedestal system, 
a reasonable model is proposed in Figure 11. The rotor 
pedestal, which includes a bearing pedestal and foundation 
column, and foundation are viewed as a series of masses 
connected by either viscous or clastic sections. The damping 
effects should be taken into account. Crook and Grantham 
used a similar model to predict the dynamic response of a 
375 MW T-G vrith a steel low -tuned foundation by the 
Prohl-Myklestad method (1 ]. They found that the resonance 
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Figure 9. Dynamic displacements or columns as functions of turbine speeds. 
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TABLE I 



Acceleration at Peak 

Steady State Test 


Startup Test 

Acceleration 

Turbine-Generator 

dB (re 10 6 g) 

dB (re lO^g) 

L. P. Turbine 

Rear Bearing 
Generate- 

51 Hz, 106.6 dB 

60 Hz, 107 dB 

Front Bearing 
Generator 

4*>,5 Hz. 97.3 dB 

60 Hz, 100 dB 

Rear Bearing 

49.5 Hz, 95.2 dB 

60 Hz, 97.2 dB 

Foundation-Pedestal Column, 9.6 M High 


5th 

40 Hz. 92 dB 

60 Hz. 92 dB 

6th 

27 Hz, 84.2 dB 

60 Hz, 84.5 dB 
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Figure l i. Theoretical mc.it ! of T-G anj ns low-tuned 
foundation 


at the journals is slightly less than 3.000 rpm. 

The experimental results obtaincJ in this investigation are 
similar. However. Crook and Gramham's calculated vibration 
amplitudes are much higher than those in this paper be- 
..a use they neglected the damping within the pedestals. The 
numerical values for the various parameters contained in the 
proposed model in Figure 1 1, except the damping coefficient 
r.i the pedestal, can lie found from Reference 1 . 

Wilson and Brebbia also calculated the mode shape for the 
case of free vibra.tou by the finite clement method. The 
calculated mode shapes arc consistent with the experimental 
results shown in Figure 3 of this paper. 

CONCLUSIONS 

At steady state, only the T-G bearing vibration at the har- 
monics of 120 Hr, resulting from a misaligned coupling, 
might relate to the deflection or displacement of the foun- 


dation pedestal. However, most 120 Hz components in the 
frequency spectra detected at the T-G bearingt are negligible. 
The deflection or displacement of the foundation pedestal 
docs not lead to unacceptable T-G vibration. 

No severe resonance at the natural frequencies of the T-G 
and its foundation system was observed during startup. The 
maximum accelerations at the resonances are smaller than or 
comparable to those at a normal operating condition. The 
resonance peaks detected at the different locations occur at 
27 Hz, 40 Hz, 51 Hz and 70 Hz. There is no peak at or near 
60 Hz (3.600 rpm). 

Furthermore, there was no noticeable change in vibration 
amplitudes ar a normally operating unit during the startup of 
the other unit. 

Based on the above observed results, it can be concluded 
that the low-tuned concrete foundations for these large 
turbine-generators provide fully satisfactory performance. 

The vibration mode shapes of the foundation pedestals 
during a steady state, and the resonance frequencies of the 
T-G and its supporting structure system obtained in this 
experiment, arc --onsistent with the previous theoretical 
studies. 
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ENVIRONMENTAL TESTING OF UHF BLADE ANTENNAS 
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Ser /ice problem! encountered on the F-4 aircraft showed that 
the environmental qualification requirements for the UHF blade 
antenna were not adequate. Although the individual environmental 
levels were stringent enough, only combined environmental testing 
could duplicate the conditions in which service problems occurred. 
Techniques for applying static side load simultaneously with vibra- 
tion and temperature were developed. The combined environment 
testing was instrumental in the rapid assessment of antenna modi- 
fications and resulted in a final configuration which proved 
satisfactory in service. A need to include combined environmental 
testing for qualification of blade antennas was established. 


INTRODUCTION 

The UHF communication blade antenna on the 
F-4 and other high performance aircraft has 
encountered periodic environmental problems 
over the years. The antenna is cantilevered 
from the aircraft moldline into the airstream. 









Difficulties on the F-4 were uncovered 
during vibration testing in 1964 and 1969. 
Again in 1974, an antenna incorporating elec- 
trical and product bility improvements failed 
during vibration qualification. Flight vibra- 
tion and airload (i.e. sideload) measurements 
were made. These verified the test levels. 
Supplemental high frequency vibration levels 
were also established. Bv early 1975, an 
antenna redesign with a tungsten powder damper 
passed all requirements and was incorporated 
on the F-4. 

However, in late 1979, complaints were 
received from the service that the antenna 
was developing cracks and some users con- 
sidered it unsatisfactory for flight usage. 

The structure of the antenna is shorn in 
Figure 1. It consists of two pieces of alum- 
inum separated by a dielectric insert at 
approximately mid-span. Structural integrity 
is provided by ’ steel pin with a tension 
preload, and the outer surface is . -ed 
with a polyurethane rain erosion ng. The 

cantilevered assembly is mounted ... ugh the 
base to the aircraft moldline. 

The service failures consistad of cracks in 
the rain erosion coating at the lower dielec- 
tric-aluminum interface (Figure 1, Item 11;. 
Disassembly of a returned antenna showed the 
dielectric insert to be disbonded and cracked 
completely across the interface. 
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Fig. 1 Sketch of UHF Antenna 


The service complaints resulted in a new 
investigation. Tests on production antennas 
showed them capable of meeting all individual 
qualification requirements. Since the individ- 
ual environments had been previously verified 
by flight data, it was concluded that the 
service problems had to be caused by combined 
environments. 

Examination of the antenna construction 
shows that a loss of preload, or sufficient 
moment to overcome the preload, would result 
in gapping, because of the Insert-aluminum 
interface has little or no tensile continuity. 
Since dielectric materials are relatively 
weak in bending, gapping would also account 
fcr the cracks in the dielectric insert and 
the erosion coating. 
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The primary environments causing moment at 
the Insert are vibration and static side load. 
The combination of the two could result In 
loads which exceed the preload along the face 
of the Insert. Since the preload and mechan- 
ical properties of the dielectric are also 
affected by temperature all three environments 
were considered important for combined testing. 

Test techniques were developed for simul- 
taneous application of side load with vibra- 
tion and temperature. These combined environ- 
mental tests were successful in reproducing 
the type of cracks reported in service. 

Subsequently, a modified antenna was 
developed which massed all individual and 
combined environmental tests. The modified 
antenna was approved for production and re'.ro- 
fit and has provided satisfactory service to 
date. 

The dynamic and combined environment por- 
tions of the antenna test program including 
the development of an improved side load 
method for combined environment testing, are 
presented below. 

REVIEW OF FLIGHT DATA AND QUALIFICATION VERIFI- 
CATION TESTING 

Flight test data were obtained on an instru- 
mented antenna in late 1974 to aid in evaluat- 
ing the vibration and airload qualification 
requirements. Data were obtained at repre- 
sentative F-4 flight conditions and maneuvers. 
Flight side loads were well below the specified 
test levels as were the vibration data. 

A composite of the reduced accelerometer 
data for the lateral direction at the antenna 
base is compared in Figure 2 to the qualifi- 
cation test spectrum. The flight data were 
reduced with 1/3 octave band techniques. 
Probability density analyses were performed 
on each 1/3 octave filtered time history to 
obtain a crest factor which was used to con- 
vert the rms acceleration level to a sinusoid 
at the center irequency of each 1/3 octave 
band. The flight data are well below the 
qualification test spectrum at all frequencies. 
The flignt data show peaks at 12 and 125 !!z 
and several above 500 Hz. These were easily 
predictable. The 12 Hz peak Is due to fuselage 
lateral bending. The 125 Hz peak is the first 
lateral bending mode of the ar.tenna, which 
has a somewhat lower frequency than the canti- 
levered bending mode due to the root flexi- 
bility provided by the aircraft structure. 

The higher frequency peaka are associated with 
higher modes of the antenna. 

Similar flight vibration data obtained at 
the antenna mid-spcn and tip were also well 
below the comparable laboratory responses. 

Due to the in-service problems, however, It 
was decided to verify that production an- 
tennas would still meet qualification test 


requirements. An untenne was taken from 
stores and subjected to all individual environ- 
ments including vibration at ambient and ele- 
vated (127°C) temperatures, shock (15g, 11 
millisecond, half-sine) and static side 
loading to 5.512 x 10* newtons per square 
meter (8 psi). The antenna successfully com- 
pleted all testing with no structural or 
electrical failures. It was evident therefore 
that the individual flight environments did 
not produce the flight problems. 



Fig. 2 Competitor of Flight Data and Qualification Tsrt 
Spactrun 

ANTENNA Via RAT ION CHARACTERISTICS 

Tests were performed on antennas to obtain 
response characteristics and mode shapes to aid 
in determining the vibration modes relevant to 
the flight failures. A typical cantilevered 
antenna frequency response plot of tip trans- 
missibility is shown in Figure 3 for the 
latere axis. 

The associated lateral vibration modes are 
presented in Figure The first mode is 
easily identified as -.antilevered bending, 
whereas the ’-emaining three are primarily 
associated with tip motion. From the response 
and mode shape data, it appeared that the pre- 
dominant mode involving high moments in the 
area of flight failures was the first lateral 
bending mode. 
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■ In order to further assess which modes were 

i important, vibration responses of a cracked 

antenna returned from service were compared 
, with those of an undamaged antenna, using side 

loads of 0, 44.5 and 89 newtons. The results 
at the qualification vibration levels are pre- 
f sented in Table 1. 


trie insert for the cracked antenna. This 
indicates gapping and a loss of stiffness 
across the Joint, resulting in a reduced 
frequency. 

These data tended to confirm a failure 
mechanise associated with a combination of the 
fundamental lateral bending mode and side 
loading. 

Table 1 Compa ri son of Latsml Resonant Frequencies of Now 
and Cracked Antennas 


TIP 

SIDE LOAD 

(N) 

INPUT 

VIBRATION 

dd 

RESONANT FREQUENCIES (Hi) | 

NEW ANTENNA 

FAILED ANTENNA 


10 

260 

236 


36 

060 

630 

0 

26 

830 

790 


26 

1040 

1010 


26 

1800-1960 

1700-1800 


10 

280 

116 

44.6 

36 

- 

630 

(10 LB) 

25 

800420 

730 


26 

1000-1060 

620 


26 

1700-1960 

1690 


10 

240 

80-90 

80 

36 

- 

630 

(20 LB) 

26 

810 

710 


26 

1020 

870 


26 

1760 

1580 
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BASE SPANWISE STATION 



The date showed no significant difference 
without the side load. However, the addition 
of side load had a dramatic effect on the first 
bending mode frequency of the cracked antenna, 
dropping it from 235 to 80 Hr, while the new 
antenna frequencies dropped only slightly. 

The higher mode resonances were changed only 
slightly by elde load. 

An examination of the first bending mode 
shapes showed evidence of gapping in the 
cracked antenna under combined load as shown 
in Figure 5. For the no load case, the mode 
shapes and frequencies are quite similar while 
for the sideload case, there is a distinct 
break in curvature in the region of the dielec- 


Fiq. 5 First Mode Comnsrison of Nsw and Cracked Airtwmat 
with and without Side Load 


COMBINED ENVIRONMENT TESTING 

In order to duplicate the flight failures, 
a new antenna was subjected to combined vibra- 
tion, side load, and tenperature environments. 
A soft spring system was used to apply side 
loading. 

A schematic of the test setup is shown in 
Figure 6. The uniform pressure was correlated 
with tip loading by equating moments at the 
lower center of the dielectric lnserc. 
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Note: Sideload applied at 50% chord and approximately 
0.41 cm (0.16 in.) inboard from tip. 


Fig. 6 Schematic of Combined Environment Test Setup 
Using a Soft Spring for Tip Side Loading 


The test parameters were based on the 
maximum ranges anticipated during flight for 
the F-4 . Testing was performed on several 
antennas, as summarized in Table 2. It con- 
sisted of resonant dwells for the critical 
first lateral bending mode at vibration levels 
of +3, 5 and 10 g's; at room temperature, 

93.3°C (200°F) , and -40°C (-40°F); with tip 
side loads of 0, 44.5, 89, and 133 newtons 
(0, 10, 20 and 30 pounds) . 

Table 2 Summary of Combined Vibration/T emparatura/Sido 
Load Supplemental Antenna Tatting 
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In general, there was a significant drop In 
frequency with Increasing side load at all 
temperatures. (Apparent Inconsistencies in 
this trend were attributed to coupling between 
the soft spring side load setup and the an- 
tenna.) These frequency reductions are Indi- 
cative of gapping at the dielectric Insert. 

During the final test point for the first 
antenna, there was a significant drop In 
resonant frequency during the dwell, indicat- 
ing probable failure. Post test examination 
revealed cracks In the polyurethane coating and 
the dielectric insert similar to those found 
in the returned service antenna. The second 
and third antennas both failed during the last 
portion of the high temperature testing, and in 
both cases the steel preload pin broke. 

Tests were also conducted to evaluate the 
effects of increasing the pin preload. The 
results showed that this helped but the dielec- 
tric insert still cracked during the tests. 

In parallel with the analysis and testing 
being conducted by MCAIR, the vendor developed 
a modified antenna which provided significant 
structural improvement. A fiberglass strap 
was bonded to the dielectric insert and the an- 
tenna base and tip as shown in Figure 7. Two 
of these antennas were provided to MCAIR for 
test. 



FIBERGLASS 

STRAP 


IH I IWI 

Fig. 7 Sketch of Modified Antenna Showing Fiberglass Strap 


The first antenna was subjected to the 
vibration qualification test spectrum and to 
supplemental high frequency testing in the 
lateral axis. The fundamental lateral bending 
mode was found to be approximately 240 Hz. 
Shock and static side loa! tests were also 
performed. The antenna passed these tests 
without any structural or electrical degra- 
dation. 
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The second antenna waa aubjactad to com- 
bined vibration/ temperature/aldeload teating 
aimilar to that applied to the unmodified 
Serial No. 0897 antenna in Table 2. It 
exhibited atable fundamental mode frequency 
characteristics with aide load at all tempera- 
turea , although mechanical coupling with the 
side load aprlng ayatem continued to be a 
problem. 

Additional combined teating waa performed 
at high temperature (93.3*C) with tip side 
loads of 133 newtons in an attempt to test 
the antenna to failure. The vlbratloi exciter 
ayatem being used was nearing its maximum 
capability of 35 g's (at 244 He) when the test 
setup waa damaged and testing was discontin- 
ued. There was still no structural or elec- 
trical damage to the antenna. 


to fit various antanna shapes and to control 
the air gap. which was approximately 0.025 
to 0.050 centimeters (0.010 to 0.020 in.) 
for the UHF antenna configurations tested. 



The structural integrity of the modified 
antenna was verified by the tests performed 
including combined environmental levels 
which equal or exceed those anticipated on 
the aircraft and which caused catastrophic 
failure of the unmodified version of the 
antenna. This modified version has been in 
use on the F-4 without field problems for more 
than tvo years. 

DEVELOPMENT OF AN IMPROVED SIDE LOADING SETUP 
FOR COMBINED ENVIRONMENT TESTING 

As mentioned above, numerous problems were 
encountered with the soft spring loading system 
in the test setup. Mechanical coupling and 
damping response effects made it difficult to 
Interpret test response changes. Also, de- 
sired bending moments and shears could only 
be approximately matched, and test efficiency 
was poor due to the need for frequent loading 
system adjustments, especially when the access 
was required to the thermal enclosure which 
then had to be restabilized. 

Several different types of springs and 
attachments were tried. Although improvements 
were made, the coupling effects and damping 
problems persisted. Therefore, another means 
was needed for applying side load. 

The system decided upon was compressed air. 

A sealed air cavity was constructed beneath 
the antenna, with a small air gap around the 
edges of the blade. Cjmpressed air was in- 
troduced into the cavity at a rate faster than 
it could leak out around the gap between the 
blade and the fixture. By regulating the flow 
rate, it was possible to hold a uniform over- 
pressure on one side of the antenna. Static 
pressures up to 5.512 x 10^ newtons per square 
meter (8 psl) were easily obtained and even 
higher pressures could be obtained if desired. 
It was also found that high or low tempera- 
tures can be obtained rather quickly by pre- 
heating or precooling the inlet air supply. 

This system is shown in Figure 8. The 
template on the top of the cavity can be cut 


Fig. 8 Combined Environment Test Setup Using Compressed 
Air for Side Loading 


Frequency response data were obtained on an 
antenne to evaluate the teat setup. The side 
loading was conducted in the same test setup 
at room temperature, both with the soft spring 
tip side load mechanism and the compressed air 
system. 

Tie no load baseline, and representative 
equivalent tip, and uniform pressure side 
loading responses are presented in Figure 9. 

As shown, the fundamental mode responses using 
the air loading are cleaner and have less 
damping (lie. higher amplitude response). 

The higher mode responses are similar with 
either loading system. 


TtFSPNmOSIJl LOADS 




1.080 tO.OSO 
FREQUENCY Hi 


Fig. 9 Frequency Response Comperiton with Spring end Air 
Side Load 


In order to fully check ov.l the setup, a 
new unmodified antenna was tested se com- 
bined environment compressed air fi. arc. 
Results are summarized in Table 3. The antenna 
lasted somewhat longer than previously tested 
unmodified antennas; however, the pin bro^e 
and the dielectric insert was completely 
severed at the 2.756 x 10^ newtons per square 
meter (4.0 psl) aids load, 93.3*C (r00*F) 
test point while peaking for the lOg resonance 
dwell. 
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CONCLUSIONS 


Tsbls 3 Summary of Combined Tatting on Original Antenna 
with Comprattod Air 
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'Antenna tailed while peaking at 10 g, pm broke, dielectric severed at lower interface 


with aluminum and antenna separated into two pieces. 
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The following conclualons were made on the 
basis of the analysis and testing program for 
the F-4 UHF antenna: 

o Blade antennas can be located so that 
significant steady air loads result. 

These may ccmblne with vibration over a 
broad temperature range. Thus, separate 
tests may be Inadequate. 

o Combined vibration/ temperature/ sideload 
tests can be developed which duplicate 
actual service conditions sufficiently 
well to duplicate service failures and 
evaluate design modifications. The 
antenna so modified in this paper has 
performed satisfactorily on the F-4 
in service to date. 

o The best method developed for applying 
side loading with vibration and tempera- 
ture Is a compressed air chamber with a 
non contact seal . This Is recommended 
for future combined environment testing 
of blade antennas. This method Is also 
applicable for standard Individual static 
side load tests. 

o Test requirements for externally unted 
blade antennas should Include combined 
vibration, temperature and side lead 
testing, either as an option dependent 
on antenna design, or as a mandatory 
test. 


The trend of decreasing resonant frequency 
with Increasing side load, typical of this 
type of antenna. Is evident in the drta, 
although the percentage drop isn't quite as 
large as for antennas tested with the spring 
system. The maximum antenna stiffness reduc- 
tion, again fpic’l.ly, occurred at high tempera- 
ture. 


Iu. test setup checked out very well and 
'i 1 - _gh and low temperature conditions were 
r-_ ny obtained. The change of side load was 
much easier to accomplish at temperature and 
the test setup was found to he satisfactory 
in all respects. 

There are several ad vantages of compressed 
air for combined loading. The uniform pressure 
produces correct shear and moment distribu- 
tions. There is no mechanical coupling, cj 
that damping and response effects due to the 
loading system itself are minimized, making 
test responses easier to interpret. The side 
load ..an be changed rapidly and accurately 
with temperature effects induced in the air 
supply, which optimizes test time. Addition- 
ally, the loading system can be used to perform 
individual static side load Lests. 
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SHOCK ISOLATION FLAT FORM FOR SEASPARROW LAUNCHER 


Paul V. Roberts 
Raytheon Company 
Missile Systems Division 
Bedford, Mass. 


* 8 ' 0 -17302 


A unique low frequency shock isolation platform is used to mount 
the NATO SEASPARROW Launcher aboard U. S. shock-hardened ships. 
The platform attenuates shock response to less than structural design 
loads which were used to minimize weight on smaller NATO ships and 
for compatibility with airborne missiles. 

The platform also isolates shipboard vibration in the most critical 
range of hunch cell and missile resonances between 14 and 30 Hz, and 
furnish -• i lightweight convenient means of installing the launcher 
aboard * .p. 


INTRODUCTION 

The NATO SEASPARROW Launch -r, 
mounted on the platform, is shown in Figure 1, 
on the aft main deck of the DD963. The 
L jncher carries eight (8) SEASPARROW mis- 
siles and operates through 360” train and -10” 
to 90” elevation. 



Figure 1 - Platform and Launcher on DD 963 

The launcher and missile supports were 
designed to specified underwater shock static 
equivalent loads of 1 5 g vertical, 9 g trans- 
verse and 6 g longitudinal to the ship. These 
loads were taken as simultaneous combined 
limit loads for structural analysis purposes. 

The Sparrow missile is designed for flight loads 
and the missile attachments to the launcher are 


designed for aircraft carry loads which are 
compatible with the launcher combined limit 
loads. 

Use of low frequency 6 tc 10 Hz mounting 
was recommended at the beginning of the NATO 
SEASPARROW program as the be9t means of 
not exceeding structural design loads due to 
underwater shock. Analyses during program 
development using the DDAM method of Nav 
Ships 2 50-423-30 confirmed the desirability of 
a low frequency platform as indicated i" the 
shock spectrum curves of Figure 2. 

The capability of the launcher to withstand 
underwater shock on the recommended plat- 
form was demonstrated by successfully with- 
standing full Mil-S-901C underwater shock 



moix** - ' I H| 


Figure 2 - Shock Spectrum versus Deck 
Frequency 



barge teats on an engineering design model 
(EDM) platform. These tests were conducted 
at the Hi Test Lab Facility in Arvonia. VA. 

The plate shock isolation platform, cov- 
ered in this paper, was verified as satisfac- 
tory from a MU-S-901C underwater shock 
standpoint, by DDAM analyses showing that it 
was superior to the EDM platform. 

In addition to underwater shock verifica - 
tion, the paper covers platform design, devel- 
opment, ease of ship installation, shipboard 
vibration isolation, weight reduction and plat- 
form structural analyses. 

PLATFORM RATIONALE 

Principal engineering philosophies upon 
which the platform design was based are out- 
lined as follows: 

1) Provide low vertical and transverse/ 
rocking natural frequencies between 
about 4 and 8 Hr to attenuate under- 
water shock. 

a) The platform basically protects 
the launcher, and missile sup- 
port'’ to less than structural 
design loads, when subjected to 
Mil-S-901C shock. 

2) Isolate Mil-Std-167 shipboard vibra- 
tion at principal launch cell and mis- 
sile resonances above 14 Hz. 

a) The 4 Hz to 8 Hz natural frequency 
range isolates vibration above 14 
Hz and avoids most severe blade 
passage excitation, which is 
above 8 Hz on most ships. 

3) Incorporate damping to limit response 
at 4 Hz to 8 Hz platform frequencies 
due to Mil-Std-167 inputs. 

a) Although excitation at shaft speed 
occurs in the 4 Hz to 8 Hz range 
on many ships, acceleration in- 
puts are only *0. 05 g to ±0. 2 g at 
these frequencies. 

b) Hydraulic dampers with C/C c = 

0. 25, limit response in the 4 Hz 
to 8 Hz range to less than *l/2 g. 

4) The platform must be convenient to in- 
stall directly to ships decks without 
auxiliary support structure. 

a) The platform ring, bolts or welds 
to a mating ring which is welded 
to the ships duck. Secondary sup- 
ports are necessary to attach the 
dampers to the ship. 

5) Toe platform must be within the 
launcher swing envelope for minimum 
deck space and be light weight for use 
aboard smaller NATO ships. 


61 The launcher must meet all opera- 
tional performance criteria when 
mounted on the platform, such as 
response times to position commands, 
servo stability and pointing accuracy. 

PLATFORM DESIGN 

The principal element of the platform con- 
sists of a 2.64 meter diameter 606I-T6 alumi- 
num plate. The plate tapers from 7. 62 centi- 
meter thickness at the inner launcher attach- 
ment interface to 2. 54 centimeter at the outer 
ship attachment ring. The plate and other 
principal platform elements are shown schema- 
tically on Figure 3. 


PLATE GASKET 



Figure 3 - NSSMS Launcher Platform 

The plate is tapered for more uniform 
stress distribution, weight reduction wherein 
the thinner outer portion comprises most of the 
area, and for shedding sea and rain water off 
the upper surface. Stresses are still highest 
at the thicker launcher attachment interface, 
but the plate thickness was limited to 7. 62 cm 
to reduce machining waste and material cost. 

Parametric studies using the ASSAS finite 
structural element computer model described 
under structural analyses in Figure 14, showed 
that the inner to outer taper was efficient from 
the standpoint of maintaining high rocking mom- 
ent stiffness for servo stability and pointing 
accuracy along with low vertical stiffness for 
shock isolation. The inner to outer taper was 
found superior to flat plates and reverse tapers 
from these stiffness standpoints, within allow- 
able plate stresses. 

In addition to the plate, the other major 
metallic element in the platform is a 1. 27 cm 
thick steel angle ring used to support the outer 
plate circumference. The mounting ring is 
welded and bolted to a mating ring welded to the 
ship* s deck. The deck attachment is described 
in more detail under shipboard installation. 
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A 3. 18 cm thick neoprene seal is used be- 
tween the lower surface of the outer plate and 
the mounting ring. In addition to water sealing, 
the neoprene seal furnishes additional down 
flexibility and provides a softer pinned as com- 
pared to clamped structural joint at the outer 
plate support. In the up direction, flexibility 
as well as sealing is provided by 48 neoprene 
bushings located below the mounting ring as 
shown on Figure 3. 

In addition to the lower neoprene bushings, 
48 silicone bushings are inserted in the plate 
from the upper surface to provide overall seal- 
ing as well as flexibility around the outer plate 
attachment bolts. The silicone material was 
used in the upper bushings to resist heat due 
to occasional rocket motor blast inpingment on 
the upper surface of the plate, during missile 
firings at high elevation angles. 

From a dynamic standpoint, the neoprene 
seal and bushings primarily furnish flexibility 
in the plane of the plate for isolation of trans- 
verse and longitudinal underwater shock. The 
elastomeric flexibility allows about 2. 1 cm of 
relatively soft transverse plate motion and pro- 
tects elements close to the launcher base such 
as the slip ring assembly, which would be more 
directly subjected to deck plane shock. 

In the principal vertical shock direction, 
the plate is much more flexible than the elas - 
tomeric elements and undergoes about 88 per- 
cent of the overall vertical deflection. 

Forty-eight shoulder bolts are used to 
maintain compression on the outer neopren? 
seal to prevent leakage during sea state conii- 
tions. Short period seal separation would 
occur at response motion peaks resulting from 
underwater shock. The launcher base ring is 
attached to the platform on a 0. 81 meter bolt 
circle by means of 2. 54 cm diameter bolts. A 
0. 32 cm thick 70 durometer neoprene gasket is 
used between the launcher base and platform 
for sealing and to relieve stresses in the 
launcher base due to radial moment loading re- 
sulting from plate deflections. 

The four hydraulic dampers shown in Fig- 
ure 3 primarily furnish vibration damping in 
the vertical direction since Mil-Std-167 input 
accelerations are much greater at the higher 
vertical frequency. Little damping is needed 
at the lower 4. 1 Hz transverse frequency since 
the Mil-Std-16? input is only ±0. 052 g. com- 
pared to * 0. 18 g at the vertical frequency. 
Damping measured experimentally during dy- 
namic tests averaged about C/C c = 0. 25 at the 
7, 6 Hz vertical frequency, compared to C/C c = 
0. 036 at the first 4, 1 Hz transverse frequency. 
The dampers incorporate force limiting at high 
underwater shock imput velocities such that the 
damper attachments were designed for a 2000 
lb. (8896N) peak force. 

DEVELOPMENT ANALYSES AND TESTS 

The following principal analyses were con- 
ducted during platform development. Results 


of these analyses are discussed in separate 
sections as appropriate. 

ASSAS computer parametric stress and 

deflection analyses 

DDAM shock response analyses described 

in later paragraphs 

NASTRAN vibration response analyses 

Elastromeric element load and stiffness 

analyses 

Structural analyses of principal elements 

Static, dynamic and servo tests were con- 
ducted during development as follows. 

Static Tests 

1) Vertical down to 63 percent load 

2) Vertical up to 49 percent load 

3) Transverae/rocking to 37 percent 
load 

Dynamic Twang Teats 

1) Vertical 

2) Rocking 

Servo Tests 

1) Response times to step commands 
and pointing accuracy 

2) Stability gain margin 

3) Response to sinusoidal frequency 
sweeps 

Results of the static stress measurements 
were used to extrapolate stresses up to full 
underwater shock design loads. Experimental 
stresses varied locally from stresses predicted 
analytically primarily due to bolt hole cut outs. 
Deflections measured in the static tests were 
used to verify stiffness analyses and for final 
DDAM and vibration analyses. 

Dynamic twang tests were conducted in the 
vertical and transverse directions to verify 
natural frequency and damping ch racteristics. 
Step release loads were suddenly applied at pre- 
dicted levels by means of shear bolts. 

The servo performance tests were con- 
ducted wich 8, 4 and 0 missile loads on the plat- 
form, and with the launcher base rigidly 
attached to a concrete deck. Results of these 
t *sts showed that response times to required 
elevation and train angles were almost identical 
on and off the platform, and that transient 
motions were well within pointing accuracy re- 
quirements. In addition, both rate and position 
loop gain margins were fully satisfactory on the 
platfo/m. 

The elevation rocking frequency resulting 
from the sinusoidal frequency sweeps, reduced 
from 6. 1 Hz hard mounted to 4, 1 Hz on the plat- 
form with an eight missile load. This reduction 
was anticipated due to inherent platform flexi- 
bility provided for shock and vibration isolation. 
Train frequencies remained essentially the 
same both on the platform and hard mounted. 

Basic dynamic information resulting from 
tests and analyses are summarized below for 
the most severe eight missile load condition. 
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Vertical natural frequency 7.2 Hr (shock) 

Vertical natural frequency 7.6 Hz (vibra- 
tion) 

Vertical damping C/C c = 0. ?5 
Vertical operational travel *4. 3 cm 
Maximum down vertical travel -10. 2 cm 

Tranaverse/rocking natural frequency 
-4. 'Hz 

Transverse/rocking damping C/C c =0. 036 

Transverse operational travel at plate 
±1. 6 cm 

Maximum transverse travel at plate 
±2. 2 cm 

Predicted operational travel at eg ±3. 8 cm 


direction only, at both docks by means of leaf 
springs or other suitable flexures. In addition 
to being heavy, the trunnion takes up consider- 
able below deck space and makes access to the 
launcher base difficult. 



El-vation and train frequencies with vari- 
ous missile loads found during the servo test 
frequency sweeps are listed as follows: 

No. Missiles Elevation Train 

0 4. 9 Hz 6. 6 Hz 

4 4. 3 Hz S. 4 Hz 

8 4. 1 Hz 6. 2 Hz 

SHIPBOARD INSTALLATION 

In addition to reducing equipment weight 
due to permitting lower shock design loads, the 
platform comprises a convenient lightweight 
method of mounting the launcher aboard ship. 

Figure 4 illustrates basic advantages in 
installing the large 2. 64 meter diameter plat- 
form on the deck compared to attaching the 
0. 89 meter diameter launcher base directly to 
the deck. The rocking stiffness R, in elevation 
directions normal to the deck, must be main- 
tained for pointing accuracy and servo stability. 
For a given rocking stiffness R the vertical 
stiffness K. given by K = 8R/D^, is inversely 
proportional to the diameter squared. With the 
smaller launcher diameter therefore, the deck 
stiffness K, for the same R, is 8. 8 times 
higher than with the platform. 

Due to the high vertical to moment stiff- 
ness without the platform, it is difficult to pro- 
vide sufficient rocking stiffness for satisfactory 
pointing accuracy and servo stability and at the 
same time be sufficiently flexible in the verti- 
cal direction to protect the launcher and mis- 
siles from underwater shock. 

Ship installations without the platform use 
additional deck structure for satisfactory oper- 
ational moment stiffness but end up being too 
stiff vertically for Mil-S-901C underwater 
shock protection. A possible approach, which 
is inherently much heavier than the 681 kg plat- 
form, is to attach the launcher to a trunnion ex- 
tending from the launcher deck to next lower 
deck. The trunnion provides high moment 
stiffness and is soft mounted in the vertical 
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Figure 4 - Deck Loads with 0. 89 m Launcher 
Base versus 2. £.4 m Dia Platform 

Figure 4 also illustrates basic advantages 
with the larger 2. 64 meter platform from the 
standpoint of much lower deck loads at the 
platform mounting ring. For a given over- 
turning moment M the peak load F per running 
unit of deck length is approximated by Fr 
4M/D-. The force F is therefore also in- 
versely proportional to the diameter squared 
such that with the smaller 0. 89 meter launcher 
stand the peak load is again 8. 8 times higher 
than with the 2. 64 meter platform. 

With the lower deck loads and amount of 
deck structure swept by the 8. 29 meter plat- 
form circumference, additional deck support 
structure has been found unnecessary. The 
ship attachment ring has therefore been welded 
directly to the deck and the matching platform 
ring has been tack welded and bolted to the deck 
ring. Secondary structure is needed for attach- 
ment of the four hydraulic shock absorbers. 

SHIPBOARD VIBRATION 

Although primarily designed for under- 
water shock ; «olation, •’he low 4 Hz to 8 Hz plat- 
form natural frequencies also isolate shipboard 
vibration in the most critical range of launcher 
and missile resonances between about 14 Hz 
and 30 Hz. 

Principal advantages of the platform from 
a shipboard vibration standpoint are illustrated 
on Figure 5. Launcher rigid body vertical re- 
sponse to the shipboard vibration inputs of Mil- 
Std-167 are shown on Figure 5 both on the plat- 
form and with rigid mounting. 

As indicated by the lower platform response 
curve, vibration isolation occurs above about 
l 1 Hz with the 7. 6 Hz platform vertical fre- 
quency and progressively improves at higher 
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frequencies. The upper curve for rigid mount- 
ing consists of required Mil-Std-167 input 
accelerations, since inputs would be applied 
directly to the launcher base with rigid mount- 
ing. 



Figure 5 - Base Acceleration with Rigid 
Mounting and on Platform-Mil -Std- 167 Input 

Reference to Figure 5 shows that peak 
Mil-Std-167 inputs occur at 25 Hr, 33 Hi and 
15 Hi. These peaks correppond to missile 
cantilever resonances between 20 Hi and 25 Hz. 
free missile resonances between 28 Hz and 33 
Hz and the launch cell vertical cantilever reso- 
nance at 14 Hz to 18 Hz. When mounted on the 
platform, inputs to the launcher are basically 
attenuated about 50 percent to over 80 percent 
from the low to high range of these resonances. 

Beefing up the launcher cell structure to 
raise the 14 Hz vertical frequency above 33 Hz, 
is questionable from a technical standpoint 
without prohibitive wei(;nt increase. Various 
methods of restraining missile body bending 
resonances in the launch cells also invariably 
lead to additional weight, complexity and inter- 
ference with misfile launch. 

Vibration response analyses were made 
using the NASTRAN computer program and the 
launcher/misscle dynamic models used for 
DDAM underwater shock response. These 
analyses were made on the plate platform dis- 
cussed herein, the EDM platform and with the 
launche • rigidly mounted. 

Results oi the NASTRAN vibration analyses 
are illustrated on Figure 6, which shows re- 
sponse of lump mass upper outboard missile 1 
to Mil-Std-167 inputs. It is seen that response 
on the plate and EDM platforms is much more 
f. .vo rable than rigid mounting and that response 
on the plate platform is most favorable. The 
effect of the higher damping used with the plate 
platform in reducing response at the 7. 6 Hz 
vertical resonance is also evident o:, Figure 6. 
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Figure 6 - Missile 1 Response Mil-Std-167 Input 

The ability of the launcher and missiles to 
withstand Mil-Std-167 vibration on the EDM 
platform was verified during tests at Ogden 
Technology Labs using a large hydraulic shaker. 
Since the NASTRAN analyses indicate more 
favorable vibration characteristics on the plate 
than on the EDM platform the launcher rnd 
missiles should satisfactorily withstand Mil- 
Std-167 conditions on the plate platform. 

Inasmuch as the plate platform is highly 
damped in the vertical direction with C/C £ a 
0. 25, calculated response at the 7. 6 Hz 
vertical frequency with the ±0. 18g Mil-Std-167 
input would be approximately ±0. 4g. Since no 
internal resonances occur in the region of 
7. 6 Hz, response of the overall launcher mis- 
sile assembly should not exceed ±l/2g at the 
platform vertical resonance. 

The transverse/rocking resonance on the 
platform of 4. 1 Hz with eight missiles, is not 
nearly as well damped as the vertical 7, 6 Hz 
resonance. Since the Mil-Std- 167 input at this 
lower frequency is only ±0. 052g however, re- 
sponse accelerations would still be on the 
order of ±l/2g with damping measured during 
dynamic tests. 

Vibration on most ships is well below Mil- 
Std-167 levels, especially during normal cruise 
conditions. On ships with shaft rotation speeds 
in the region of 4 Hz, low level transverse 
vibration could occur on upper launcher ex- 
tremeties up to about 0. 25 cm double amplitude. 
The peak ±l/2g Mil-Std-167 response level, 
corresponding to about 1. 5 cm double amplitude 
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would only be anticipated during severe maneu- 
vers such as high speed turns and backing down. 
Similar levels could also occur due to ship 
damage, especially to the shafting and blades. 

Hard mounted vibration tests conducted at 
the Raytheon Sudbury Facility, wherein Mil- 
Std-167 levels were applied at the launcher 
base, verified the launcher and missiles to with- 
stand Mil-Std-167 levels rigidly mounted. Two 
series of tests were conducted. The first up Co 
9 Hx verified the launcher and missiles for 
installation on AOR and aE type ships having 
blade passage frequencies below 7 Hz. The 
second series of tests used auxiliary missile 
snubbing devices for overall Mil-Std-167 verifi- 
cation. 


each of the four ends of the beam structure 
were oriented at 45 deg in order to reduce 
vertical compared to lateral and longitudinal 
stiffness and to partially decouple translational 
and rotational modes. 

A comparison of characteristics of the 
production plate and EDM platform, which 
principally effect underwater shock isolation, 
is tabulated below. The vertical and trans- 
verse/rocking frequencies are lower on the 
plate than on the EDM platform and the plate 
platform is more highly damped. All of these 
factors are more favorable from a shock isola- 
tion standpoint. 

pi$C COM 


UNDERWATER SHOCK VERIFICATION 


VERTtCAL FREQUENCY 7.2 Hi 


7.8 Hi 


Verification of the shock isolation platform 
to protect the launcher and missiles from Mil- 
S-901C shock, was based on DDAM analyses 
showing that the platform would furnish better 
shock protection than an engineering design 
model, EDM, platform on which the launcher 
passed full underwater shock tests. 

The EDM platform is shown in the photo- 
graph of Figure 7, as set up for underwater 
shock barge tests at the Hi Test Lab Facility 
in Arvonia, VA. The platform structure con- 
sisted of two longitudinal and two lateral, 

27. 6 cm depth by 3. 66 meter length, wide 
flange steel I-beams about 0.91 meters apart. 
The launcher was bolted to a 3. 18 cm thick 
steel plate welded over the I-beam center sec- 
tion. The four ends of the structure were boxed 
in by means of 30. 5 cm depth channels welded 
to the I-beams. 


VERTICAL C/C c 
ROCKING FREQUENCY 
VERTICAL K 

TRANSVERSE K 
ROCKING k 
WEIGHT 


0.25 0.1 


4.1 Hi 

13.8 (101 6 N/m 
(79 000 Ibv'l") 

7.2 (10l 6 N/m 

10.6 (10) We 

7037 kg (15,500 


7.4 Hi 

25 7(10) 4 N/m 
(147. 000 Ibs/in.) 

16.5 00) 6 N/m 
45.2 (10) Wm) 
8989 kg 09,800 lb) 


The 1500 lb. (681 kg) plate platform is 
much lighter than the 5800 lb. (2633 kg) EDM 
platform as indicated above. Stiffnesses in 
principal directions, also listed above, show 
the same trend as natural frequencies, but can- 
not be compared directly, especially ir. trans- 
verse/rocking modes, since they are effective 
at different locations. 


Twenty Nichols Engineering sandwich 
mounts were installed between the four ends of 
the platform and deck by means of 0. 91 meter 
length steel angle fittings. The five mounts at 



Figure 7 - EDM Launcher Set-Up On Barge 


Results of shock isolation analyses de- 
scribed in the next section under D i^M Analy- 
ses, verify lower shock transmis ,r with the 
plate compared to EDM platform at all locations 
particularly on the missiles and launcher ex- 
tremities. Since Mil-S-901C shock tests were 
successfully passed on the EDM platform, the 
plate platform was therefore shown to be accept 
able from an underwater shock standpoint, 

DDAM ANALYSES 

The Dynamic Design Analysis Method 
(DDAM) is specified by the Navy in Nav 
Ships 250-423-30 for analytical verification of 
underwater shock integrity of shipboard equip- 
ment. The DDAM analytical procedure was 
therefore used for both shock and vibration 
response analyses during development of the 
NSSMS launcher as well as for the shock isola- 
tion platform. 

The method was initially programmed on a 
UNIVAC 1 108 computer for launcher natural 
frequency and mode shape analyses and for 
shock response using a separate program for- 
mat. Both of these programs were converted 
to a CDC 6700 computer for analyses of the 
shock isolation platform. 
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The dynamic models used in the DDAM 
analyses consisted of a sixteen degree-of- 
freedom vertical model, and a 25 element 
transverse/rocking model. The vertical 
model is shown in Figure 8 as an example. 

The 25 elements in the transverse model in- 
clude 9 transverse, 9 rotational and 7 missile/ 
cell vertical elements. The eight missiles 
were separate elements in the vertical model 
but similar port and starboard missiles were 
combined into four elements in the transverse/ 
rocking model. 
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Figure 8 - Vertical Shock Model 

Primary concentration in the DDAM 
analyses was in the vertical direction. Since 
the plate platform is much softer in trans - 
verse/rocking directions than the EDM plat- 
form, more favorable shock isolation character- 
istics are predictable in the transverse direc- 
tion. Calculations also indicated response 
accelerations well below the 6g minimum trans- 
verse unidirectional load allowed by the DDAM 
procedure. 

Influence coefficients were calculated for 
the launcher and missiles during EDM launcher 
development. A structural compute • program 
designated STIFF was used to calculate cell 
structural deflections and was found very 
accurate without joint degradation factors. 

The missiles were treated separately as dis- 
tributed bodies on their supports and then 
lumped into single elements for equivalent 
loading. Elements such as the stand, pedestal, 
truss, gear system and bearings were analyzed 
separately and utilized appropriate joint reduc- 
tion factors. Influence coefficients were also 
updated based on experimental natural fre- 
quency and mode information obtained during 
EDM launcher shock and vibration tests. 

The accuracy of the dynamic model and 
influence coefficient mass inertia matrices in 
representing the stiffness and inertia of the 
launcher, was demonstrated by experimental 
frequencies, found during hard mounted launcher 


vibration tests, which agreed with analytically 
predicted frequencies within 0. 5 Hz. 

A total of 11 frequencies with correspond- 
ing modes were used in the DDAM analyses 
such that full launcher dynamic mass was repre- 
sented. In addition to the 79, 000 lbs/in. (13.8 
(10)'6 N/m) plate platform stiffness found dur- 
ing experimental tests, stiffnesses of 17.5(10)6 
N/m, 21. 9(10)6 N/m and 26. 3(10)6 N/m were 
used in the analyses for parametric purposes. 
The EDM platform which was compared analyti- 
cally by the DDAM method had a stiffness of 
147, 000 lbs/in. (25. 7(10)° N/m). Inputs speci- 
fied by NRL 1396 of Reference 2, were used in 
the DDAM analyses. 

Figure 9 shows a comparison of predicted 
response at the launcher base (element 16 in 
Figure 9) and the inboard and outboard missiles 
(elements 5, 6 and 2, 3 respectively) plotted 
versus the first mode frequency of the launcher 
system on the EDM barge test platform and on 
the plate platform. The results indicate that 
the plate platform provides better attenuation of 
underwater shock than th< EDM platform par- 
ticularly on the missiles and launch cells. The 
dependence of shock isolation on the fundamen- 
tal frequency is inherent in this data, and was 
further supported by the parametric studies 
with various platform stiffnesses. 


PLATE EDM 



7.0 7.5 8.0 8.5 9.0 

FREQUENCY (Hx) 


Figure 9 - DDAM Response 
DDAM LOADS 

Resulting design loads at each element ob- 
tained from the DDAM analyses, on both the 
EDM and plate platform, are listed in Table 1. 
Loads obtained from the EDM platform barge 
tests, described in the next section, which were 
used for structural analyses, are also listed in 
Table 1. 
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Table 1 


DDAM versus Barge Test Vertical Loads 


ELEMENT 

WEIGHT 

EDM 

DDAM 

PLATE EDM 
DDAM BARGE 

KG 

G 


G 

1 

152 

10.5 

9.2 

8.0 

2 

282 

14.2 

11.5 

10. 1 

3 

282 

14.2 

11.5 

10.1 

4 

229 

7.7 

7.2 

4.4 

5 

282 

12.5 

11.4 

9.3 

6 

282 

12.7 

11.5 

9.3 

7 

438 

7.2 

7.2 

7.8 

8 

152 

10.5 

9.2 

8.0 

9 

282 

14.2 

11.5 

10.1 

10 

282 

14.2 

11.5 

10.1 

11 

229 

7.7 

7.2 

4.4 

12 

282 

12.5 

11.4 

9.3 

13 

282 

12.7 

11.5 

9.3 

14 

438 

7.2 

7.2 

7.8 

15 

1927 

7.3 

7.3 

7.8 

16 

535 

7.1 

7.0 

8.9 

PLATFORM LOAD = 

EWG 

9.6 

8.8 

8.3 


W 

The comparative DDAM loads in Table 1 
show about nine percent lower shock loads on 
the plate than the EDM platform, substantiating 
the superiority of the plate platform from an 
underwater shock standpoint. The DDAM 
structural design loads also show good agree- 
ment with structural loads obtained from the 
EDM barge tests. The EDM barge loads 
averaged about 86 percent of the DDAM loads 
on the EDM platform and loads showed con- 
sistent trends between elements. 

A total down vertical load of 8. 8g referred 
to the overall 14,000 lbs launcher weight was 
obtained from the DDAM analyses. The load 
was obtained from the weight data in Table 1, 
by multiplying each element weight by the 
design C load and then dividing by the 14, 000 
lbs (6356 kg) total weight. 

Resulting overall loads from the DDAM 
analyses are 8. 8g vertical and 6g transverse. 
The 6g transverse load is the least allowed by 
DDAM in the transverse direction and was used 
since preliminary checks indicated lower loads 
from the analytical procedure. 

The DDAM procedure does not require 
summation of vertical and transverse loads 
since as noted previously, peak shock re- 
sponses in the two directions at various loca- 
tions do not occur simultaneously. A conserva- 
tive procedure with respect to the 8, 8g vertical 
and 6g transverse loads is to determine stress 
for the 8. 8g unidirectional peak as well as for 
the 6g peak and then use an SRSS resultant 
stress such that: 


= V (a 8.8g )2 T ( *6g )2 


The peak stress using this approach was close 
to the stress using combined EDM platform 
barge test loads discussed in the next section. 

BARGE LOADS 

Underwater shock loads on the EDM 
launcher and missiles during Mil-S-901C 
barge tests were used for structural analyses 
of the plate platform. This was considered 
conservative since it was shown by DDAM 
analyses that about 9 percent lower shock loads 
would occur on the plate platform. 

The Mil-S-901C test requires that the 
launcher with 8 missiles be mounted on the 28 
x 16 ft (8. 5 x 4. 9 m) floating shock platform as 
shown on Figure 7, and then be subjected to a 
27. 2 kg explosive charge suspended at a water 
depth of 7. 32 m for each of five shots. Hori- 
zontal distances to the charge from the side of 
the barge along its transverse centerline, are 
18. 3, 12.2, 9. 1, 7.6 and 6. 1 meters respec- 
tively for each of the five shots. 

Shock inputs are progressively more 
severe at the closer distances with the closest 
6. 1 meter standoff being most severe. This is 
shown on Figure 10 by vertical velocity histo- 
ries measured on the barge for the five shots. 
Velocities at the launcher base were sub- 
stantially less than on the barge, varying from 
26 percent at the 18. 3 meter standoff to 57 per- 
cent at the closest 6. 1 meter standoff. 

Acceleration data were measured versus 
time through the underwater blast response 
period at 17 locations including the launcher 
base, the upper pedestal structure, extremities 
of cell structure both toward and away from the 
blast and both inboard and outboard missiles. 

Since shock response peaks occur at dif- 
ferent times at various locations in the vertical 
direction and at still different times in the 
transverse direction, detailed analyses of the 
shock response data are required to determine 
maximum simultaneous loads on the platform. 

Acceleration time histories measured on 
the launcher and missiles in the vertical direc- 
tion for the most severe 6. 1 meter standoff 
charge, are plotted on Figure 1 1. Similar 
acceleration time histories in the transverse 
direction are plotted on Figure 12. It is evident 
from these data that both vertical and trans- 
verse peak accelerations are reached at the 
launcher base and pedestal first, and that peaks 
on the cell structure and on the missiles occur 
considerably later. Vertical and transverse 
peaks at the same locations also generally 
occur at different times. 

The approach used to determine the maxi- 
mum simultanecus vertical plus transverse 
load on the EDM platform, from the data on 
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Figures 11 and 12, was to determine the overall 
load acting on the platform at each instant 
through the blast period. The instantaneous g 
load at each location was multiplied by the 
effective weight at that location for all launcher 
missile weight elements. The sum of the ele- 
ment weights, times their accelerations was 
then divided by the total weight W to obtain the 
overall load acting on the platform based on the 
launcher missile weight. The element weights 
used in this analysis are the same as listed in 
Table 1 for the DDAM analyses. 

G n 

— n n 

<- = w 

Results of the barge test loads analyses 
are shown on Figure 13, in terms of the over- 
all launcher vertical and transverse g loads 
with time, during the peak shock load period. 
The maximum load combination is a peak 
vertical shock load of 4. 5g which acts simul- 
taneously with a transverse load of 4. 6g. Add- 
ing the Ig static weight acting down to the 4. 5g 
dynamic load gives the following simultaneous 
underwater shock loads on the EDM platform: 

Vertical 5. 5g 
Transverse 4. 6g 


Since the plate platform would undergo 
lower shock loads than the EDM platform, based 
on DDAM analyses, the above EDM experimen- 
tal loads were considered conservative for 
structural design of the plate platform. 

STRUCTURAL ANALYSIS 

The ASSAS computer model of Figure 14 
was used for early parametric studies to opti- 
mize the platform plate and for structural 
analyses of the final configuration. Elements 
were taken in six layers throughout the plate as 
well as in the launcher base and outer mounting 
ring. Smaller elements were taken in the 
region of plate, launcher and ring interfaces. 

The effect of gaskets, such as between the plate 
and launcher base, was considered by softening 
elements immediately adjacent to the plate. 

Deflection analyses indicated that the 
tapered plate, 7. 62 cm at the center and 2. 54 cm 
at the outer periphery provided the soft vertical 
and sufficiently high moment stiffnesses desired. 
A sketch of plate deflection patterns under down 
vertical and rocking moment loads is shown in 
Figure 15. 

Stress analyses of the platform were con- 
ducted under the 5. 5g vertical and 4. 6g trans- 
verse rocking loads developed in the previous 
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Figure 11 - Barge Test Vertical g' s 

section under barge loads. Calculated radial 
and tangential str' ss distributions versus plate 
radius, are shown in Figure 16, due to the 
axial and moment loads. It is seen that tangen- 
tial stresses are highest, for both type loads 
and that peak stresses occur in the region of 
the launcher bolt circle. 

The computer model considered axisym- 
metric geometry, and did not specifically take 
the plate bolt holes or localized introduction of 
loads through the bolts, into consideration. 

Strain gauge data, measured during static te3ts 
in critical areas of the launcher bolt circle, 
were therefore used for plate structural verifi- 
cation under maximum underwater shock loads. 
Experimental stresses were less than calculated 
at the bolt centerline but were higher just for- 
ward and aft of the bolts due to load transfer re- 
sulting from the bolt clamping action. 

In addition to the plate, structural analyses 
were made on other principal platform elements 
such as the outer steel angle ring, the launcher 
base, 2. 54 cm launcher hold down bolts, l, 27 
cm shoulder bolts and damper attachments. 
Stiffness and bearing load analyses of the elasto- 
meric seals and bushings were also conducted 
and found satisfactory. Analyses were also 



Figure 12 - Barge Test Transverse g' s 

made to ■ erify that peak loads due to other 
environments such as green water, gun blast, 
vibration and sea state conditions did not ex- 
ceed underwater shock loads. 

With the peak 5. 5g vertical and 4. 6g trans- 
verse underwater shock loads, considering 
elastic material behavior, a peak combined 
stress of 306, lOOPa at the 14. 5 inch (35. 8 cm) 
radius was obtained from the platform static 
test strain gauge data. This point, noted on 
Figure 17, is below the 317, 200Pa elastic 
allowable and well below the 427, 500Pa modulus 
of rupture allowable considering some plastic 
deformatio- of the 6061 T6 material. 



Figure 13 - Barge Test Load Summary 
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Figure 14 - AASAS Computer Model 



Figure 15 - Deflections Due Down and 
Transverse Loads 
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Figure 16 - Combined Load - Moment and Axial 

It is evident from Figure 17 that stresses 
due to the unidirectional 8. 8g vertical and 6g 
transverse DDAM loads would be below the 
317, EOOPa allowable. The SRSS resultant of 
308, 900Pa using the DDAM 8. 8g and 6g stresses 
is close to the stress using the combined barge 
loads. 

The 317, 200Pa plate allowable of Figure 17 
is based on material certification data obtained 
from Aluminum Company of America for the 



Figure 17 - Maximum Stress versus Allowable 

6061 T6 plate material used in the platform. 

The certification data indicates minimum 
stress allowables of (7 U =332, 300Pa ultimate, 
cry = 302, OOOPa yield and 13 percent to 15 per- 
cent elongation. For shock applications in- 
volving short period dynamic loads, a stress 
allowable c r, between yield and ultimate is 
allowed as follows. The factor F = 0. 5 for 
typical plate cross sections. 

<J A = a y + F (<r u - cr y ) 

From a practical standpoint, the higher 
allowable of cr^ = 427, 500Pa could be used for 
the plate configuration. This higher allowable 
is based on the modulus of rupture approach 
wherein major unyielded areas of the plate pre- 
vent significant permanent deformations due to 
localized yielding in small surface areas. 
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